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ELECTROSTATIC SOUND WAVE MODES IN A PLASMA* 


F. W. Crawford 
Microwave Laboratory, W. W. Hansen Laboratories of Physics, Stanford University, Stanford, California 
(Received April 17, 1961; revised manuscript received May 31, 1961) 


The dispersion relation for ion waves in an 
infinite plasma is!»* 
2 
wo... 
we = 1+ = RT ) (1) 
pi p e 





where my, is the positive-ion mass, £§ is the phase 
constant, T, is the electron temperature, zero 
ion temperature is assumed, and is Boltz- 
mann’s constant. The constant y depends on the 
assumptions made concerning collisions and will 
lie between 1 and 3.? The ion plasma frequency 
is given by 

Wy; = (ne? /eqm 4)”, (2) 
where Np is the ion density, e is the electronic 
charge, and €, is the permittivity of free space. 

In the long-wavelength limit w/8 + (yk Te/my)”, 
giving rise to the so-called electrostatic sound 
wave mode of propagation’ with velocity inde- 
pendent of n,. A low-frequency cutoff dependent 
on tube dimensions might be expected. 

Ion waves have often been suggested as a source 
of low-frequency fluctuations occurring in the 
parameters of de discharges.* To test this hy- 
pothesis, dc mercury-vapor discharges at a 
pressure of ~1yu have been studied in tubes of 
the type shown in Fig. 1. Each contained an 
oxide-coated spiral cathode, mounted transverse- 
ly in a molybdenum cup. One tube of uniform 
diameter was available, and eight others of 
various column diameters, but roughly similar 
construction in the enlarged section containing 


the cathode assembly. The dimensions were 
comparable to the electron=-neutral mean free 
path. Spectra of anode voltage fluctuations were 
observed to contain peaks varying slowly with 
current, and typical results are shown in Fig. 1. 
Similar spectra were picked up on a ring slid- 
ing along the outside of the glass. Attempts to 
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FIG. 1. Experimental discharge tube, and typical 
spectra of low-frequency fluctuations. 
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measure phase shift have indicated that the 
fluctuations must travel along the column at 
velocities > 10° cm/sec over most of its length, 
but that over the first few centimeters the phase 
velocity is of the order of 10’ cm/sec.‘ 

Although it is doubtful whether longitudinal ion 
wave modes have ever been observed, it is pos- 
sible that radial modes, generated near the 
cathode, and near cutoff within the column, could 
explain these results. Their theory may be de- 
rived from the cylindrical-polar form of the wave 
equation, and if such solutions are valid the fre- 
quencies of higher order modes should bear 
fixed numerical relations to the frequency of the 
lowest mode. For a mode in which the potential 
varies as cosn@, and is assumed to have a node 
at the wall, 8 will be determined by the zeros of 


Jn(pd/2). 

Data on the relations between the frequencies 
of peaks observed in Tubes I-IX are shown on 
Fig. 2(a) over a wide range of current and num- 
ber density (~ 10°-10"° ions/cc). Most of the 
frequency ratios are distributed about the aver- 
age values 1.68 and 2.34 which are close to the 
ratios suggested by our simple theoretical con- 
siderations. The occasional second harmonics 
observed could easily arise through nonlineari- 
ties. 

Assuming that radial electrostatic sound waves 
are involved, we may substitute for the constants 
in the expression for the frequency, 


f= (2.405/nd)ykT,/m .)”, 










































w&|2 
-~* 
BESSEL 
- J, ~2.405 (f.),5. f 
«x FUNCTION { ° of) +S. SEONt,) f =1.59 Soe 30 
rs ZEROS J, ~3.832 (f,) — —_ 
Ww 
e 
ra 35! 2.34 
4 
“ 2r e . & *—*- 2.00 
a ° ° 
—e—_—e S, _ & gy 980 Fe 
a a ieee % 1.68 
o ee 
(a) 
” 
o I 1 
e 
< 
a 
10° 
s CONSTRICTED 
a 
4 FIG, 2. (a) Frequency ratios 
§ of peaks in fluctuation spectra. 
10° af ‘ (b) (fod) variation with column 
. Te= 35000 °K current density. 
2 Te=!7800 °K y=2.0 
Y =1.7 TUBE d cm 
UNIFORM o°t-8 1.01 - 1.07 
TUBE 4-M-WI}.95-2.00 
*= Ml 2.18 
(b) v= wi 2.56 
10° oom 4.81 
1.0 10 100 1000 


(ma/cm?) 


CURRENT DENSITY 











VoLUME 6, NUMBER 12 


PHYSICAL REVIEW LETTERS 


June 15, 1961 





to obtain the result 
Sod = 490(/T_)”? cm/sec, (3) 


where f, is in cps, d is incm, and Tg is in °K. 
Figure 2(b) shows (f,d) plotted against column 
current density. It has been assumed that the 
plasma column fills the tube and no correction 
has been attempted for sheath thickness. Lang- 
muir probe measurements of Tg were made at 
the points marked. Substitution of these values 
in Eq. (3) gives values of y in the range expected. 
The temperatures were measured nearer to the 
anode than the cathode so that the values may not 
be appropriate to the actual generation region. 
The results suggest that an electrostatic sound 
wave mechanism may be operating, and that it 
is enhanced by the presence of a constriction. 
The Bessel function relations are appropriate to 
the first and second radially varying modes, and 
the first azimuthally varying mode, though it is 
difficult to see why an azimuthally varying mode 


should be excited so strongly and detected so 
easily. To determine the ratios more accurately, 
modification of Eq. (1) is required to include the 
axial and radial drift velocities of the ions and 
electrons. 

It is a pleasure to acknowledge many discus- 
sions of the problems involved in this work with 
Dr. G. S. Kino, and private communications from 
Dr. R. L. Moore and Dr. I. Alexeff on the radial 
ion wave mechanism. 
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ENHANCED DIFFUSION AND OSCILLATIONS IN WEAKLY IONIZED PLASMAS 


Jean-Frangois Bonnal, Georges Briffod, and Claude Manus 
Section d’ Ionique Générale, Service de Physique Appliquée, Centre d’ Etudes Nucléaires de Saclay, 
Gif-sur-Yvette (Seine-et-Oise), France 
(Received May 8, 1961) 





The present study deals with observations on 
enhanced diffusion performed on a weakly ionized 
plasma. 


Experiments have been performed with a Reflex- | 3 


type discharge in H,, N,, and Ar, the pressure, 
density, and magnetic field conditions of which 
are variable. The discharge works in continuous 
mode. The diffusion is measured by means of 
the escape flux of the ions transversally to the 
magnetic field B,. The measurement is made 
with a probe external to the plasma, negatively 
polarized with regard to the plasma potential. 
The variation of the escape flux as a function of 
the magnetic field for different values of the pres- 
sure parameter is shown in Figs. 1 and 2. 

All curves show three separate regions: 

1. B< By: The escape flux decreases with the 
magnetic field and increases with the pressure 
according to the classical laws of diffusion con- 
firmed experimentally by several researchers!~* 
(Fig. 1). 

2. Bo< B< By: Anew diffusion mode appears 
in which the escape flux increases with the mag- 
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FIG. 1. Escape flux of ions transversally to the 
magnetic field; 0<B <1000 gauss. These results have 
been obtained with a cold-cathode P.I.G. discharge, 
working with hydrogen, of length 110 mm, diameter 
30 mm, voltage between anode and cathode 400 volts, 
current 200 ma, 
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FIG. 2. Escape flux of ions transversally to the 


magnetic field; 1000 < B <4000 gauss. 


netic field. This has to be connected with the 
Lehnert experiments. The influence of the pres- 
sure on the position of the curves is just opposite 
to that of the previous mode: The escape flux 
grows as the pressure decreases. (This mode 

is accompanied by a considerable noise emission.) 

On each curve, the maximum value of the flux 
is reached for B= By. As the pressure de- 
creases, both this maximum and the critical 
magnetic field are shifted toward lower field 
values. 

3. B> By: The escape flux decreases again 
along with the noise level. But the order of the 
curves with respect to the pressure is maintained 
(Fig. 2). 

Regarding the dependence of the critical field 
on the discharge radius and the atomic mass of 
the gas (H,, N,, Ar), the laws are similar to 
those pointed out by Lehnert and Hoh.*~® 

Finally, the critical magnetic field is inde- 
pendent of the density of the charged particles. 

The intense emission of noise by the discharge 
above the critical field is detected by an external 
antenna for the high-frequency component and is 
observed on the oscillogram of the anode-cathode 
voltage for the low-frequency components. The 
high-frequency spectrum is particularly wide 
and reaches values higher than the plasma fre- 
quency. 
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FIG. 3. High-frequency noise level. 


The relative variation of the noise level with 
the magnetic field at 1000 Mc/sec is shown in 
Fig. 3. The general trend of these curves and 
their dependence on the pressure are similar 
to those of Fig. 2. 

To sum up: (1) The noise appears only in the 
case of B> B,; (2) its amplitude increases as 
the pressure decreases; and (3) its maximum 
value shifts toward weak magnetic field as the 
pressure decreases. 

The amplitude of the low-frequency noise 
spectrum follows roughly a variation identical to 
the variation shown in Fig. 3. 

The following remarks can be made: 

(a) Even though the length of the discharge and 
the pressure range are entirely different from 
the conditions corresponding to the experiments 
carried out by Lehnert, there appears a critical 
magnetic field and an enhanced escape flux.’ 

(b) The high-frequency noise spectrum is very 
important and is connected with the value of the 
escape flux. These different observations, 
mainly the existence in the spectrum of high- 
frequency noise of components corresponding to 
frequency higher than the plasma frequency, has 
suggested a theoretical treatment based on a 
microscopic approach. The instabilities are 
related to the existence of anisotropies in the 
velocity distribution function, anisotropies due 
mainly to the existence of an electrical field in- 
side the plasma. This treatment is quite differ- 
ent from that developed by Hoh® and Kadomtsev 
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and Nedospasov,° and will be published in the 
near future. 
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OBSERVATIONS OF de HAAS—van ALPHEN OSCILLATIONS IN p-TYPE PbTet 


P, J. Stiles,” E. Burstein, and D. N. Langenberg 
University of Pennsylvania, Philadelphia, Pennsylvania 
(Received May 3, 1961) 


We should like to report the results of a study 
of de Haas—van Alphen (dH-vA) oscillations in 
the magnetic susceptibility of p-type PbTe by 
pulsed magnetic field techniques. Although ex- 
tensive studies of the dH-vA effect in metals 
have been made, as far as we know this is the 
first time that susceptibility oscillations have 
been observed in a semiconductor. Oscillations 
in transport properties such as magnetoresistance 
and Hall effect have, of course, been observed 
in semiconductors as well as in metals. Because 
studies of magnetoresistance and Hall effect oscil- 
lations require electrical contacts and involve the 
current as an additional parameter, the meas- 
urement of susceptibility oscillations by the im- 
pulse method’ has the advantage of greater sim- 
plicity. We feel that such measurements will have 
wide application in the study of the band structure 
of semiconductors. 

In each of the four p-type PbTe samples studied, 
more than 30 oscillations were observed. Three 
of the PbTe samples were cut from a single boule 
as close to each other as possible. Hall effect 
measurements made by Allgaier” yielded a hole 
density of 3 x 10'*®/cm® for this boule. The fourth 
PbTe sample had a hole density of about 1 x 10'*/ 
cm*. These samples were circular cylinders 
with the cylinder axis parallel to one of the three 
principal symmetry directions. The magnetic 
field was applied along the axis of the cylinder. 

A typical photograph of the oscillations in dm /at, 
where m is the total magnetic moment of the 
sample, is shown in Fig. 1 and the successive 
maxima and minima are plotted versus B™ in 





FIG. 1. dH-vA oscillations in p-type PbTe. The 
double trace shows the oscillations for increasing 
and decreasing magnetic field. The magnetic field 
varies from 7.3 kgauss at the left to 16.4 kgauss at 
the right. 


Fig. 2 for two different orientations. Table I 
lists the periods measured in the three symmetry 
directions. 

It has been suggested by earlier workers that 
the energy surfaces at the edge of the valence 
band in PbTe are ellipsoids of revolution ori- 
ented along the (111) directions.’ Allgaier has 
recently suggested on the basis of his magneto- 
resistance and Hall effect studies at 77°K and 
room temperature that an additional valence 
band maximum appears to be involved. Our data 
indicate a band maximum at k=0 and band maxi- 
ma at the centers of the {111} faces of the Bril- 
louin zone, although it is not clear that the ad- 
ditional valence band maximum we see is the 
one suggested by Allgaier. The extremal cross 
sections associated with the k=0 maximum vary 
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FIG. 2. Successive maxima and minima in dm/dt are plotted as integers and half integers vs B-', with 
n=3.910'® holes/cm? and T=4.2°K, (a) The plot for the longer period observed with B|| {100). (b) The plots 
for all three periods observed with B]| (111). Twenty-four more maxima and minima observed in the series 


with the smallest slope are not shown here. 


with orientation by about 25%. Assuming para- 
bolic bands, the energy surfaces at the {111} 
zone faces are ellipsoids of revolution with a 
longitudinal to transverse mass ratio of 6.4+0.3. 
From the dependence of the oscillation ampli- 
tude on magnetic field and on temperature, we 
derive a transverse effective mass for the el- 
lipsoids of (0.042 + 0.006)m, (m,=free electron 
mass) and an equivalent broadening temperature, 
T,, of about 7°K to 10°K. T, includes the appar- 
ent broadening due to hole density variations in 
the sample as well as collision broadening. We 
feel that the large value of Ty, is primarily due 
to variations in hole density. On this basis the 
value of T, observed in these experiments cor- 
responds to a 10% variation in hole density, a 
reasonable value for such materials. 

If the complete geometry of the Fermi surface 
is known, one can determine the total number of 
holes from the total volume enclosed by the sur- 


fade in k space. We estimate the number of holes 


in the three crystals cut from a single boule to 
be 3.9x10"*/cm', assuming four ellipsoids and 
a sphere at K=0 of cross section equal to the 
average of the three measured cross sections. 
This is reasonably close to the density of 

3 x 10?*/cm* obtained from Hall effect measure- 
ments. At this concentration the Fermi level 
is 0.035+ 0.005 ev below the top of the ellip- 
soidal band. By comparing the results from 
the two samples with different hole densities, 
and assuming parabolic bands and our value of 
the ellipsoid transverse effective mass, we ob- 
tain an estimate of the effective mass of the 
central section, m* =(0.12+0.02)m,. We also 
estimate that the band edge at K=0 is higher in 
electron energy than the band edge at the Bril- 
louin zone face by 0.002+ 0.002 ev. 

It is of interest to consider the factors deter- 
mining the sensitivity of the impulse method for 
observing the dH-vA effect in semiconductors 
and to compare this sensitivity with that achieved 


Table I, Experimental dH-vA periods. 





Orientation of the 
magnetic field 


Period* 
(gauss~! x108) 


Measured hole 
density (cm-) 





Ellipsoids k=0 section 
400) 13.3 6.07 9.5 x10'? 
4.00) 5.37 2.63 
410) 7.50 2.42 
3.70 3.9 x10'® 
4111) 8.39 2.15 
4.18 





a 
There is an error of 2% associated with each period, primarily due to the calibration of the magnetic field. 


Calculated from dH-vA periods. 
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in observations of the effect in metals. The ex- 
perimentally observed signal is proportional to 
dm /dt = V(dm/dH )(dH /dt), where V is the sample 
volume and M is the magnetization. In semi- 
conductors the carrier density is orders of mag- 
nitude less than it is in metals, the oscillation 
period is much longer, and dM/dH may be 
smaller by a factor of 10°. On the other hand, 
because of the longer period in semiconductors, 
larger field inhomogeneities due to eddy current 
shielding can be tolerated. Furthermore, the 
eddy currents are reduced by higher resistivi- 
ties. Both V and dH/dt may therefore be larger 
for a semiconductor than for a metal and the 
smaller sensitivity due to the small dM/dH can 
be partially compensated. The signals observed 
in the present experiments are 10-100 times 
smaller than typical signals observed in experi- 
ments on metals. Nevertheless the sensitivity 

is quite adequate. In fact, the sensitivity of the 
method is such that the only limitation to observ- 
ing the effect in semiconductors is the theoreti- 
cal one that the materials must have E,- > k(T+ T,), 
where Ey is the Fermi energy measured rela- 
tive to the edge of the band in question. Thus, 

in addition to the requirement that the sample 

be degenerate (E> kT), one needs uniform car- 
rier densities and long collision times (E> kT}). 
The requirements of a long collision time 

Ey > h/t) can be expressed in terms of the 





f 
i 
i 
4 








density and mobility of the carrier as n™°y > 10%5 
volt sec for the case of a single spherical en- 
ergy band. There are a number of semicon- 
ductors which readily satisfy these conditions, 
such as the lead salt series PbS, PbSe, and 
PbTe,* which we are now investigating. 

The authors wish to express their thanks to 
R. S. Allgaier, W. W. Scanlon, and B. B. Hous- 
ton for many enlightening discussions and for 
supplying the crystals investigated. We would 
also like to thank L. Indyk and B. N. Taylor for 
their valuable assistance in setting up the ex- 
periment. We also wish to express our ap- 
preciation to Dr. D. N. Stevens for his interest 
and encouragement in supporting this work. 
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GROUND STATE OF AN ISING FACE-CENTERED CUBIC LATTICE 


A. Danielian 
Wheatstone Physics Laboratory, King’s College, London, England 
(Received May 4, 1961) 


Many antiferromagnetic compounds have a 
face-centered cubic (fcc) magnetic structure. 
Although in the case of nearest neighbor inter- 
actions the ground state of such systems is 
clearly degenerate, no exact information has 
been available regarding either the magnitude 
of its degeneracy or any of its thermodynamic 
properties. In the case of an antiferromagnetic 
Ising triangular lattice, it has been shown that 
the entropy per spin is finite’ and the zero-field 
susceptibility infinite? at T=0°K. In this Letter 
we show that, assuming nearest neighbor Ising 
interactions, simple topological and energy argu- 
ments lead to a complete classification of the 
ground-state configurations. Hence, a precise 
enumeration of the states is possible and exact 
thermodynamic results are deduced at T =0°K. 

We consider an antiferromagnetic system 
consisting of N spin moments on an fcc lattice. 
Each spin is capable of existing in one of two 
possible states denoted +, -, and interacts with 
its nearest neighbors (n.n.) only. The interac- 
tion energy is +J for a pair of like spins (++ or --) 
and -J for a pair of unlike spins (+- or -+). 

First, the energy of the ground state of the 
system is determined. An fcc lattice of N sites 
can be subdivided into N tetrahedra, such that 
each n.n. interaction bond is found on one and 
only one tetrahedron; hence there are a total of 
6N such bonds. In the ground state, each one 
of these tetrahedra must be in its own ground 
state, i.e., two of its spins in the state + and 
two in the state -. This ground state of the tetra- 
hedron is denoted by the symbol w, and the cor- 
responding energy is -2J. Hence the energy of 
the ground state of the fcc system is -2NJ. It 
then also follows that every tetrahedron in the 
lattice is in state w. 

We next consider the fcc lattice built up of 
layers of triangular lattices. In Fig. 1 sucha 
layer (1) is shown-—the circles at the centers 
of the triangles denote sites on layer (2) which 
is immediately under layer (1), and points at 
the centers of the remaining triangles denote 
sites on layer (3) which is immediately above 
layer (1). Each triangle of (1) is the base of a 
tetrahedron, the vertex of which is a site either 
on (2) or (3). Also, each site on a given layer is 
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the vertex of a tetrahedron, the base of which is 
a triangle on an adjacent layer. 

Since in the ground state of the fcc system 
each tetrahedron is in the state w, it follows 
that each triangle of spins in a given layer 
should be in one of the states (++-) or (--+). 
Therefore any triangular layer should also be 
in one of its ground states. The number of such 
states is of the order 2°*”, where n is the num- 
ber of spins (or sites) on a triangular lattice.! 

Thirdly, we note that if the configuration of 
the ground state of layer (1) is determined, then 
the configurations of both layers (2) and (3) are 
also determined. This obtains because each tri- 
angle has one of the configurations (++-) or (--+), 
and since it is also required that each tetrahedron 
be in the state w, the state of the spin at the ver- 
tex of the tetrahedron is determined, and this 
vertex is in the adjacent layer. It follows that 
the ground configurational state of the whole fcc 
system is determined by the ground configura- 
tional state of any one of its triangular layers; 
hence we can deduce that its degeneracy is 
< 2°47 = »~N*5, We now proceed to show that 
the strict inequality holds. 
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FIG. 1. Three successive layers of a face-centered 
cubic lattice: layer (1) denoted by triangles; sites on 
layer (2) denoted by circles; and sites on layer (3) are 
at the centers of the remaining triangles. 
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FIG. 2. The a* and f* clusters. Notation as for 
Fig. 1. 


Figure 2 shows two clusters, each having a 
central spin with its 12 nearest neighbors. Such 
clusters contain 8 tetrahedra: 6 having the tri- 
angles of (1) as bases, the remaining two having 
the central spin for their vertices and their bases 
on layers (2) and (3), respectively. The condition 
that each tetrahedron must be in the state w re- 
duces the number of basic types of such clusters 
to four. Two of them, a* and g*, are shown in 
Fig. 2; the other two, a and g”, are obtained 
simply by reversing all the spins. 

It follows that in the ground state, each + 
spin is at the center of an a* or g* cluster and 


each - spin is at the center of an a or 8° cluster. 


The possible ground states of a triangular layer 
can therefore be obtained by grouping these four 
clusters. The number of possibilities is re- 
stricted by a property of the a cluster. If a par- 
ticular site—chosen as origin-—is the center of 

an a cluster, then the adjacent clusters along 

the 00 axis must also be a clusters. Hence all 
clusters along the 00 axis must be a clusters. On 
the other hand, adjacent clusters alongaa or bb 
axes may be either a or 8. If one therefore 

looks at the triangular layer along the 00 direc- 
tion, one can have a row (parallel to 00) of a 
clusters followed by another row of either a 
clusters or 6 clusters. Thus the triangular layer 
will consist of rows, each consisting entirely of 

a or g clusters (e.g., Fig. 1). The number of 
arrangements is clearly 2F where R is the num- 
ber of rows. 

The degeneracy of the ground state of the fcc 
system is therefore exp(AN”* 1n2), and it follows 
that in the limit N= the entropy per spin tends 
to zero as T+ 0°K. 

In the above scheme all triangular layers have 
identical configurations and since all possible 
configurations have zero moment, it follows 
that the susceptibility is zero at T=0°K. 

I am grateful to Dr. M. F. Sykes for drawing 
my attention to this problem and for valuable 
discussions. 





1G. H. Wannier, Phys. Rev. 79, 357 (1950). 

2M. F. Sykes (to be published), based on the theory 
developed by M. E. Fisher, Phys. Rev. 113, 969 
(1959). 





SUPERCONDUCTIVITY AT HIGH MAGNETIC FIELDS AND CURRENT DENSITIES 
IN SOME Nb-Zr ALLOYS* 


T. G. Berlincourt, R. R. Hake, and D. H. Leslie 
Atomics International, Division of North American Aviation, Canoga Park, California 
(Received May 24, 1961) 


Zero electrical resistance at unusually high 
magnetic fields and current densities has recent- 
ly been reported in the compound Nb,Sn,',? and in 
the alloy Mo—25 at. % (atomic percent) Re.* The 
present results on Nb-Zr alloys were obtained 
during the course of a systematic investigation 
of resistive superconducting transitions in cer- 
tain Ti, Zr, Nb, Hf, and Ta-rich binary alloys 
in magnetic fields up to 30 kgauss. The data 
show that cold-worked Nb-rich Nb-Zr alloys 





display zero electrical resistance at current 
densities as high as 10° amp/cm? at 30 kgauss 
and at 4.2°K.*»> As far as we are aware, only 
the specially prepared Nb-clad “Nb,Sn” cores 
of Kunzler et al.’ have exhibited greater zero- 
resistance current-carrying capacities above 
10 kgauss. The present measurements also 
reveal some interesting high-field supercon- 
ducting effects in cold-rolled Nb-Zr alloys: 
(a) a marked dependence of the magnitude of 
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the maximum zero-resistance current upon the 
direction of the externally applied transverse 
magnetic field relative to the rolled surface of 
the specimen, and (b) an anomalous increase of 
the maximum Zero-resistance current with in- 
crease of magnetic field strength over certain 
field regions in certain alloy specimens.® 
Figures 1 and 2 show curves of critical cur- 
rent density, J ¢ VS magnetic field strength, 
H, for some Nb—25 at.% Zr and Nb—12 at.% Zr 
specimens for different temperatures and trans- 
verse field orientations. Preliminary data for 
Nb-38 at.% Zr and Nb-88 at.% Zr specimens 
have also been obtained and will be reported at a 
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FIG, 1. Critical current density, J, for restoration 
of the first observable resistance vs applied transverse 
magnetic field, H, for Nb-—25 at.% Zr alloy specimens 
1 to 9 at different temperatures and field orientations 
parallel (H ||) and perpendicular (H) to the rolled faces 
of the specimens. For specimen 1, which was ground 
to thickness, Hy,,, and Hyj, designate field directions 
such that J. is, respectively, a maximum and a min- 
mum at 30 kgauss and 4.2°K. For source material, 
fabrication method, and cross-sectional dimensions 
of the specimens see Table I. 
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later date. The source materials, method of 
fabrication, and cross-sectional dimensions 

of the specimens are indicated in Table I. For 
each alloy specimen the data were normally 

taken by slowly increasing the current through 

the sample at constant field strength until the 
least detectable resistive voltage, Vy was ob- 
served. For the data of Figs. 1 and 2, V,=0.25 
uv. Separate tests using higher sensitivity ap- 
paratus indicated that J, vs H curves are vir- 
tually the same for 0.005 pv< V,< 0.25 pv. The 
current required to produce V, divided by the 
cross-sectional area of the specimen was defined 
to be the “critical current density,” J,. (In cases 
where J, 22x10* amp/cm’, the onset of resistance 
was accompanied by the almost simultaneous burn- 
out of the specimen.) The specimens were copper- 
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FIG, 2. Critical current density, J c’ for restoration 
of the first observable resistance vs applied transverse 
magnetic field, H, for Nb—12 at.% Zr alloy specimens 
10 and 11 at different temperatures and field orientations 
parallel (H ||) and perpendicular (H) to the rolled face 
of the specimen. For specimen 10, which was ground 
to thickness, Hy, and Hy)j,, designate field directions 
such that J. is, respectively, a maximum and a mini- 
mum at 30 kgauss and 4.2°K. For source material, 
fabrication method, and cross-sectional dimensions of 
the specimens see Table I. The line L and points P, 
Q, R,S are explained in the text. 
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plated over 0.6 cm at their ends and were then 
soft-soldered between copper current lugs. The 
copper potential leads (spaced about 0.7 cm apart) 
were clamped against the specimens with beryl- 
lium-copper spring clips. 

Aside from the large high-field current densi- 
ties observed for some Nb-Zr specimens, the 
wide range of J, values (=2 x 10° - 10° amp/cm?) 
shown in Fig. 1 for specimens of the same nomi- 
nal alloy composition is surprising. Further in- 
vestigation will be required to isolate the impor- 
tant variables affecting the high-field supercon- 
ducting resistive characteristics. Micrographic 
examination of the various source materials (see 
Table I) under polarized light indicated that all 
the samples were predominantly bcc as might be 
expected from an inspection of the phase diagram, 
considering that the source material alloys were 
all arc-quenched. A small amount of second phase 
material (probably hcp) was apparent under the 
microscope in all the source materials. Compar- 
ison of the J. vs H curves for specimens fabri- 
cated from the same source material, but with 
different amounts of cold-work (see Figs. 1, 2, 


7 





and Table I) suggests that the more intensive the 
cold-working, the higher the value of J c Cold- 
working has previously been observed to increase 
the “critical fields’ and the high-field critical 
currents®»%s?° in pure Nb. 

The data of Figs. 1 and 2 also show that J, val- 
ues for cold-rolled specimens depend markedly 
upon the direction of the transverse magnetic 
field relative to the rolled face of the specimen. 
Polar plots of J c YS transverse field orientation 
for the cold-rolled specimens of Figs. 1 and 2 are 
roughly elliptical in shape with maxima and min- 
ima at field directions, respectively, parallel 
(A ||) and perpendicular (HL) to the rolled faces 
of the specimens, regardless of (a) the position 
in which the specimens were soldered into the 
current lugs, or (b) the cross-sectional geometry 
of the specimens. J. vs H data for these two ex- 
tremal field directions are shown for several cold- 
rolled specimens in Figs. 1 and 2. The ratio 
JH \l)/J_(H4) at 30 kgauss and 4.2°K varies be- 
tween 2 and 3 for the cold-rolled specimens. A 
less prominent orientation effect is also observed 
in specimens No. 1 (Fig. 1) and No. 10 (Fig. 2) 


Table I. Source, fabrication method, and cross-sectional dimensions of Nb-Zr alloy specimens. 








Atomic Cross-sectional 
Specimen Source percent b dimensions 
number material Zr Fabrication (cm) 
1 a.m.° button L1 25 ground to 0.019 cm 0.019 x0.022 
2 a.m. buttonl1 25 c.r.? from 0.100 to 0.025 cm 0.023 x0.025 
3 a.m. casting M2 25 c.r. from 0.033 to 0.015 cm 0.015 x0.023 
4 a.m. casting M2 25 c.r. from 0.091 to 0.030 cm 0.018 x0.030 
5 CBC wire® 25 None 0.025 diam 
6 CBC wire 25 None 0.025 diam 
7 CBC wire 25 e.r. to ~0.007-cm strip 0.0074 x0.066 
8 CBC wire 25 c.r. to ~0.007-cm strip 0.0069 x0 .037 
9 CBC wire 25 c.r. to ~0.007-cm strip 0.0069 x0.038 
10 a.m. buttonL2 12 ground to 0.025 cm 0.025 x0.025 
11 a.m. buttonL2 12 c.r. from 0.100 to 0.025 0.023 x0.025 
12 a.m. buttonL2 12 c.r. from 0.100 to 0.025 0.022 x0.023 





,Caleulated from weights of starting materials. 


The method of producing the desired strip thickness is given. Specimens of small width were sheared from the 


oStrips so produced. 
a.m. means arc-melted. 


c.r. means cold-rolled with no intermediate or final anneal. 
Wire produced by the Chase Brass and Copper Company. An arc-melted ingot of diameter 0.89 cm was cold- 


drawn to 0.025-cm diameter wire with no intermediate or final anneal. 
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which have not been cold-rolled. The orientation 
dependence may be associated with preferred 
crystalline orientation and/or the configuration 
of the dislocation network. The orientation de- 
pendence observed here, in fields where flux pen- 
etration of the specimen is doubtless nearly com- 
plete, is not to be confused with the well-known 
field orientation dependence of the bulk specimen 
superconducting transition, which is determined 
strictly by the specimen geometry. Indeed, for 
most of the present specimens the nearly square 
cross sections (see Table I) insure that the field 
sees nearly identical specimen geometry whether 
set parallel or perpendicular to the rolled face of 
the specimen. 

The J, vs H curves for Nb-—12 at.% Zr speci- 
mens at 4.2°K (Fig. 2) show anomalous peaks at 
fields of about 25 kgauss. At 1.2°K the J, values 
increase and the peak apparently moves to a higher 
field value. The “peak effect” can be observed in 
a rather striking way by increasing the magnetic 
field at a constant current such as to pass through 
the anomalous J.-H region, e.g., along line L of 
Fig. 2. As the field is increased a resistive volt- 
age first appears at point P, reaches a maximum 
of 125 pv at point Q, and then disappears at point 
R. The appearance and disappearance of resistive 
voltage is also observed upon lowering the field 
from point S. A peak very similar to that shown 
in Fig. 2, and occurring at the same field value, 
was observed in another Nb-—12 at.% Zr specimen 
(No. 12, Table I) in two different runs, in between 
which the specimen was resoldered to the current 
lugs. Anomalous positive slope in Je vs H curves 
has also been observed in cold-rolled Ti—29 at.% 
Ta™ and cold-rolled Nb.° An unusual decrease in 
the resistance with increasing magnetic field over 
narrow field regions has been previously observed 
in superconducting transitions of Nb® and of Ti- 
29 at.% V,"" and is probably related in origin to 
the “peak effect” in the J, vs H curves of Fig. 2. 

The observation of large high-field critical cur- 
rents in flexible Nb-—25 at.% Zr wire suggests the 





use of such wire for superconducting solenoids 
generating fields possibly well in excess of 30 
kgauss. Investigation of this prospect and further 
studies of high-field superconducting phenomena 
are in progress in this laboratory. 

We wish to thank D. M. Sellman and P. Q. Sauers 
for valuable assistance in the electrical measure- 
ments, H. Conrad for metallurgical advice, and 
C. G. Rhodes for micrographic examinations. We 
are also grateful to J. Port and A. Springmeyer 
of the Chase Brass and Copper Company for their 
expert fabrication of Nb-Zr wire, and to H. Reiss 
for his enthusiastic encouragement of this work. 





*This research was supported by the U. S. Atomic 
Energy Commission. 
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ASPHERICAL 3d ELECTRON DISTRIBUTION IN BODY-CENTERED CUBIC METALS 


Frank Stern 
U. S. Naval Ordnance Laboratory, White Oak, Silver Spring, Maryland 
(Received March 15, 1961) 


We have estimated the asphericity of the 3d 
electron atomic scattering factor for x-ray and 
polarized neutron diffraction in body-centered 
cubic metals on the basis of a simple energy- 
band model. The results agree with the neutron 
observations on iron,’ which indicate a departure 
of the 3d charge distribution from spherical 
symmetry. 

At a lattice site having cubic symmetry, atomic 
d levels are split into a doubly degenerate set 
called e g and a triply degenerate set called tg ge 
Electrons in these states have an aspherical 
charge distribution and an atomic scattering 
factor which is not a smooth function of the Bragg 
angle.* In metals the overlap of wave functions 
of neighboring atoms leads to energy bands which 
have no degeneracy in general except for wave 
vectors in directions of high symmetry. The 
charge density will still have cubic symmetry, 
and can be expanded in the form 


p)=Lp WY, ¢), (1) 


where y, 6, and @ are spherical polar coordin- 
ates and the Y7(@, ¢) are cubically symmetric 
combinations of spherical harmonics of order J, 
normalized to 47 on the unit sphere.* If the total 
charge density in the crystal can be written as 

a sum of charge densities centered on the in- 
dividual lattice points,® the atomic scattering 
factor will be 


f(u,B,y)= Xi. fo,(rh (urs rr“dr, (2) 


where yw, 8, andy are the spherical polar co- 
ordinates of the change in the wave vector upon 
diffraction, i.e., 1 =47sinép, g/A, and j;(ur) 
is the spherical Bessel function of order /. For 
cubically symmetric combinations of d states of 
pure e, symmetry, for example, the charge 
density has the form p(f) = p,(r)[1+ (9/21)“Y,@,y)], 
and Eq. (2) reduces to the corresponding expres- 
sion given by Weiss and Freeman.® 

Determination of the 3d charge density in a 
metal requires a summation of charge densities 
of all the occupied states in the 3d band. A rough 
estimate, but one which probably displays the 
essential features of the asphericity, is to use 
the charge density near the peak of the radial 
wave function, since this will give the dominant 





contribution to the atomic scattering factor. We 
have in that way estimated the asphericity for 
metallic iron, using the results of a cohesive 
energy calculation® in which energy values and 
wave functions were found for each of the five 
3d bands for 13 inequivalent points of high sym- 
metry in the Brillouin zone. The total volume of 
the Brillouin zone was divided among the 13 
points, and the weighted sum of the charge densi- 
ty peak heights of the occupied states was calcu- 
lated. The ferromagnetism of iron was taken in- 
to account by displacing the spin-up and spin- 
down bands with respect to each other (“up” 

here is the direction opposite to that of the mag- 
netization). If we assume that the charge density 
in iron is contributed by seven 3d electrons and 
one 4s electron per atom,® and that the observed 
moment of 2.2 Bohr magnetons per atom is en- 
tirely contributed by an excess of 3d electrons 
with spin up, we find that there are 4.6 3d elec- 
trons with spin up, and 2.4 3d electrons with 

spin down. Under these assumptions, the total 
charge density near the peak at 0.4 Bohr radii 

is found to have 44% e,, symmetry and 56% to g 
symmetry. A spherically symmetric charge 
distribution would be 40% eg. The net spin-up 
charge density near the peak is 52% e 2” which 
we will compare with the experimental value, 
50%.» § 

Two corrections, neither of which has a large 
effect on the results for iron, must be made to 
the calculated asphericity. Exchange polariza- 
tion contracts the wave functions of spin-up elec- 
trons relative to those of spin down. As an upper 
limit to this effect, we take the result of Wood 
and Pratt’ that the ratio of spin-up to spin-down 
charge density is 1.14 near the peak in the 3d 
radial wave function. The result is to change 
the fraction of the net spin-up peak having e g 
symmetry from 52% to 51%. The exchange 
polarization effect is small in all cases studied. 

A second correction is the explicit crystal 
field shift. Callaway and Edwards? find that the 
é, States are raised relative to the fg, states by 
0.094 Z-0.128Z ev, depending on the choice of 
wave functions, where Z is the net charge near 
each nucleus. We estimate an upper limit to Z to 
be 2.7, of which 1.0 comes from the uniformly 
distributed 4s electron, and 1.7 is an estimate of 
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the 3d charge density which is spread out.* Thus 
an upper limit to the amount by which the e 
states are raised is 0.35 ev, which is about 4% 
of the 3d bandwidth. If we take this shift into 
account for iron, the net spin-up charge density 
peak is 56% e g If both the exchange and crystal 
field corrections are made, the charge density 
peak is 54% e g 

For comparison with experiment we use the 
average of the corrected and uncorrected values 
described above. Thus the calculated net spin- 
up charge density peak in iron has 53% e g Sym- 
metry, and the calculated total charge density 
peak has 43% e g Symmetry. The former value 
is to be compared with the estimate® from polar- 
ized neutron diffraction measurements’ that the 
net spin-up charge distribution in iron is 50% e 
Pickart and Nathans® find in ordered Fe,Al that 
the net spin-up charge, for those iron atoms 
which have only iron atoms as nearest neighbors, 
is about 60% e g The x-ray form factor measure- 
ments of Batterman, Chipman, and DeMarco’® do 
not appear to allow an estimate of the asphericity 
of the total 3d electron charge distribution to be 
made. 

It is tempting to try to extend the calculation 
made for iron to other body-centered cubic met- 
als, and in particular to disordered ferromag- 
netic alloys." Our calculation can be so extended 
if the relative order of energy levels is deter- 
mined primarily by the bcc structure, and only 
secondarily by the particular transition elements 
occupying the lattice sites. Results based on 
this assumption are shown in Table I. The mini- 
mum asphericity of the net spin-up charge peak 
occurs near the alloy having the highest moment, 
while the maximum occurs for an atomic number 
slightly less than that of iron. 

The asphericities calculated using charge-den- 
sity peak heights differ only slightly from those 
that would be found from density-of-states consid- 
erations. That makes possible a simple descrip- 
tion of the results. In nonmagnetic material the 
lowest states of both spins are filled until all va- 
lence electrons are used up. This fixes the posi- 
tion of the nonmagnetic Fermi level. In magnetic 
material the asphericity of the net spin-up charge 
density is determined by the symmetry character 
of the states in the neighborhood of the nonmag- 
netic Fermi level. In iron these states have pre- 
dominantly e, symmetry, ® which accounts for 
our results. This qualitative method can be ex- 
tended to face-centered cubic metals. Recent 
band structure calculations for nickel and fcc 


676 





Table I. Calculated asphericity of the 3d charge 
density peak. The columns represent body-centered 
cubic transition metals and alloys.* The numbers of 
3d electrons per atom having spin up and spin down 
are given by N, and N_, respectively. The third row 
gives the percentage of the total charge density peak 
with €g symmetry. The fourth row gives the corre- 
sponding percentage for the net spin-up peak height. 
The first number of each pair is calculated directly 
from the band structure, > while the second number 
includes corrections for exchange polarization and 
crystal field effects. The best value to use in compar- 
ing with experiment is the average of the two percent- 


ages. 





N4-N. 1.5 1.8 2.1 2.2 2.5 2.2 
N++N. 6.3 6.6 6.9 7.0 7.5 7.8 
X ray 44,37 42,41 45,43 44,43 40, 39 39, 39 
Neutron 58, 42 49,48 55, 56 52,54 44,46 48, 47 





a 
See reference 11. 
See reference 6. 


iron” suggest that the nonmagnetic Fermi level 
corresponding to the ten valence electrons of 
nickel lies in a portion of the band having states 
of predominantly tg, symmetry, while the non- 
magnetic Fermi level for the nine valence elec- 
trons of cobalt lies in a region of more nearly 
e,, character. These estimates are consistent 
with the results of polarized neutron diffraction 
experiments for nickel and cobalt.!»*»* 

The author thanks Dr. E. Callen and Dr. R. A. 
Ferrell for helpful comments on the manuscript. 
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VARIATION OF THE ELASTIC MODULI AT THE SUPERCONDUCTING TRANSITION 


G. A. Alers and D. L. Waldorf 


Scientific Laboratory, Ford Motor Company, Dearborn, Michigan 
(Received May 22, 1961) 


Very recently there has been considerable inter - 
est in the changes in the Debye 6,’ and the zero- 
point energy” associated with the transition from 
the normal to the superconducting state. Of much 
longer standing is the problem of observing the 
discontinuity in the bulk modulus associated with 
the zero-field superconducting transition and re- 
lating it to measurements of the pressure depend- 
ence of the critical field.* A calculation of all 
these quantities can be made if the changes in all 
the elastic moduli of a given superconductor are 
known. This Letter presents measurements of 
the change in all three elastic moduli of the cubic 
superconductors Pb, V, and Nb. Tantalum was 
also studied but will not be discussed because it 
showed impurity-sensitive anomalous behavior 
in both the normal and superconducting states. 

The measurements were made using 10-Mc/sec 
sound pulses in a modified “sing-around” appara- 
tus designed and built by Forgacs* to have a sen- 
sitivity of one part in 10’. Changes in the three 
elastic moduli C; = $(Cy4+Cy9+2Cyqq), C=Cyg, 
and C’= $(Cy, -Cy9) were directly determined. 
The Ta and Nb single crystals were purchased 
from the Linde Company and were nominally 99.9% 
pure. The vanadium sample was prepared in this 
laboratory to a purity of 99.98%. The lead crys- 
tal was grown from high-purity lead to which 0.02% 
Bi was added in order to keep the electron mean 
free path at helium temperatures much less than 
the wavelength of the sound waves used. With 
these purities, no difficulty with attenuation 
changes at the superconducting transition was 
encountered. Corrections for dimensional changes 
at the transition are too small to be considered in 
all cases.° No correction for a deviation from 
adiabatic sound wave propagation conditions within 
the samples was necessary for V and Nb, but for 
Pb a correction of at most two parts per million 


M(n)-M(s) 


near 4°K might be necessary for the longitudinal 
mode and the bulk modulus. 

The actual measurements of the difference in 
the elastic moduli between the normal and super- 
conducting states are shown in Figs. 1, 2, and 3. 
It is to be noted that the moduli in the supercon- 
ducting state are always less than in the normal 
state. The dip in the C modulus difference for 
niobium near 2°K is probably associated with im- 
purities since vacuum outgassing decreases its 
depth. The sag in AC; /C, for Pb around 4°K can 
be eliminated by a reasonable correction for the 
deviation from adiabatic conditions for the sound 
wave. But the drop in AC; /C; below 2°K would 
remain. 

The method of measurement also yielded the 
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FIG. 1. The difference in the elastic moduli of 
lead between the normal and the superconducting 
states as a function of temperature. 
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FIG. 2. The difference in the elastic moduli of 
vanadium between the normal and superconducting 
states as a function of temperature. 


temperature dependence of the elastic moduli in 
both the normal and superconducting states. In 
the normal state the moduli varied according to 
the relationship 


M =M,(1 -a@T), (1) 


as T approached zero. The values of @ are given 
in column 3 of Table I. Since the elastic moduli 
are second derivatives with respect to a strain 
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FIG. 3. The difference in the elastic moduli of 
niobium between the normal and superconducting 
states as a function of temperature. 


of the total energy of the solid, it seems reason- 
able that this T? dependence can be attributed to 
the conduction electrons which have an internal 
energy proportional to T?. 
Temperature-dependent elastic moduli may be 
expected to make the Debye @ temperature-de- 
pendent. Daunt and Olsen? have shown that a tem- 
perature-dependent @ can produce a large contri- 
bution to the specific heat from the zero-point 


Table I. The adiabatic elastic moduli at 4.2°K and the Debye temperature 9) computed from them for the super- 
conductors Pb, V, and Nb. The other columns are quantities derived from the present measurements as is de- 


scribed in the text. 











Modulus a 6 — Py - Hs) 4 KAe H,?/81 4B/B (B/41)(8H,/ep) 
Element [10'! dyne cm~*] [10~' deg K~*) [°K] (°K) [10-* ev/atom] [10-" ev/atom) [ppm] [ppm] 
Pb C= 6.99 34 105 4.6x10~ 0.45 4.7 4.0 3.2 
3.6 
= 1, e 
Cc 94 57.5 4s 
C’ =0.506 69.2 4.7 
Vv C, =22.18 3.3 400 1.5x10~ 15 6.2 <4 2.1 
C = 4.60 10.3 
C’ = 5.65 4.6 
Nb 1, = 22-50 2.7 277 1.6x10™* 15 17 0.7 2.7 
C = 3.09 12.8 
Cc’ = 6.01 2.3 





678 
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lattice vibrations. They have used this contribu- 
tion to explain the anomalously low specific heat 
observed in superconductors.® Unfortunately, 
their argument neglects other energy terms which 
necessarily compensate for the zero-point energy 
variations. If the specific heat is derived from 
the entropy, which is more easily specified than 
the total energy, it is easily shown that a temper- 
ature-dependent 6 simply introduces a correction 
factor of the form (1 -2d1n@/d1nT) to the usual T° 
lattice specific heat. Using the present measure- 
ments to estimate this correction factor, it follows 
that the contribution to the specific heat from a 
temperature-dependent @ is too small and of the 
wrong sign to explain the observed specific heat 
of superconductors.® 

Even though the zero-point energy does not ap- 
preciably affect the specific heat, it is an impor- 
tant part of the total lattice kinetic energy. Ches- 
ter’ has shown that the difference in the lattice 
kinetic energy between the normal and supercon- 
ducting states should be comparable to the total 
free-energy difference between the two states. 

If one assumes that the relative changes in sound 
velocity observed at 10 Mc/sec are the same at 

all frequencies, then the change in zero-point 
energy between the normal and superconducting 
states can be calculated. Column 6 of Table I 
gives this calculated change in zero-point energy, 
while column 7 shows the total free-energy change 
at 0°K as determined from the critical magnetic 
field. The two energy changes are indeed quite 
comparable. Thus, if all the frequencies are 
changed by a fractional amount equal to the changes 
observed at 10 Mc/sec, then the zero-point energy 
contribution to the lattice kinetic energy must not 
be neglected. 

It is interesting to note that the modulus most 
affected by the transition to supérconductivity is 
the smallest shear modulus (C for Nb and V and 
C’ for Pb and Ta). Since the change in modulus 
is related thermodynamically to the change in the 
critical field with stress,® the present results 
would indicate that a shear stress would be most 
effective in changing the critical field. Pippard® 
has shown that the critical field in cubic super- 
conductors must depend upon a higher power of 





the shear stress than the first. Thus measure- 
ments of the variation of the critical field with 
shear stress should not only be large but also 
nonlinear. This effect has been observed recently 
by Seraphim and Marcus.’° 

Column 8 of Table I gives the discontinuous 
change in the bulk modulus observed at the super- 
conducting transition in zero magnetic field. It is 
related to the pressure derivative of the critical 
field by 


4B/B =(B/4 m)[dH _/aPP at T=T (2) 
Column 9 gives the right side of Eq. (2) as com- 
puted from the dH /dP data of several authors.?>" 
The agreement is excellent for the soft supercon- 
ductor Pb but very poor for vanadium and niobium. 
The failure in the case of the hard superconductors 
is probably associated with impurity effects and 
demands further study. 

The authors wish to thank Dr. A. W. Overhauser 
for his encouragement to make the measurements 
and for his important contributions to the interpre- 
tation of the final data. Mr. R. L. Forgacs of this 
laboratory also deserves special recognition for 
his design of the electronic apparatus which has 
far outperformed its original specifications. 
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FIG. 2. Effect of twist on the magnetic transition 
curve of a single-crystal rod of Sn oriented with 
tetrad axis nearly at 45° to the crystal axis. 


hence the transition curve is broadened but its 
mean position is not displaced. This behavior 
is observed in a reversible way, as shown by 
Fig. 2. 

We wish to thank A. Burgess and P. Roland 
for assistance with the measurements, and 
Y. Budo for assistance in the design and con- 
struction of the apparatus. 
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for example reference 9. 

“The magnitude of the second order effect at 200 atm 
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it may be estimated from the measured value of 
(9°T,./ap*)T, =7. 8 X10~"° deg/atm’, which gives 
(0°H/ap*) 7, = 2 x10~* gauss/atm*, If the observed 
small values of AH, are in fact second order effects, 
the first order coefficient 8; would be even smaller. 

Grenier’ s measurements! have been extrapolated to 
T;, and account has been taken of the angle of the 
stress with the tetrad axis [see (6)]. 

®Two perpendicular directions define the shear, the 
stresses being tangential to the planes perpendicular 
to these directions. 
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VALENCE SPIN-ORBIT SPLITTING AND CONDUCTION g TENSOR IN Sit 
L. Liu 


Institute for the Study of Metals and Department of Physics, University of Chicago, Chicago, Illinois 
(Received May 15, 1961) 


Roth et al. have calculated conduction g tensors 
in semiconductors, and have obtained excellent 
agreement with experiment for Ge*»? and InSb.* 
In the absence of crystal wave functions the cal- 
culations for these crystals and Bi* were based 
on a two-level approximation to the effective- 
mass sum rule. Then the g tensor could be es- 
timated from observed spin-orbit splittings and 
effective masses. Roth also analyzed Si in the 
two-band approximation, but her original esti- 

1 § mate of the isotropic g shift was too large by a 

factor of 3. Elliott pointed out® that because of 

| the proximity of the conduction band edge along 
A to the point X, the crystal spin-orbit splitting 
should be small (at X it is zero by symmetry). 
Our calculations show that the valence spin-orbit 
splitting at A is actually too small by'a factor of 
3 to explain the experimental g shift. Yafet re- 
alized® that in Si the term considered by Roth 
might give the wrong sign for the isotropic g 
shift; our calculations have verified his conjec- 
ture. Finally, recent spin resonance experiments 
by Feher and Wilson’ on strained Si samples have 
shown that 6g,, which is negligible in the two- 
band model, is actually greater than og | . There- 
fore the two-band approximation is inadequate 
for Si, and a detailed calculation using crystal 
wave functions is required. 

Several quantum-mechanical caiculations of 
the g factor of electrons in crystals have been 
reported using cellular methods.*»® Here we 
use the orthogonalized plane wave (OPW) method, 
which is well suited to Si. We begin with the 
spin-orbit splitting of the valence band; the spin- 
orbit splitting was first calculated from OPW 
wave functions by Cohen and Falicov.*° 
OPW crystal wave functions are divided into a 

“smooth” plane wave part and a core part: 





= Da" K+K)E+K)") + Doe" lo, ge) 
K t (1) 


where @ labels irreducible representations, 








waves, and be” is determined from the re- 
quirement that aa be orthogonal to the core 
orbital ¢,. In caneatinn (de lHg.9, |¥i"), where 
Hg , is the spin-orbit operator, we find that 96% 
of the matrix element comes from core-core 
terms. The spin-orbit splitting of the p valence 
bands along k =(A00) is therefore conveniently 
expressed in terms of the 2 core orbitals. For 
instance, at k= 0, 


a )= 31a i 


To5, ‘l2p,)1, @) 


2 x 
(2p, an 


where 


A. =(2/V3)a(111)b 


5! 9 ptt) +V2 a(002)b 


9 p 002) tees 
(3) 


The parameters a and b are taken from the band 
calculation of Kleinman and Phillips.“ Sixty-five 
plane waves were used in the expansion (3). 

Consider first the 2p spin-orbit splitting 
3| (2 py |\H,.6.|2p,) |. Tomboulian and Cady” 
completed the identification of the x-ray emis- 
sion lines 2 p** +2s and 2p”* +2s for the second 
row of the periodic table. Their values for the 
2p spin-orbit splitting for Mg and Si are listed 
in the first line of Table I. By invoking Slater’s 
rule that the missing electron gives an extra 
screening charge 0.3e, we can use Tomboulian 
and Cady’s values for the spin-orbit splittings 
to estimate the 2 psplitting in neutral Mg and Si 
(second line of Table I). Calculated values for 
Mg ** and Si are listed in the third line. It will 
be seen that the discrepancies are quite similar 
for Mg and Si. These discrepancies are not un- 
derstood; in the calculations quoted below, the 
corrected experimental value for the 2p spin- 
orbit splitting is used. 

From (2) and (3) we now find Ag 9 (Tg5) = 0.042 
ev, in good agreement with the infrared value 
0.0441 + 0.0004 ev of Zwerdling et al. A plot 
of A. o As) is given in Fig. 1. 

Starting from the effective-mass (k- p) Hamil- 
tonian and treating H; 9, as a perturbation, Roth* 


(k +K)* is a symmetrized combination of plane has derived the following formulas for the g ten- 
sor in Si: 
4 1 by uy x vz vz 
=6g = quae A A 
08 = 8, Rem U EF (4,10, 145 A 1H, !45 445 Ip, 144) 
u,v Op Ov 
8 x iv iv vz vz 4 
+Re—— >> gz (4, l#, baci 1 (a Ip, 14, (A Ip, 1A,); (4) 
u,v Op Ov 











VoLUME 6, NUMBER 12 


PHYSICAL REVIEW LETTERS 


June 15, 1961 





+t 1 uz pz v2 
+ReT] ls p—e— 64, |P 145 (aS Ie IAS (AS 

















dg = og 
4 E 
ye my Op Ov 
4 1 
u,v Op O 
E 
u,v Ou Ov 
03 « 
3 ma 
§ \ 
— .02 \ 7 
coal 
o ‘ 
q ~*~ 
a Oo! 7, — 
k=O k=(1,0,0)2 


FIG. 1. Sketch of the spin-orbit splitting of the 4; 
valence band of silicon. The dashed curve represents 
interpolation from the calculated results. Note that 
although there is a spin-orbit splitting for the valence 
state I5’ at k=0, the lower split level is associated 
with 4,’. So as far as 4s band is concerned, the split- 
ting is equal to zero at the zone center. 
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FIG. 2. Sketch of the energy bands of silicon along 
[100] axis of the Brillouin zone, after Kleinman and 
Phillips. Superscript ¢ is used to denote the 2p core 
states. 
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Figure 2 shows the relevant energy levels. The 
conduction band edge is believed"® to be at A 
=0.85(27a~"). Note the 2p core levels. 

Because of the selection rule (X4”|H, 9 .|X4”) 
=0, the most important matrix elements of Hg 
are (A5” WH 9/45”) ~(X,4" lHg.9.!X4"") #0. From 
Table II we see that terms of this form, in which 
one of the levels is a 2p core state, make the 
largest contribution to 6g. Comparison with 
Feher and Wilson’s values is made in Table III; 
the agreement is excellent. 

According to Roth*® the “one-valley” interband 
term H, is responsible for donor spin-lattice re- 
laxation for H ||(100). Roth’s expression for H, is 


H, =4 Aste, . oF, + “8 ) +cyclic term} , (6) 


where A is given in terms of the deformation 
potential E, 2 by 

E 
45 (q/!P 14,,)¢4,,| ag My) 


A im E,,7E,, 





eM Y 
(4, IP IAs (ae 1H | 1AD,) 
i. z 
+(A, 1H, 5 145 Kas Ip, 145,)}- (7) 


We have calculated all the matrix elements inA 
except that of the deformation potential. From 
Feher and Wilson’s measured value’ of A (0.44 
+ 0.04), we find that (Ag lEyz |4,) =23 ev. For 
comparison, the intraband shear deformation 


Table I. Spin-orbit splitting of 2p core state. 











Mg Si 
Exp. (ev) 0.27 0.72 
Corr. exp. (ev) 0.22 0.60 
Cale. (ev) 0.17 0.49 
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in the two-band approximation. 


Table II. Relative magnitude of contributions tog tensor. The first term in og i is the only one contributing 








Term in Eqs. (4) and (5) 


Relative magnitude with respect to 




















involving: first term in Eq. (4) 
6 4, ipl.” a,” | > Me, a,*) ? (lp, | 4) 1 
(,l9,145”), ed ere "lp ,14,) 0.7 
(a, 12,145)» ala, 1a.) (as 1p,14,) 0.7 
(a, l2,14,°). wo \H, Z14,%), (.*1p,14,) -2.6 
4, I2, 4.) ; a.” | H, ol a,”), (a*l2, 4) -2.6 
By aylel4g, Aslelas, gL, O14) —~ 
(12,1459. (71m, 14,9, (4, "1p,14) ae 
(a, 19,1457), (71H, 214,” "lp ~0.2 
(a, 12,145), lM, 21a), (‘lp 14) “ 
(310,185), (71a 18,9, ,'lp,|4) <a 
Table III. g tensor for Si. 114, 90 (1959). 
M. H. Cohen and E. I. Blount, Phil. Mag. 5, 115 
6g, 6g, (1960). 
R. J. Elliott, Phys. Rev. 96, 280 (1954), and 
— aw Ratan Way. wae taal . 
Exp. -0,0028 -0.0040 'G. Feher and D. K. Wilson, 1960 Prague Conference 








potential matrix element is Eg =(A, IEyy 1A) 
=7 ev.’” Evaluation of this discrepancy requires 
a detailed theory of deformation potentials. 

The author wishes to express his gratitude to 
Professor J. C. Phillips for suggesting the 
problem and for many helpful comments, to 
Professor M. H. Cohen and Dr. L. M. Falicov 
for interesting discussions, and to Dr. L. Klein- 
man for showing him the Si wave functions. 
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SUPERCONDUCTING CHARACTERISTICS OF SUPERIMPOSED METAL FILMS 


Paul H. Smith, Sidney Shapiro, John L. Miles, and James Nicol 
Advanced Research Division, Arthur D. Little, Inc., Cambridge, Massachusetts 
(Received May 31, 1961) 


The electrical properties of superimposed films 
of superconductors and normal metals have been 
examined in a series of experiments employing 
the techniques of resistance measurement, per- 
sistent currents, and electron tunneling. The re- 
sults indicate that the normal metal, e.g., Ag in 
a Ag-Pb pair, does in fact go superconducting 
and that the system behaves as does a supercon- 
ductor with a transition temperature determined 
both by the constituent materials and by their 
relative thicknesses. It is argued that the obser- 
vations are not due to alloying, diffusion, or sim- 
ilar effects. The experiments are compared with 
those of Meissner’ on the contact resistance be- 
tween crossed plated wires of normal and super- 
conducting metals. 

Measurements of resistance, critical temper- 
ature, critical current, and critical field vs 
temperature were made on long thin films pre- 
pared by vacuum deposition on glass substrates 
at pressures lower than 2x107° mm Hg. Typi- 
cally these samples were made by superimposing 
along the length of an Ag film, 1 cm long and 0.3 
mm wide, a longer Pb film, 0.15 mm wide, to the 
ends of which potential and current connections 
were soldered. Several such samples were de- 
posited simultaneously on the same substrate 
with, if desired, different thicknesses of either 
metal. There was a negligible delay and no ex- 
posure to the atmosphere between successive 
depositions. Table I shows the thicknesses of 
two typical sets of films of “open-faced” and 
“closed” sandwich construction, illustrated in 
Fig. 1(b). In general, the longitudinal resistance 
of the samples was 1 or 2 ohms at room temper- 
ature and lower by a factor of 4 at a temperature 
just above T,. 


Table I. Total thickness in angstroms. 











Sample No. Lead Silver 
1A-1 484 +30 110 +30 
1A-2 490 +30 362 +30 
1A-3 518 +30 886 +30 
1A-4 598 +30 1840 +30 
35-1 508 +30 ove 
35-2 417 +30 7100 +30 
35-3 487 +30 7140 +30 





686 





Resistance was measured by passing a low- 
frequency current through the sample and dis- 
playing its voltage-current characteristic on an 
oscilloscope. Critical temperatures were de- 
termined to within a few millidegrees by the 
first appearance of zero slope, corresponding 
to zero resistance, at the origin of the trace. 
Isothermal critical fields were determined sim- 
ilarly by applying an external magnetic field, 
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FIG. 1. (a) Critical current vs temperature for Ag- 
Pb samples of different geometry and silver thickness. 
The repeat data were obtained one week after the initial 
experiment. (b) Schematic cross section of the samples 
showing geometry and relative thickness. The critical 
current for sample 35-1 (pure lead) was 0.35 ampere, 
and was virtually independent of temperature. 
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tangential to the films and normal to the direc- 
tion of current flow. Critical currents were 
obtained by noting the value of the instantaneous 
current at which the characteristic deviated from 
zero slope. Near the critical temperature 
critical currents as small as 25 microamperes 
were measured. 

Figure 1(a) shows the variation with tempera- 
ture of the critical currents of the samples 
described in Fig. 1(b) and Table I. The progres- 
sive depression of the transition temperature 
(I e* 0) with increasing thickness of Ag for vir- 
tually constant thickness of Pb is to be noted; 
in particular the low transition temperature, 
1.87°K, of samples 35-2 and 35-3. A progres- 
sive reduction in the critical currents, ata 
given temperature, also occurs with increasing 
thickness of Ag, but variation in total thickness 
complicates comparison between samples. All 
the curves show a linear relationship between J 
and T, except for the gradual approach to the T 
axis at small values of J,. The latter may be 
indicative of intermediate-state effects. The 
reproducibility of the observations may be judged 
from the two sets of data points shown in Fig. 1. 
One set was obtained after cooling the specimens 
to 80°K within 10 minutes of deposition; the sec- 
ond set, in one case, was obtained one week 
later after storage at room temperature. No 
obvious change occurred in the properties of 
the films. So far, only preliminary data on the 
critical fields have been obtained but, as might 
be expected, the critical fields for the superim- 
posed films are lower than for pure Pb films at 
the same temperature. 

The systematic change of the transition tem- 
perature with increasing Ag on a constant thick- 
ness of Pb, the reproducibility, and the very low 
transition temperatures (e.g., the transition 
temperatures of samples 35-2 and 35-3 were 
approximately one-fourth that of pure Pb) make 
it unlikely that the effects are secondary in na- 
ture and the result of impurities, strains, dif- 
fusion, or alloying. All metals used were 
99.999% pure. The nature of the effects was 
unchanged by the type of construction (“open- 
faced” or “closed” sandwiches) or by the se- 
quence of deposition. The data did not change 
with the passage of time. Finally, the Pb-Ag 
phase diagram is well known? and the supercon- 
ducting properties of the alloy system have been 
studied.* The solid solubility of Ag in Pb is ex- 
tremely small (a maximum of 0.2 at.% at 300°C) 
and that of Pb in Ag is modest (a maximum of 


2.8 at.% at 600°C, dropping to 0.8 at.% at 300°C). 
Further, any solid alloy of Pb and Ag is a mix- 
ture of the two solid solutions. Allen has shown 
that the transition temperature of Pb is substan- 
tially unaltered by the addition of Ag up to a Ag 
concentration of about 80% whereupon the mix- 
ture ceases to be superconducting. His results 
suggest that the superconductivity of the alloy 

is due to the presence of a continuous phase of 
Pb containing some Ag in solution; when this 
phase is not present in sufficient amount to form 
a continuous path through the specimen the alloy 
remains normal. It follows, therefore, that the 
maximum amount of Ag in Pb has a negligible 
effect on the transition temperature of Pb, where- 
as Ag containing the maximum amount of Pb re- 
mains normal. 

The results reported in the preceding para- 
graphs indicate that a pure “normal” metal film 
in immediate contact with a film of a supercon- 
ductor will depress the transition temperature 
of the superconductor, the depression being 
greater for increasing thickness of the normal 
metal. To determine if the “normal” metal ac- 
tually becomes superconducting, two additional 
sets of experiments were performed. 

In the first, Ag films of various thicknesses 
were sandwiched between Pb films in a persis- 
tent current loop, arranged so that the current 
flowed through the Ag. Since Sn was present 
elsewhere in the loop, persistent currents were 
observed only below 3.7°K. In one loop contain- 
ing about 2800 A of Ag between two Pb films, 
each about 2000 A thick, a persistent current 
of 2.9 amperes at 3.64°K was observed. Although 
a change in current of a few percent would be 
easily detected, no change was observed ina 
period of 13 hours. From prior experience in 
the use of the persistent current loop and the 
associated measuring techniques,’ the inductance 
of the loop is known to be less than 107° henry. 
Accordingly, the upper limit to the loop resis- 
tance can be set at 107*® ohm, whereas the com- 
puted resistance of the Ag (using known area, 
thickness, and bulk resistivity) is greater than 
10~? ohm. Thus unquestionably a substantial 
reduction in the resistance of the Ag was ob- 
served. In a further experiment, a Au film was 
sandwiched between a Pb film and a Sn film. In 
this case, a persistent current could not be gen- 
erated above 3.57°K, well below the transition 
temperatures of both Sn and Pb; thereafter the 
critical current increased smoothly to more than 
2 amperes at 3°K. 
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In the final set of experiments, electron tunnel- 
ing’ from Pb and Sn films through a Formvar 
dielectric layer into the Ag side of a Pb-Ag sam- 
ple was observed. These preliminary experi- 
ments have demonstrated the presence of an elec- 
tron energy gap in the Ag, with a low-temperature 
limiting value of 0.00016 ev in a sample consisting 
of 5700 A Ag superimposed on 1500 A Pb and hav- 
ing a transition temperature above 4.2°K. 

Following measurements of the contact resis- 
tance between crossed plated wires, Meissner? 
has reported the depression of the transition tem- 
perature of a superconductor by an adjacent nor- 
mal metal and the possibility of passing a super- 
current through a thin normal metal film between 
bulk superconducting metals, but the conditions 
under which these effects occurred did not pre- 
clude the possibility of oxide barrier layer forma- 
tion or physical penetration of the metal films. 

He therefore found difficulty in drawing clear-cut 
conclusions regarding the relationship between 
the various observations or their fundamental 
importance. On the other hand, the experiments 
reported in this Letter demonstrate that when a 
film of a normal conductor is placed in direct 
contact with a film of a superconductor, the 


combination exhibits the attributes of a super- 
conductor, including the ability to pass a super- 
current through the otherwise “normal” metal, 


the presence of an energy gap throughout the 
sample,® and a transition temperature which 
depends upon the film thicknesses and the mater- 
ials involved. 

This work has been greatly helped by the advice 
of R. S. Davis and A. E. Slade and the assistance 
of R. M. Sawdo to whom grateful acknowledgements 
are made. 





'H. Meissner, Phys. Rev. 117, 672 (1960). This 
paper contains references to earlier work. 

2M. Hansen, Constitution of Binary Alloys, (McGraw- 
Hill Book Company, New York, 1958), 2nd ed., p. 40. 

33, F. Allen, Phil. Mag. 16, 1005 (1933). 

‘arthur D. Little, Inc., Progress Report, October, 
1960 (unpublished), p. 12. 

5J. Nicol, S. Shapiro, and P. H. Smith, Phys. Rev. 
Letters 5, 461 (1960). 

€g. Shapiro and P. F. Strong, Bull. Am. Phys. Soc. 
6, 121 (1961); P. H. Smith, S. Shapiro, and J. L. Miles, 
Bull. Am. Phys. Soc. 6, 121 (1961). 

I, Giaever, Phys. Rev. Letters 5, 147, 464 (1960); 
I. Giaever and K. Megerle, Phys. Rev. 122, 1101 (1961), 

’Tunneling experiments are being continued to obtain 
simultaneous measurements of the energy gap on oppo- 
site sides of a Ag-Pb sample, this information being 
pertinent to the question of a position-dependent energy 
gap, proposed by R. H. Parmenter [Phys. Rev. 118, 
1173 (1960)]. 








VoLuME 6, NUMBER 12 


PHYSICAL REVIEW LETTERS 


June 15, 1961 





SUPERCONDUCTIVITY IN THE NEIGHBORHOOD OF METALLIC CONTACTS 
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In the previous Letter, Smith, Shapiro, Miles, 
and Nicol’ confirm earlier reports by Meissner? 
and others of a change in the superconducting 
properties of thin metallic films in contact with 
thin films of other metals. Parmenter® has con- 
structed a theory of such contacts, but a bounda- 
ry condition he employs is yet to be justified 
from more fundamental considerations. In this 
note we should like to present a simple micro- 
scopic theory of superconductivity in such con- 
tact neighborhoods based on a modification of 
the parameter [N(0)V ] which occurs in the BCS* 
expression for the energy gap: 


-1/N(0)V 


€9 = 2tw), e (1) 


If two metallic samples—one a superconductor, 
the other not—are placed in contact, the proper- 
ties of the entire material change from that of a 
superconductor in one material to that of a nor- 
mal metal in the other. The range of the inter- 
action between electrons that produces the super- 
conducting state—the interaction due to phonon 
exchange and that due to the screened Coulomb 
repulsion—has been estimated to be about 107° 
cm.° This might suggest that at a contact sur- 
face the change from superconducting to normal 
properties would occur in this very short dis- 
tance. However, due to the large coherence 
distance between zero-momentum pairs, the 
superconducting correlation can extend deep 
into a volume where the interaction between the 
electrons is in fact zero. In this respect the 
situation is similar to that of the deuteron whose 
wave function extends large distances beyond 
the range of the nuclear potential. This creates 
the possibility that thin films of differing metals 
deposited on one another profoundly influence 
each other’s superconducting properties. 

To be specific, we consider two thin metallic 
films in contact over the plane x=0. The left- 
hand film (called 1) has a thickness ¢,, while the 
right-hand film (called 2) has a thickness ¢,. In 
this situation the electron-electron interaction 
is a function not only of momenta and the rela- 
tive coordinate y, but also depends upon the ab- 
solute position of the two electrons in the x di- 


rection, x, and x,: 


V(r,k,...)=v if x, andx,< 0 


=0 if x, or x, are larger than zero. 


(2) 


Due to this electron-electron interaction there 
is a nonzero matrix element Vp-p for scattering 
from a two-electron state labelled by k to one 
labelled by k’. This matrix element, summed 
over all k’ and averaged over ek in the interaction 
region, yields [N(0)V],,, in (1), which determines 
the energy gap and the transition temperature. 
In a detailed treatment, of course, the energy 
gap might be a function of direction as well as 
position. 

The essential observation made here is that 
this average will be decreased if the electron 
normalization volume is increased while the 
electron-electron interaction acts over only a 
part of the volume. This should result in a de- 
crease of the transition temperature of a super- 
conductor in contact with a normal metal. At 
the same time under the proper circumstances, 
the same argument implies that a normal film 
in contact with a superconductor can itself be- 
come a superconductor. 

To illustrate, consider the extremely simple 
case of two metals in perfect contact (no oxide 
barrier between them) with the same Fermi 
energy and the same effective mass. In this 
case, if the film to the right is normal usually, 


(WOW hya= 7 VOW], (3) 


where [N(0)V], is the interaction constant for a 
pure specimen of metal 1 while [N(0)V],,, is 

that for films 1 and 2 in contact. Due to the ex- 
ponential dependence of the energy gap on N(0)V, 
under the above conditions even the thinnest films 
of normal material would produce drastic altera- 
tions of the energy gap in a thin superconducting 
film. 

However, matters are not quite this simple. 
The differing Fermi momenta in two metals 
produce refraction and, for some angles of inci- 
dence, total internal reflection. More important, 
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under usual experimental conditions a chemi- 
sorbed oxygen layer is almost certain to form 
between the two metals. This will create a 
potential barrier of the order of several tenths 
of an electron volt for a distance of several 
angstroms. Such a barrier will tend to separate 
the two materials, as will any mechanical bar- 
rier or separation. 

We therefore expect that the factor t,/(t, + ¢,) 
gives the greatest reduction of the effective 
interaction and that the actual reduction factor 
should have the form, crudely, 


t,/(t, + Bt), 0<B<1 (4) 


where 8 is determined by the barrier between 
films, the difference in effective mass and well 
depth—all of the effects which prevent electrons 
from freely moving from one film to the other. 
Preliminary calculations indicate that a reason- 
able value of 8 will at least crudely reproduce 
the data of Smith et al.’ 

The arguments presented above have as a 
necessary converse the implication that the 
contact region of nonsuperconducting materials 
should become superconducting when in contact 
with superconductors. The effective penetration 
of electrons from one region to another is limited 
among other things by the electron mean free 
path; the further superconducting electrons pene- 
trate into the “normal area” the smaller the en- 
ergy gap should be. However, as there should 
be only one transition temperature for an entire 
sample, one might expect in a lead-silver con- 
tact that the energy gap would vary spatially, 
reaching its minimum at the outer silver surface, 


690 


in spite of the fact that the transition temperature 
remains high. 

The ideas discussed here have many experi- 
mental consequences. It would be of great in- 
terest to measure the critical temperature as 
a function of film thickness when specimens 
have been placed on one another in a high vacuum 
to reduce the surface layer. The influence of 
different effective masses and Fermi momenta 
in the two neighboring specimens as well as of 
the purity of the nonsuperconducting film on T, 
is of interest. Also the variation of T,, with 
surface layer would be interesting, especially 
in the light of recent tunneling experiments. 
More detailed theoretical investigations of these 
and related questions are being pursued at pres- 
ent. 

The author wishes to thank Dr. P. H. Smith, 
Dr. S. Shapiro, Dr. J. L. Miles, and Dr. J. Nicol 
for communication of their results prior to pub- 
lication. He also wishes to thank Professor 
J. Dillon for a very helpful conversation. 
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It has generally been assumed that in those 
regions of the nuclear periodic table where there 
is a considerable prolate deformation, a smooth 
systematics of rotational and vibrational charac- 
teristics would be obtained as protons and neu- 
trons are added. Recently considerable evidence 
has been collecting to show that the smooth trends 
in these systematics are not always obtained. It 
is the purpose of this paper to enumerate these 
deviations in the systematics, to classify them 
according to their position in the nuclear period- 
ic table, and finally to give a tentative explanation 
involving structural changes expected on the basis 
of the Nilsson level systematics. 

In the case of octupole vibrations (that is, bands 
with K =0 andJ=1-, 3-, 5-,...) experimental evi- 
dence at first seemed! to indicate a random rela- 
tionship between the energy of the octupole vibra- 
tion and the neutron or proton number. Recent 
experimental evidence supplied by multiple Cou- 
lomb excitation? and exceedingly careful decay 
scheme work® together with the earlier work* 
determines the scheme of octupole vibrational 
energies versus neutron number for the differ- 
ent isotopes shown in Fig. 1(a). 

It is immediately obvious that a sharp dis- 
continuity or singularity in octupole vibrational 
energies occurs at 142 neutrons. Figure 1(b) 
shows a similar discontinuity in the hindrance 
factor observed in the a decay from the ground 
state of even nuclei to octupole vibrations of 
nuclei in the vicinity of 142 neutrons.® A second 
effect involving 8 and y vibrations has recently 
been pointed out by Durham, Rester, and Class.° 
They have demonstrated a discrepancy in the 
nuclear strength parameter, p, between Th?” 
and U*** and a similar discrepancy in the ratio 
of £0 to E2 transition probabilities from £8 vi- 
brational states between these two nuclei. It 
should again be noted that Th?* has 142 neutrons. 

If one considers systematics of the ground- 








FIG. 1. (a) Energy of the octupole vibrations in the 
actinide region as a function of neutron number. (b) 
The hindrance factor in a decay from the ground state 
of even-even nuclei to the octupole vibrational state 
in the daughter nucleus. The neutron number is that 
of the daughter nucleus. 
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state rotational band instead of vibrational states, 
the effects are much smaller. But, as was first 
pointed out by Hollander,’ the energy of the first 
rotational state goes through a minimum at 142 


neutrons in U and at 146 neutrons in Pu. Further- 


more, a kink in the ground-state to ground-state 
alpha-decay systematics to nuclei of 144 neutrons 
has been observed by Chetham-Strode.® 

One is tempted to propose that the discon- 
tinuity in octupole vibrations at 142 neutrons is 
the result of a change from permanent to meta- 
stable pear shape just at this number of neutrons. 
This proposal has been seriously considered, 
but must be rejected for a number of reasons. 
For example, the expected increase in moment 
of inertia, #?/29, of the rotational band super- 
imposed on the octupole vibration as one goes 
to an unstable pear shape does not occur. Per- 
haps the most significant reason is the need for 
a more general structural change which can ex- 
plain the deviations in the systematics not only 
for vibrational bands, but also for the ground 
state and other nuclear parameters. It seems 


possible that such an explanation may result from 


Table I. 
the ground state are listed in the body of the table. 
and Nilsson, # 


Level assignments for odd neutrons with neutron number ? 139. 
For energies and level assignments see, e.g., Mottelson 


the reconstitution of the core of the nucleus in 
the region of neutron number 142. 

In Fig. 2(a) the Nilsson levels indicate that 
the i ,..|606] level, which was emptied in the 
process of achieving a deformation, may again 
be expected to fill at about 142 neutrons. These 
structural changes (both the emptying and the 
refilling of the i ,,. level) are shown schemati- 
cally in Fig. 2(a) as a dotted line. Such a recon- 
stitution of the core might be expected to change 
the deformation of the nucleus but should not 
change its general prolate shape. Consequently, 
for most nuclear properties one would expect a 
small discontinuous effect. However, for nu- 
clear vibrations determined by a collection of 
intrinsic states which are suddenly drastically 
altered, the effect might be large. 

Indeed, an attempt® to calculate the energy of 
the octupole vibration does suggest a “knee” in 
the curve at 142 neutrons. It remains to be seen 
if more sophisticated calculations can reproduce 
the major discontinuity indicated in Fig. 1(a). 

Perhaps most interesting is the fact that such 
a structural change makes two predictions which 


Energies of the levels in kev above 











Neutron 3/2+ 5/2- 5/2+ 7/2- 1/2+ 5/2+ 7/2+ 9/2- 

number (631) (752) (633) (743] (631) (622) (624) [734] Isotope 

139 \° Th?” 

185 0 390 Th”! 

141 (313) 0 (295) ((400)) us 

0 Pu 

143 0 0.08 u% 

0 145 Pu? 

0 (145) ut 

145 (512) 392 0 286 Pu? 

0 Cm! 

(0) y239 

147 0 172 Pu! 

0 Cm*48 

(0) Pu%3 

((631)) 255 0 394 Cm 





®3. R. Mottelson and S. G. Nilsson, Kgl. Danske Videnskab. Selskab, Mat.-fys. Skrifter 1, No. 8 (1959). 
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FIG. 2. (a) Nilsson neu- 
tron levels for the actinide 
region. The dotted line 
represents the emptying of 
the i;[606] level as addi- 
tional neutrons are added 
with the consequent de- 
crease in deformation. 

(b) The measured defor- 
mation as a function of 
neutron number in this 
vicinity. '° 
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are amenable to test. In the first place, the re- 
filling of the i ,,,.[606] level should produce a 
small but definite decrease in the deformation. 
This may have recently been observed by Bell, 
Bjgrnholm, and Severiens’® and is reproduced 
in Fig. 2(b). A more direct test of the structur- 


al change is shown by the dotted line in Fig. 2(a). 


If the refilling of the i ,,. level comes with the 
addition of two neutrons in going from 139 to 141 
neutrons, then, as indicated in Fig. 2(a), one 
expects to see the $+[633] ground state repeated 
both for 139- and 141-neutron nuclei. Table I 
gives the assignments for odd-neutron levels in 
the actinides.’ The only exception to regular 
filling (see arrow in Table I) is when one goes 
from 139 to 141 neutrons. In both cases the 
ground state is the #+[633] state. Table I shows 
that the energy differences in neutron levels for 
various nuclei differ and therefore that no single 
Nilsson level diagram is adequate to explain all 
nuclei in this region. Consequently one might 
expect that this refilling of the 7 ,,,. level in 
another set of nuclei might occur at a slightly 
different neutron number, e.g., in going from 
141 to 143 neutrons. 

There are, then, a large number of structural 
changes which occur as one successively adds 
neutrons and protons to a nucleus. There is the 
sudden drastic change from a spherical to a de- 
formed nucleus. This occurs as one goes from 
88 to 90 neutrons” and from 86 to 88 protons." 
There is the more gradual return from a de- 
formed prolate nuclear shape to the spherical 
shape in the vicinity of the Os and Ir isotopes." 
Finally, it seems probable that even within the 
region of prolate deformation, there are fairly 
drastic structural changes brought about both 
by neutrons and by protons in both the partial 
breakup and reconstitution of the core of the 
nucleus. Whereas the evidence here has been 
presented exclusively for the refilling of neutron 
levels in the actinides, there is some evidence 
for a similar effect in the rare earths. Each of 
these structural changes ought to result in for- 
bidden transitions between nuclei with different 
structures. Experimental observation of such a 
nuclear Franck-Condon principle in the region 





of 142 neutrons would represent additional evi- 
dence for this proposal. 

The suggestion presented here uses the older 
model in which two particles are placed in each 
of the Nilsson levels until all are used. Levels 
are therefore either full or empty. The diffuse- 
ness of the Fermi surface and the consequent 
partial occupation of many levels in the vicinity 
of the Fermi surface, as suggested by the more 
recent model, should smear out to some extent 
the effect of core reconstitution. 
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The interest in the shapes of resonances in the 
cross sections of nuclides for interactions with 
s-wave neutrons has now extended over a period 
of twenty years. In spite of the long and rather 
intense study of this problem in which the shapes 
of resonances attributed to scattering, fission, 
and radiative capture have been investigated ex- 
tensively, the shapes associated with some proc- 
esses have not been observed. One of these is 
the resonance capture which results in a single 
mode of radiative decay of the compound nucleus, 
i.e., decay by emission of the gamma ray to a 
particular energy level of this nucleus. If this 
can be considered to be a single-channel process 
in the compound nucleus model, the theory says 
definitely that the partial radiation widths which 
correspond to each mode of de-excitation should 
belong to a statistical population drawn from a 
X’ distribution with 1 degree of freedom, and 
the wave amplitudes associated with each mode 
should add coherently to give the cross sections 
of the particular process.! In other words, the 
process should be similar to the elastic scatter- 
ing of neutrons in these particulars. Data al- 
ready presented’ indicate that the observed par- 
tial widths are consistent with the x? distribution 
with 1 degree of freedom. 

Here we are concerned with the second point. 
Coherent addition of the wave amplitudes is usu- 
ally made apparent through interference effects. 
However, in contrast to the situation for scat- 
tering resonances, resonances which are caused 
predominantly by capture have shown none of the 
asymmetries associated with interference effects, 
the explanation being that the interference effects 
cancel because of the addition of very many inter- 
ference terms of random signs. It is therefore 
apparent that an interference effect associated 
with radiative capture can be observed only if 
the study selects a single mode of decay, or only 
a very few. 

The aim of the present experiment was to meas- 
ure the relative intensity of individual capture 
gamma rays as a function of neutron energy in 
the vicinity of the Pt! resonances at 11.9 and 
19.6 ev, both of which are known’ to have J =1. 

The measurements were made with the Ar- 
gonne fast chopper® operated at a speed of 


15000 rpm. The sample, a piece of normal 
platinum 6 in. x12 in. x0.080 in., was mounted 
25 m from the chopper at an angle of about 22° 
to the beam. The over-all time-of-flight reso- 
lution of the system was 0.16 psec/m. The cap- 
ture gamma rays were detected by a Nal(T1) 
scintillator, 6 in. long by 8 in. in diameter, 
whose output was fed to an analog-to-digital 
converter. All of the digitalized pulse heights 
and the time differences corresponding to the 
energies of the captured neutrons were recorded 
on magnetic tape in the ANL three-parameter 
analyzer.* In order to be able to correct for the 
effects of possible gamma-ray summing, two 
experiments were performed, one with the crys- 
tal 54 in. from the center of the sample, and 
another with the crystal 10 in. from the sample. 
However, both this series of experiments and 
others® performed for the purpose of determin- 
ing partial radiation widths show that for plati- 
num no such correction is necessary. 

For the two Pt!®° resonances being considered, 
the easiest transition to observe is the one di- 
rectly from the initial state to the ground state. 
The relative intensity for this transition may be 
studied by means of the ratio R,/Ry, where R, 
is the intensity of pulses (in a given time channel) 
for which the height is in the neighborhood of the 
binding energy and R., is the intensity of pulses 
in a broad band at a much lower energy. It is 
assumed that, because of the complexity of the 
capture y-ray cascade, Ry is a measure of the 
rate of neutron capture. The experimental values 
for the ratio R,/R, are shown in Fig. 1. The 
point (indicated by x ) shown for E = 0.0253 ev 
was the result of a separate experiment per- 
formed with the rotor of the chopper stationary. 
Although the quantity of interest in this experi- 
ment is the relative strength of the ground-state 
transition compared with the total capture in 
resonances of Pt!™ with total angular momentum 
J=1, the ratio R,/R,, is a measure of the frac- 


tion of this transition in the total capture in all 
isotopes of platinum. However, over the energy 
region dominated by the 11.9- and 19.6-ev reso- 
nances, the difference between these ratios is 
not significant; the other isotopes make a signifi- 
cant difference only for the experimental point 
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Along with the data shown in Fig. 1 are three 
curves computed on the basis of a simple multi- 
level formula® applicable to the present case in 
which [« D, Curve A shows the shape obtained 
on the assumption that no interference exists. 


FIG. 1. The relative strength 
of the 7.92-Mev ground-state 
transition in Pt'® following neu- 
tron capture in platinum as a 
function of neutron energy. The 
open circles are the result of an 
experiment in which the distance 
from sample to detector was 
103 in., the solid circles from 
one in which this distance was 
53 in. The indicated errors are 
standard statistical errors. 
Curve A is the shape of the 
variation expected if no inter- 
ference occurs between levels 
of the same J. Curve B shows 
the shape expected if the com- 
ponents associated with the res- 
onances at 11.9 and 19.6 ev 
(indicated by arrows) interfere, 
and the contributions from other 
levels are ignored. CurveC is 
the expected shape if all of the 
levels listed in Table I are in- 
cluded and interference is al- 
lowed between ail levels of the 
same Jd. 


Curve B was computed for the case of inter- 
ference between the 11.9- and 19.6-ev reso- 
nances with no others included. Curve C shows 
the shape obtained if all resonances of the same 
J are allowed to interfere. The parameters used 


Table I. The resonance parameters used in computing the curves of Fig. 1. The partial radiation widths Ty: 


although of the order of magnitude of the absolute partial radiation width T 


only with respect to each other; the absolute values are of no importance here. 


is proportional to the wave amplitude. 


0° must be considered meaningful 
The quantity in the last column 











* * 
0 V/2_p--1)1/2 
Ey (ev) J r (ev) .. (ev) To (ev) lol \jE|“*r-] 
-13.6° 0 0.0302" 0.100" 0.0 0.0 
vind a a b 
-13.0 1 0. 0092 0.100 <0. 00001 +3.2 
11.9 i 0.00344 0.1124 0.00505° +43.5 
19.6 a 0, 0016" 0.100% 0. 00095° -15.9 
67.4! ts 0.0044 0.100° 0.0039° -42.5 











*Chosen so that the contribution to the thermal capture cross section would be consistent with the measured 


value for Pt'®, 
bAssumed. 
See reference 2. 


John R. Waters, Phys. Rev. 120, 2090 (1960). 


“See reference 5. 
See reference 7. 
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in the computation are listed in Table I. A few 
remarks about some of these parameters are 
relevant. Since the known positive energy levels 
do not account for the measured thermal capture 
cross section for Pt!™, a bound level needs to be 
included in the computations to obtain results 
consistent with the experimental data. On the 
assumption that [,,=0.100 ev, a set of param- 
eters E, and l,,° was selected for this level, for 
each possible value of J, which provides the cor- 
rect amount of capture at thermal energies. A 
value of the partial radiation width I, and the 
sign of the wave amplitude were selected for 

the set with J =1 that gave a proper fit to the 
data. (An equally good fit can be obtained with 
the opposite sign and correspondingly different 
Lyo:) Of course, [9 = 0 for J=0. No obvious 
choice can be made between these solutions, so 
all sets of parameters have been included in the 
table. There are known to be two resonances in 
Pt!*® near 68 ev.” Studies of the capture gamma 
rays that result from neutron capture in these 
resonances suggest that only the higher energy 
resonance of the pair hasJ=1.° An area analy- 
sis of some unpublished transmission data, under 
the assumption that fr, =0.100, leads to the esti- 
mate I, = 0.136 ev for this resonance. 

On the basis of the excellent agreement be- 
tween the data and curve C and the obvious dis- 
agreement between the data and the curve A, it 
can be concluded that interference does occur 
between radiative transitions associated with 
neutron resonances of the same J. Curve B 
merely strengthens the argument by showing 
that the dominant features of the data (the maxi- 
mum between the resonances and more important- 
ly the minimum above) are accounted for by con- 


sideration of only the well-known nearby reso- 
nances and are not strongly dependent on the 
parameters of the less well-known distant reso- 
nances. 

Since the number v of channels available for 
this process of radiative capture to a particular 
final state has been the subject of much debate 
during the past few years, it is of interest to 
emphasize that the present results are consist- 
ent with a single-channel process. Furthermore, 
the improvement in the quality of the fit brought 
about by the inclusion of the 68-ev resonance 
strengthens the argument that y=1, and certain- 
ly requires that v be a number very near one, 
if not equal to it, since the interference between 
three neighboring levels would be very unlikely 
if the process could proceed through many chan- 
nels. 





*Work performed under the auspices of the U. S. 
Atomic Energy Commission. 
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Phys. Rev. 104, 483 (1956). 

*L. M. Bollinger, R. E. Coté, and T. J. Kennett, 
Phys. Rev. Letters 3, 376 (1959). 

3L. M. Bollinger, R. E. Coté, and G. E. Thomas, 
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Recently a T=1 resonance in the A-7 system, 
called Y,*, has been observed with a mass of 
1385 Mev.*® Two types of resonances have been 
predicted that might relate this observation to 
other elementary-particle interactions: (1) P,, 
resonances in the A-7 and =-7 systems predicted 
by global symmetry,’ corresponding to the (3,) 
resonance of the 7-N system, (2) a spin-3 Y-7 
resonance resulting from a bound state in the 
K-N system.®"° The position and the width of 
the observed Y,* resonance agree with both 
theories, but since the spin and parity have 
not yet been determined, we cannot distinguish 
between the two theoretical interpretations. 

Global symmetry™ (including a phase-space 
factor) predicts a branching ratio R = (Y,**+ 2°+2*)/ 
gy, = A° +n*) oo as ae >t +7°)/(Y,** ~ A°+n*) 
= 4(0.225)~5%. The K-N bound-state model sug- 
gests values of R considerably larger than 5%. 
However, when nonzero effective ranges are 
taken into account,’* R can become quite small, 
especially if the 2-A parity should be odd. 

To investigate these possibilities, we have 
continued our study of K~-p interactions at 1.15 
Bev/c in the Lawrence Radiation Laboratory 15- 
in. hydrogen bubble chamber by studying events in 
which a = is observed. The total cross sections 
for these interactions are shown in Table I; only 
statistical errors are indicated. The separation 
of D*+2* +7° and 5*+2* +27° events was diffi- 
cult because many of the latter events will also 
fit the first hypothesis. The numbers given in 
Table I and in the mass plots below were cor- 
rected to account for this ambiguity. The cor- 
rection factor was estimated by using our 5*+7 
+n*+77 events. 

To obtain an upper limit for the branching 
ratio, R, we combined the 2-27 events into 
different charge states of the 2-7 system. All 
charged = hyperons were observed; however, 
in the £° cases only two-thirds of the events 
were observable because of the neutral decays 
of the A°. Furthermore, we had estimated that 
about one-third of the £°r*n~ events also fitted 
a An*n™~ interpretation and had been included in 
already published data.’ Consequently each D°n*n- 
event was given a weight of 2.25. The resultant 
mass:spectra are shown in Fig. 1. In the cases 


¥$ 
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of ({-7)* and (Z-m)~ there appears to be no excess 
of events in the region of M=1385 Mev. Using 
the number of (A-7*) and (A-7~) events with 1355 
Mev<M,_,< 1415 Mev from reference 1, and as- 
suming that all charged =-n7 systems in the same 
mass regions of Fig. 1 are Y,*, we obtain Rmax 
<8%. This treatment yields an unrealistic upper 
limit, since there is no evidence of any peaking 
above background. The =++2*+27° events pos- 
sibly misidentified as r+ +2* +7° (or vice versa) 
do not fall into the mass band used in this anal- 
ysis. 

We conclude that the £/A branching ratio R for 
the strong decay of the T=1 Y,* is at most a few 
percent and is consistent with zero. This result 
agrees with the value of R obtained by Bastien 
et al.’* As indicated above, this value of R does 
not rule out either the global symmetry or the 
K-N bound-state model of the Y,* resonance. 

Next, we wish to report the results of our study 
of the three reactions: 


(1) 
(2) 
(3) 


K~+p22U*+n +040", 
K~+p2Ll"4+n*4nt40-, 


K~+p22l°+n°+nt+n-. 


Reactions (1) and (2) are readily identified and 
measured, but reaction (3) cannot be identified 
unambiguously. Accordingly, we discuss first 
the results pertaining to reactions (1) and (2). 

In a search for possible 2-7 resonances, we have 


Table I. Cross sections for the =-producing interac- 
tions at 1.15 Bev/c. 





No. of events Cross sections 





Reaction (uncorrected) (mb) 
K~+p>2-+2" 87 1.40 +0.16 
> rt+a 84 1.34 +0.18 
> ot4+n7 +70 57 0.97 +0.16 
+p +n n9 54 0.83 +0. 20 
+ D4 atta 27 0.97 +0. 20 
> ten + n+ 7! 13 0.18 +0. 06 
+> D-+0t4+ 794+ 779 9 0.12 +0. 05 
~~ ptente nt 19 0.19 +0. 06 
+r +n +n tnt 13 0.12 +0. 05 
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FIG. 1. Mass plots of the charged and neutral 2-7 
systems, including curves representing phase-space 


distributions. (a) Mass of (Z-7)~, from the reactions: 


i +p>Dotnm +n* or D7 +7940", (b) Mass of (Z-7)" " 
from the reactions: K~+p-> 2°+n* +77 or Z Tat + am : 
(c) Mass of (E-7)°, from the reactions: K~+p-> 5° +17 
+r or D7 +n' +0", 





plotted in Fig. 2 histograms of the invariant mas- 
ses of the = and each of the three pions in reac- 
tions (1) and (2). Figure 2(b) refers to the © and 
pion of like charge; Fig. 2(a) to the © and each of 
the pions of unlike charge. For this reason twice 
as many events appear in Fig. 2(a) as in Fig. 2(b). 
The plotted curves are mass distributions ex- 
pected on the basis of a uniform phase-space 
population. The histogram of Fig. 2(b) agrees 
with the phase-space curve, but the = and unlike- 
charged pion distribution appears to exhibit an 
anomaly, indicating a concentration of events 
with a 2-7 mass of about 1405 Mev. A more 
detailed investigation of the distribution of our 
events in the four-particle phase space suggests 
that it is quite unlikely that we are dealing with 
a statistical accident. In fact, in practically 
every event one of the neutral 2-27 systems has 
a mass which lies in the resonance region. If 
one interprets the observed distribution as a 
resonance, its peak corresponds to a mass of 
1405 Mev, and its full width at half maximum is 
about 20 Mev after unfolding experimental errors. 
To investigate further the possibility of a 2-7 
resonance, we studied the 39 two-prong events 
associated with a A that did not fit the K-+p+A 
+n++m7 or K~+p+2°+n*+727 interpretations. 
These events could be: 


K~+p+2D°+nt+a-+7°, (3) 
K~+p7A +nt+7+7°, (4) 
K~+p2l°+nten 40° +70°, (5) 
K~+poA +nt+n74+0°+7°, (6) 


Identification is very difficult because only reac- 
tion (4) is sufficiently overconstrained to permit 
a kinematical fit. Furthermore, most of the 
events that are actually examples of reaction (3) 
fit hypothesis (4), but generally with a larger 

x” value. 

Of the 39 events, 16 had xy”? >2 when kinemati- 
cally fitted to the one-constraint hypothesis (4). 
Most of these events are probably due to reaction 
(3), since a priori only 17% of the events due to 
reaction (4) should have x?2>2. Also, only one 
example of the reaction K>+p>A+at+n~ +t +77 
and no examples of K~+ p+2+t+a*+2*+27+7° 
were observed; thus reactions (5) and (6) are 
probably rare. Even though a kinematical fit 
to hypothesis (3) is impossible, one can obtain 
the invariant mass of the D°-7° system from the 
incident K~ momentum and the measured momen- 
ta of the two charged pions. However, since no 
kinematic constraints can be imposed on such 
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events, the experimental errors will, in general, 
be larger than for fitted events and fluctuate more 
widely. Therefore, the data are better repre- 
sented by ideograms. 

Figure 3(a) shows the ideogram of the mass 
distribution of the 16 events with y?>2. The 
three events with M <1320 Mev can be interpreted 
as the tail of the x” distribution of reaction (4); 
the four events with M > 1450 Mev are probably 
due to reactions (5) and (6). The remaining nine 
events fall into a narrow band centered at about 
1386 Mev and are most probably due to reaction 
(3). The plotted curve is the mass distribution 
of £°-2° systems based on phase space and nor- 
malized to nine events. It is worth emphasizing 
that due to our experimental errors of ~20 Mev, 
the actual peak is expected to be somewhat sharp- 
er than shown. 

Figure 3(b) shows the corresponding distribu- 
tion for the events with y?<2. In order to permit 
a direct comparison with Fig. 3(a), again only the 
measured momenta of the charged pions were used 
to obtain the mass ideogram. The measured dis- 
tribution appears to agree with that expected from 
phase space for An°n+n- events. No anomaly at 
M ~ 1390 Mev is observed. Thus there does not 
appear to be any evidence of the T=1 reso- 
nance in the An°n+ a” data. Furthermore, if one 
fits all 39 events to the An°n+2a~ hypothesis and 
then calculates the A-1°, A-a+, and A-1~ masses 
from the fitted values, there is still no evidence 
for the Y - resonance. In particular, the peak of 
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(b) “ 
FIG. 2. Mass plots of the neu- 
= tral and doubly charged 2-7 sys- 
tems. 
| A. 
1600 


( Mev) 


Fig. 3(a) vanishes. Thus we cannot attribute the 
observed peaks in the mass distribution shown in 
Figs. 2(a) or 3(a) to the Y,* resonance, especially 
in view of the low =/A branching ratio discussed 
above. Because of this and the selection criterion 
used in isolating the events of Fig. 3(a), the nine 
events probably represent a £°-7° resonance 
linked by charge independence to the 2*+-1~ and 
~~-1* resonance already discussed. 

It is easy to show that the branching ratio B 
=N50,0/(N>+,- +N>5 -7+) uniquely determines the 
isotopic spin of the resonance. For T=2, 1, or 
0, we have B=2, 0, or 3, respectively. Neglect- 
ing possible backgrounds, and correcting for 
neutral decays and escape of the A hyperons in 
the £°-1° case, we have 8=0.6+0.2. Hence, the 
isotopic spin of the indicated resonance is zero, 
and we will call it Y,*. 

One difficulty of our interpretation of the data 
is the difference in mass of 19 Mev+ 6 Mev be- 
tween the two peaks of Figs. 2(a) and 3(a). How- 
ever, since there are two identical pions in the 
charged = cases and not in the £° cases, it is 
possible that the effect of Bose statistics could 
cause a shift of the peaks. Also, from the fact 
that the charged 2 can resonate with either of 
the two unlike-charged pions, one would expect 
interference effects between the two resonant 
amplitudes. Another possibility is interference 
between the resonance and nonresonant back- 
grounds. Both these interferences might alter 
the observed positions of the peaks. Electromag- 
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FIG. 3. Ideograms of the missing mass for the 39 
events in which a A and two charged pions were observed 
and neutral pions were also produced. (a) Events with 
x?>2 for the Ar°s*n hypothesis. The superimposed 
curve is the phase-space distribution for the £°n"1*n~ 
reaction normalized to nine events. (b) Events with 
x?<2 for the Arn’ x” hypothesis. The superimposed 
curve is the expected Arn*x phase-space distribution. 


netic mass differences in the 2+-7~, Z--n*, and 
D°-7° have a negligible effect on the position of 
the peaks if one assumes transition probabilities 


that are proportional to the decay momentum. 

The Y,* could also be produced in the events 
in which the final state consists of a = and two 
pions. If it is produced, it should appear in the 
(=-7)° mass plot given in Fig. 1(c). No significant 
peak is observed; however, the number of events 
in this region of the mass plot is uncertain be- 
cause of the difficulty of correcting for D*1*1°n° 
production. The absence of the Y,* in this final 
state (Y,* +7) could be easily understood if the 
interaction took place mainly through the T =0 
initial channel. The equality of the cross sections 
for 2-27 reactions producing = of all charges 
(Table I) is consistent with this possibility. A 
stronger test of this possibility is provided by 
the interaction K° +p (a pure T=1 state) currently 
being studied by Martin et al.* 

We believe that our data for = and three pions 
are most naturally interpreted by invoking a T=0 
x=-n resonance. However, both because of the 
small number of events involved and the complex- 
ity of the final state, we cannot regard the evi- 
dence as conclusive. Evidence for a (=* -1*) 
resonance has been obtained by Eisenberg et al., 
who have studied K~-meson interactions in emul- 
sion and find a peaking in the ([-7)° mass spec- 
trum at 1405 Mev."* This peaking could be attri- 
buted toa Y,*. In addition, Schult and Capps 
have recently invoked a T=0 resonance at a mass 
of about 1410 Mev to explain the hyperon branch- 
ing ratio in low-energy K~-d interactions.*® 

Dalitz and Tuan have shown that the (b~) solu- 
tion for the scattering lengths in K-N low-energy 
interactions will result in a 2-n7 resonance in the 
T=0 state.* Recent values for the zero-energy 
K-N scattering lengths obtained by Dalitz’® using 
the data presented by Alvarez at the Kiev Confer- 
ence,” indicate that this resonance will be at 
141543 Mev, with a half-width (['/2) of about 
20 Mev. If this explanation of the T=0 resonance 
is correct, it should have J= 4; the observed T=1 
resonance could be the resonance predicted by 
global symmetry with J=3."* Dalitz has pointed 
out that the values of the (a~) solution given in 
reference 9 are consistent with both a T=1 and 
a T=0 Yn resonance’®; both of these resonances 
should then have J= $. 
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84. Pais has pointed out to us that if the Y,* should 
turn out to be the resonance predicted by global sym- 
metry, the question arises whether the existence of 
Y,)* could have anything to do with the global symmetry 
as well. This is certainly not the case because, if the 
Y* is related to the global symmetry hypothesis, then 
there should be a corresponding T=} 7-N resonance 
with Q~160 Mev. Thus the existence of a Y)* may in- 
dicate that the assumption of global symmetry is wrong. 
However, another possibility is that this symmetry 
could be valid in the P-wave but not the S-wave pion- 
baryon interaction. 

19R, H. Dalitz (private communication). 





SIGMA DECAY MODES OF PION-HYPERON RESONANCES* 


Pierre Bastien, Massimiliano Ferro-Luzzi,f and Arthur H. Rosenfeld 


Lawrence Radiation Laboratory and Department of Physics, University of California, Berkeley, California 
(Received May 12, 1961) 


During a study of Kp interactions in the Law- 
rence Radiation Laboratory 15-inch hydrogen bub- 
ble chamber, we have analyzed a total of 249 three- 
body reactions of the type 


(1) 


K~+p2 lt+n7+n°, 
and 
(2) 


at incident K~ momenta of 760 and 850 Mev/c. Re- 
actions (1) and (2) are interesting in view of the 
existence of the resonant pion-hyperon state Y,*,' 
which is now known to influence strongly reactions 
such as 


K~+p2l~+nt+n°, 


(3) 


It is of interest to compare the £* production via 


K~+p+2A+at+a™. 
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reactions (1) and (2) with the A production via (3) | 
to obtain the £/A branching fraction R for Y,*. | 
Another interesting feature of reactions (1) and 
(2) is related to the fact that the dominant decay 
mode Y,* +A+z is accessible only to the Y,* of 
isotopic spin 1. However, pion-hyperon resonances 
in other isotopic spin states such as Y,* and Y,* 
can decay into +7. Alston et al. have already 
reported evidence for a singlet resonance Y,*, and 
we present more Y,,* data in the second half of 
this Letter. Additional data on Y*° + D++z2* have 
also been reported by Eisenberg et al.° 
On the subject of the branching ratio of Y,*, both 
current theories agree in predicting values of R 
largely undetermined but generally “small.’*** In 
particular, global symmetry favors R of the order 
of a few percent with an upper limit of 25%,* while 
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FIG. 1. Upper histogram: Sum of the mass spectrum 
for the =*z° and =~7° systems based on 49 reactions (1) 
and 76 reactions (2) at Px =850 Mev/c. Lower histo- 
gram: Sum of £*r~ and =~7* masses from the same 
reactions, The dashed curve represents phase space 
normalized to all 125 events. The solid curve allows 
for the subtraction of the shaded areas in the Y* peaks. 
For clarity, only the phase space for reaction (2) is 
shown; the difference from that for reaction (1) is 
negligible. The shaded area for charged =7 ranges 
from 1380 to 1410 Mev, which are the limits seen ex- 
perimentally in the Y,**->A+7* mass distribution’; 
the neutral shading runs from 1380 to 1430, as sug- 
gested by Alston et al.” for the Y9*. 


the KN “bound-state” model yields generally larger 
values of R,* even if no real lower limit exists 
because of the present uncertainty concerning the 
AZ parity and the validity of the zero-effective- 
range theory 50 Mev below the KN threshold. 
Our kinematics program, “Kick,” easily sep- 
arates the three-body reactions (1) and (2) from 
the topologically similar two-body reactions K~ 
+p+2Z*+n*, and from the four-body reactions 
K~+p+Z*+n* +7°+2°. A Dalitz plot of the 125 
fitted three-body events produced at 850 Mev/c 
shows the now familiar highly populated bands 


corresponding to a mass of 1385 Mev for the <7 
system. Figure 1 shows the mass histograms of 
these events. A small sketch at the upper right 
represents the scatter diagram in mass space, 
from which these histograms were projected. The 
shaded area in the conspicuous M* peak is easily 
attributed to the decay mode Y** +D*+7°. (In 
fact, as we shall discuss below, we are forced to 
conclude that it is only a statistical fluctuation.) 
From the same hopeful point of view one can ex- 


plain the broader peak at the right of the M+ his- 
togram as the Y*°+2++n7 band projected side- 
ways (see sketch). The dashed lines in Fig. 1 
represent phase space, normalized to the total 
number of events. They are seen to be a rather 
poor fit to the histograms, thus displaying more 
quantitatively the extra population in the Y* bands. 
In the neutral histogram of Fig. 1 some suggestion 
of a peak appears at 1500 to 1540 Mev, which is 
where r. is expected. However, this is not sup- 
ported either by our data at 760 Mev/c or by the 
data for reactions (1) and (2) of Alston et al. at 
1150 Mev/c,? nor is it borne out by the configura- 
tion of events on any of the Dalitz plots. For the 
rest of this Letter, we shall discuss only the peaks 
near 1385 Mev and assume that they can be ex- 
plained in terms of Y,* and perhaps the singlet 
Y,* which, according to Alston et al., shows up 
in the mass region of 1380 to 1430 Mev, partly 
overlapping Y,*.? Note that the charged peak can 
come from Y,** only, whereas the neutral peak 
is the sum of contributions from Y,*° and Y,*. 
We shall discuss the unique charged peak first and 
determine from it the Y,* branching fraction, R. 
In the 850-Mev/c column of Table I, the number 
of events in the shaded area have been used to 
compute R*, which is defined as 


_n(¥ ** +L *n°) +n(¥*" + 2D -7°) 


Rt = ' 
n(Y**+An*) +n(¥*~ +A) 





Charge symmetry requires 


n(¥ **+ + 5° *) +n(¥ *~ + 2°) 
n(Y**+An*) +n(Y* +An~) ~ 


Values of R° derived at various momenta’? from 
the study of the reaction K~+p+2°+n*+a~ are 
also given in row E of Table I. The 124 events 
at 760 Mev/c have been treated in the same way 
and are displayed in Fig. 2. At this momentum 
neither the Dalitz plot nor the histograms show 
such a convincing peak structure; the solid lines 
in Fig. 2 represent phase space and are seen to 
fit the data adequately. Nevertheless some real 
peaking may be present, particularly because at 


Rt=R°= 
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Table I. Data from K~+p—~ D*+n° +7 
P x(lab) (Mev/c) 760 850 1150" 
A. Total number of events 124 125 111 
B. o(2~n*r) + o(Z*n-7°) (mb) 0.9 +0.2 1.7 +0.3 1.8 +0,2 
y,**-- Dt+r° 
G. Events in shaded area 6+7 20+7 -3+5 
D. R* (%) 2+2 15 +6 -2+3 
E. R%from 2°rtn~) (%) 1+5 -5 +5 unreliable” 
Y*!—> 5 *an* 
F. Events in shaded area 1728 16 +8 9+10 
G. o(K” +p— Y*°+ 7% (mb) 0.12+0.06 0.22+0,11 0.15 +0,17 
H. o(K +p— Y,**+2*) (mb) 2.4 +0,15 1.9 +0,2 3.1 +0,4 








“From Alston etal., reference 2. 


The data of row E are plagued with the difficulty of A vs 2° separation; the fraction of ambiguous events increases 
quickly with beam momentum, and at 1150 Mev/c only an upper limit R°<8% can be stated. 


760 Mev/c the Y* bands in the mass scatter dia- 
gram are less well separated, so that they are 
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FIG, 2. Data from 55 reactions (1) and 69 reactions 
(2) at P K= 760 Mev/c. The same considerations as for 


Fig. 1 apply here. 
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less easily resolved and are also subject to the 
interference effects already well known to occur 
in reaction (3) at the same beam momentum.® 
The five values of the branching fraction R shown 
in rows D and E of Table I can now be combined 
to give an average of (2+2)%. This average must 
be used with some reservation, since the five in- 
put data are not very consistent; their x? value is 
8 (its average expected value is 4), and the prob- 
ability for this fluctuation is only about 10%. How- 
ever, it is consistent with the value =(6+ 2)% ob- 
served from the reaction K~ +d + 2~+7°+p by 
Levine et al.” Therefore, we are inclined to be- 
lieve that there is substantial evidence that R is 
less than a few percent and conclude that the 850- 
Mev/c charged peak is a fluctuation. If R only is 
a few percent, we can show now that very few of 
the events in the neutral peak can be attributed 1 
to Y,*°. The argument is as follows: In the pro- 
duction reaction K~+p+Y*+7, the J=1 channel 
cannot produce any Y,*°, and the J/=0 channel 
gives equal numbers of Y,*°, Y,**, and Y,*~. 
From this it can easily be shown that in order to 
explain the neutral peaks as Y,*°+2*++n*, we 
must assume that the production is dominantly 
in the J=0 channel and R~10%. However, Dalitz 
and Miller have given good arguments in favor of 
the predominance of the /=1 channel,° and further- 
more we have just concluded that R is only a few 
percent. : 
We then interpret the neutral peaks either as a 
statistical fluctuation or as evidence for the sing- 
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FIG. 3. Combined histogram of 305 =*n’° pairs at 
760, 850, and 1150 Mev/c. The neutral pairs from 
Figs. 1 and 2 have been added to those in Fig. 1(c) 
of Alston et al. (1150 Mev/c).” The dashed curve is 
the sum of the dashed phase-space curves of our Figs. 
1 and 2 and Fig. 1(c) of Alston et al.? The solid curve 
is the sum of our solid curves (in which we assume the 
existence of a neutral peak containing 17 events at 760 
Mev/c and 16 at 850 Mev/c) and a corresponding curve 
that assumes 9 events at 1150 Mev/c, as shown in 
Table I. 


let Y,* reported by Alston et al.? The excess 


events in the neutral peaks at each of the three 
momenta are listed in row F of Table I. Thea 
priori probability for this fluctuation is actually 
quite small; however, such probabilities can give 
illusory results, so we prefer to illustrate the 
situation through Fig. 3, which is a combined 


mass histogram of all the currently available 
D*n*n° data. It does seem to constitute mild sup- 
port for Y,*. In this connection two disconcerting 
facts should be noted: (a) The Y,* production 
cross section o(Y,*) (row.G of Table I) is rather 
small in comparison with that of Y,* (row H). 

(b) If Y,* shows up so strongly in the two small- 
statistics, four-body final states reported by 
Alston et al.,” it is disappointing that the larger 
three-body samples do not confirm it more 
strongly. 

This work is part of a study of low-energy Kp 
interactions carried out in collaboration with J. P. 
Berge, O. Dahl, J. Kirz, D. H. Miller, J. J. Mur- 
ray, R. D. Tripp, and M. Watson. We thank L. W. 
Alvarez and many members of his group for their 
help. In particular, we want to acknowledge close 
collaboration with the group who ran the 1150-Mev/ 
c experiment: M. H. Alston, L. W. Alvarez, 

P. Eberhart, M. L. Good, W. Graziano, H. K. 
Ticho, and S. G. Wojcicki. 





*Work done under the auspices of the U. S. Atomic 
Energy Commission. 

tNational Academy of Sciences Fellow. 

‘For a list of references on Y;*, see the review ar- 
ticles by M. H. Alston and M. Fferro-Luzzi, and by 
R.K. Adair, Revs. Modern Phys. (to be published). 

2M. Alston, L. W. Alvarez, P. Eberhard, M. L. 
Good, W. Graziano, H. K. Ticho, andS. G. Wojcicki, 
this issue [Phys. Rev. Letters 6,698 (1961)). 

*Y. Eisenberg, G. Yekutieli, P. Abrahamson, and 
D. Kessler (private communication). 

‘Dp. Amati, B. Vitale, and A. Stanghellini, Phys. Rev. 
Letters 5, 524 (1960). 

5R. H. Dalitz, Phys. Rev. Letters 6, 239 (1961). 

*R. H. Dalitz and D. H. Miller, Phys. Rev. Letters 
6, 562 (1961). 

™R. D. Levine, J. J. Murray, D. H. Miller, and 
J. Schwartz (private communication). 
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INVERSE PHOTOPRODUCTION REACTION 7° +p+y+n IN FLIGHT 






G. Gatti, P. Hillman,* W. C. Middelkoop, T. Yamagata,? and E. Zavattini 
CERN, Geneva, Switzerland 
(Received May 31, 1961) 


The cross section for the photoproduction of 
pions off free neutrons has in the past been ex- 
tracted, somewhat precariously, from analysis 
of photoproduction off deuterium.'»? We have 
now measured the absolute value of the cross 
section for the process 1~+p+y+n, from which 
the free-neutron photoproduction cross section 
is immediately calculable using the principle of 
detailed balance. The experiment was performed 
at 90° c.m. angle and a pion energy of 72 Mev, 
which corresponds to a photon c.m. energy of 182 
Mev, 53 Mev above threshold. The cross sec- 
tion was measured by detecting the neutrons in 
a plastic scintillation counter (10 cm thick and 
17.6 cm diameter) whose efficiency has been 
measured. 

In order to discriminate against the much 
more probable charge-exchange reaction 
1” +p+n°+n, the energy of the neutrons was 
measured by time of flight. The time zero was 
given by the incident pions. In order to reduce 
this and other backgrounds even further, a 
coincidence was also required with the y ray, 
which was detected in a large lead-glass total 
absorption Cerenkov counter. All three pulses 
(incident pion, neutron, and photon) were dis- 
played on a travelling wave oscilloscope and 
photographed. The film traces were read on 
IEP, a track-measuring instrument of CERN, 
after which the CERN Ferranti “Mercury” com- 
puter was used to analyze the results and to 
apply various small corrections. The final time- 
of-flight spectrum is shown in Fig. 1. The three 
peaks originate from (1) the photons from neutral 
pion decay, (2) the neutrons from the inverse 
photoproduction reaction, and (3) the neutrons 
from the charge-exchange reaction. 

The beam intensity was monitored by an argon- 
filled ionization chamber® calibrated at low in- 
tensity against a scintillator telescope. The 
beam composition was measured by a glycerol 
Cerenkov counter and found to be (35.7+ 1.0) % 
pions, (10.7+ 1.0)% muons, and (53.64 1.5)% 
electrons. 

The neutron counter efficiency was measured 
in a separate experiment. The same liquid 
hydrogen target as used in the main experiment 
was placed in a beam of neutrons from 190-Mev 
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proton bombardment of a bakelite target inside 
the cyclotron. The protons knocked out of the 
hydrogen were detected in a telescope whose 
energy threshold and channel width had a spatial 
variation, such that the corresponding scattered 
neutrons nearly reproduced the spatial and 
energy distributions of the neutron incident on 
the neutron counter in the main experiment. 

The neutron counter was placed in this flux 

and the fraction of counts in the proton counter 
which coincided with counts in the neutron 
counter was measured. In this way the loss of 
neutrons which were scattered out in the target 
and target walls in the main experiment was 
automatically taken into account. The neutron 
counter efficiency, at a bias equivalent to about 
5-Mev recoil protons and for an average neutron 











Neutron time of flight in nsec 


FIG. 1. Neutron time-of-flight spectrum. The 
peaks are (1) the photons from decay of neutral pions 
from the charge-exchange reaction, (2) the neutrons 
from the inverse photoproduction reaction with an 
average kinetic energy of 26.9 Mev, and (3) the neu- 
trons from the charge-exchange reaction with an 
average kinetic energy of 13.4 Mev. An average 
time-of-flight path of 1.45 m was used. 
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energy of 26.9 Mev, was found to be (13.6+ 0.7) %. 
The final cross section of 70+7 yb sr is 
based on 505 events (after an empty-target back- 
ground subtraction of 80 events). This 10% total 
error is mainly due to a 6% statistical error, 
a 5% error due to lack of separation between 
the peaks (2) and (3), an estimated 5% error in 
the calibration of the neutron counter efficiency, 
and a 3% error in the measured pion flux. Using 
the principle of detailed balance, a value of 
20.1+2.0 ub sr™ is obtained for a, . This is 
the 90° differential cross section divided by the 
statistical factor W =qq,(1+k/M)~*, where q 
and q, are the pion momentum and total energy, 
k the photon energy (all c.m.), and M the nucleon 
mass (in units u,=h=c=1). Figure 2 shows 
our result together with those for a,* (the proc- 
ess y+p+n*+n), and the values of a,~ obtained 
by Adamovich et al.’ from analysis of the proc- 
ess y+d+1~+p+p. Also shown are the predic- 
tions of dispersion relation theory as given by 
Chew, Goldberger, Low, and Nambu‘ and as 
calculated by Hamilton and Woolcock® with 
parameters f?=0.08 and N‘~) =0. The 17/n* 
ratio at our photon energy becomes 1.34+ 0.15, 
where we used an a,* value of 15.0 10-*° cm? 
sr~! as derived from values of Beneventano et al.° 
and Barbaro et al.” at nearby energies. This 
value of 7~/n* is in good agreement with the 
CGLN dispersion theory. However, it is in- 
teresting to note that the experimental values 
of both a,” and a,* deviate from the theoretical 
curves in the same manner around our photon 
energy. A strong p-wave m-7 resonance, which 
was introduced to explain the nucleon form fac- 
tors®»® and for which an indication was recently 
found from pion-proton inelastic scattering,’° 


is expected to contribute to the photomeson 
production and to modify the CGLN dispersion 
relations. Particularly, it has been shown by 
Ball"! and by Kawaguchi et al.” that the 77~/7* 
ratio at threshold is sensitive to the strength 
of 1-7 interaction, but more experiments on 
a, at lower energies are required to clarify 
this point. 

The guidance and encouragement of Professor 
G. Bernardini throughout the experiment has been 
essential to this work. The warm support of 
Professor P. Preiswerk is acknowledged with 
gratitude. We thank Dr. Goldschmidt-Clermont 
who kindly arranged the use of IEP. The help 
of G. Sicher and B. Smith is gratefully acknowl- 
edged as is the competence of the cyclotron, IEP, 
and hydrogen liquefier staffs. 
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FIG. 2. Values of the 7 and 1* photoproduction 
coefficients a) and a)t, respectively. © - present 
experiment; A.-Adamovich et al., reference 1; e- 
Beneventano et al., reference 6; ()- Barbaro et al., 
reference 7; a-J. G. Rutherglen, J. Walker, D. Mil- 
ler, and G. M. Patterson, Proceedings of the 1960 
Annual International Conference on High-Energy 
Physics at Rochester (Interscience Publishers, New 
York, 1960); O-J. E. Leiss, C. S. Robinson, and 
S. Penner, Phys. Rev. 98, 201 (1955); and as re- 
ported by G. Bernardini, Proceedings of the Ninth 
Annual Conference on High-Energy Nuclear Physics, 
Kiev, 1959 (unpublished); #-G. M. Lewis and R. E. 
Azuma, Proc. Phys. Soc. (London) 73, 873 (1959). 
The theoretical curves in both figures are the dis- 
persion relations of Chew et al., reference 4, as cal- 
culated by Hamilton and Woolcock, reference 5, with 
f?=0.08 and N™=0, 
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COUPLED S- AND P-WAVE SOLUTIONS FOR PION-PION SCATTERING 


B. H. Bransden* and J. W. Moffatt 


CERN, Geneva, Switzerland 
(Received March 27, 1961) 


In a recent paper by Moffat,’ a set of coupled 
equations was derived for the inverse partial- 
wave amplitudes on the basis of the partial-wave 
analysis of Chew and Mandelstam.? We present 
here the results of a numerical iteration of these 
equations when coupled S and P waves are treated, 
and all but two-pion exchange mechanisms are 
neglected. Up to and including P waves there 
occur two S-wave subtraction constants, a 
=A,/(- 4%), one P-wave subtraction constant, 
a,'=A,(- 4), anda constant £,, generated by 
the P-wave threshold behavior. 

By using crossing symmetry, Chew and Man- 
delstam’ have deduced that so long as D and high- 
er waves are small the following five conditions 
hold at v=- 4: 


34,°= 34,7 = -A, (1) 
34,2’ = - $a," =a,", (2) 
a.” - $a,°" = 27a,'+18a,". (3) 


All other derivative conditions which may be 
deduced at the symmetry point depend mainly 
on the omitted higher partial waves. 

The system of equations for the inverse ampli- 
tude was programmed for an iterative solution 
on the CERN mercury computer, supplying ), 
the effective pion-pion coupling constant, as a 
parameter. The constants a,' and é, are de- 
termined by imposing the two crossing-symmetry 
conditions (2). The solutions, as functions of 
the single parameter A, are used to determine 
the quantities entering (3). The extent to which 
the condition (3) is satisfied is a test of the con- 
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sistency of the theory, and, as we shall see, this 
does not impose any limitations on the values of 
» employed. If the solutions are to be consistent 
with the Mandelstam representation, no complex 
poles may occur in Aj;(v). It has been verified 
that no such poles occur for the solutions pre- 
sented in this Letter. 

Seven or eight interation cycles were found to 
be sufficient to determine the inverse amplitudes 
to better than 1%, when starting from given val- 
ues of , anda,’. By varying £, anda,’ in a sys- 
tematic manner until Eqs. (2) were satisfied to 
1-2%, these constants and the calculated phase 
shifts for a given A were found to be of the same 
order of accuracy. In a preliminary calculation® 
the conditions (2) were only satisfied to 15%, 
but this accuracy has been found to be insuffi- 
cient to obtain a quantitatively correct set of 
phase shifts. 

The numerical calculations have shown clearly 
that the range of \ is limited. For each negative 
» (corresponding to attractive S-wave forces) 
with |A| <0.5, a single positive value of ~, can 
be found so that conditions (2) are satisfied. 
These solutions are characterized by a very 
small P-wave phase shift and correspond to those 
computed by Chew, Mandeistam, and Noyes.‘ 
However, these solutions with positive ¢, com- 
pletely fail to satisfy condition (3) and do not, 
therefore, represent cases of physical interest. 
For negative values of ,, other solutions are 
found having a negative value of ¢,. For values 
of |\| 20.45 a bound state in the J=0 state occurs 
(cot6,° becomes negative near zero energy), but 
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FIG. 1. (-5A)[v/(v+ 1)? coté,° as a function of pv, 
for different values of A. Note that coté,’ for \= -0.45 
becomes negative near v=0. 


for smaller values of |\| this disappears and 
the solutions are characterized by large P-wave 
phase shifts. The solutions with negative ¢, are 
all found to satisfy the second derivative condi- 
tion (3), and therefore represent acceptable 
physical solutions. 

Solutions of the coupled equations for positive A 
(repulsive S-wave forces) with A> 0.25 fail to sat- 
isfy conditions (2). For A<0.25 the iterative pro- 
cedure adopted becomes unstable and we have been 
as yet unable to determine whether solutions exist. 
The calculated values of -5a[v/(v + 1)?” cot6,°, 
-2a[v/(v+1)}” coté,?, and [v*/(v + 1)}” coté,* 
are displayed in Figs. 1, 2, and 3 for -0.45<) 
<-0.325, and the corresponding values of a,* 
and é, are shown in Table I. Results for smaller 
values of |A| and full details of the calculations 
will be published elsewhere. It should be noted 
that the J=0 phase shifts are much larger than 
I=2 phase shifts at low energies, and that the 
S-wave scattering lengths become smaller as 
the resonance position moves towards smaller 
energies. (This can be observed for the S-wave 
scattering lengths a), a2, which have been in- 
cluded in Table I.) For |,|< 0.45, cot6,’ passes 
through zero from positive to negative values, 
this giving rise to a P-wave resonant cross 
section. The resonance positions vp and the 

















FIG. 2. (-2A)[v/(v+1)]” cots,” as a function of pv, 
for different values of A. The curve for A= -0.325 is 
close to that for A=-0.35 and is not shown in the inter- 


ests of clarity. 
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FIG, 3. [v3/(v+1)]”* coté,! as a function of v, for 
different values of A. 


709 








VotuME 6, NuMBER 12 


PHYSICAL REVIEW LETTERS 





June 15, 1961 





Table I, Calculated values of a,', €;, resonance parameters vp, I, and S-wave scattering lengths @, a>. 








A a,! Ey VR r Qo Q, 
-0.45 -0. 0529 -21.2 eee oe -68.0 +1.24 
~0.40 -0.0442 -22.5 7.7 6.0 +33.7 1.05 
-0.375 -0.0413 -22.9 4.8(5) 2.5 16.4 0.96 
-0.35 -0. 0397 -23.4 3.8 1. 7(5) 10.8 0.87 
-0.325 -0.0373 -24.1 3.0 1.0 7.5 0.79 





corresponding total widths [ are also included 
in Table I [the total width fr is related to the 
reduced width y by F =(vp*/(vp+1))“y]. 

The P-wave resonance suggested by Frazer 
and Fulco® was at v, =1.5 with a reduced width 
y =0.4, but more recently Bowcock et al.® and 
Frautschi’ have found that on the basis of 7-N 
scattering the resonance position Vp should be 
further out at vp =4.6 with a reduced width y 
=0.2 and y=0.4, respectively. The value v 
= 4.6 is given by our calculated phase shifts for 
|A| =0.37 with an associated reduced width 
y=0.5. The calculated S-wave scattering lengths 
do not agree with the analysis of 3-pion 7 decay,*”® 
but as this analysis is based on the assumption of 
small S-wave scattering lengths (i.e., a,, a,<« 1), 
the calculated values are not necessarily in dis- 
agreement with experiment. 

The authors wish to record their gratitude to 
CERN for its hospitality while this work was in 
progress, and to thank the staff of the CERN 
Mercury computer for their cooperation. 
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ERRATUM 





INELASTIC SCATTERING OF 18.9-Mev NUCLE- 
ONS FROM THE 9.6-Mev STATE OF C™. E. 
Bradford and B. A. Robson [Phys. Rev. Letters 
6, 550 (1961)]. 


The calculation giving Fig. 1 contains an error. 
The results of a corrected calculation for a, 
=0.917, a,=-0.399, a,=0 are shown in the 
following diagram. The conclusions of the Let- 
ter remain unchanged. 
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FIG. 1. Inelastic nucleon scattering from the 9.6- 
Mev state of C’?, The points are the experimental re- 
sults of Peelle for 18.9-Mev protons. The theoretical 
curves are for 18.9-Mev neutrons. The curves for 3~ 
and 1 are normalized at 60° and 50°, respectively. 711 
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ABSTRACTS 


In this section are printed the abstracts of Articles 
that have been forwarded to The American Institute of 
Physics for publication in THE PHYSICAL REVIEW. 
In quoting information obtained from this section be- 
fore the appearance of the corresponding Article, ref- 
erence should be made to “Physical Review (to be 
published)” rather than to this Journal. 








NONLINEAR, TIME-DEPENDENT, PLASMA 
OSCILLATIONS. David Montgomery, University 
of Wisconsin, Madison, Wisconsin (Received 
April 6, 1961). 


The Laplace-transform technique employed by 
Landau to solve the problem of the first-order 
motions in an unbounded, rarified, electron plas- 
ma is modified to solve the problem to arbitrari- 
ly high order. The transforms of the nth-order 
contributions are expressible in terms of convolu- 
tion integrals involving only terms up to order 
n-1. The method is applied to second order for 
the case of square-integrable disturbances. 


ANOMALIES IN IONIZATION COEFFICIENTS | 
AND IN UNIFORM FIELD BREAKDOWN IN AR- 
GON FOR LOW VALUES OF E/p. D. E. Golden 
and L. H. Fisher, Department of Physics, New 
York University, University Heights, New York, 
New York (Received March 21, 1961; revised 
manuscript received May 4, 1961). 

Prebreakdown ionization currents in argon 
have been measured in uniform fields for low 
values of the ratio of field strength to pressure 
E/p [5 to 12 volts (em mm Hg)"']. Currents ob- 
tained with varying electrode separation d at con- 
stant E/p and constant p could not be analyzed to 
yield values of the Townsend coefficients a/p 
andy. Currents obtained with varying p at con- 
stant E/p and constant d could be analyzed to 
yield values of a/p andy, but such currents 
yielded coefficients which depend ond. The de- 
pendence of the values of a/p ond is attributed 
to the production of highly excited atoms by reso- 
nance radiation at some distance from the posi- 
tions where the electrons lose their energy; these 
highly excited atoms then produce molecular ions 
and electrons in collisions with ground-state ar- 
gon atoms. The secondary mechanism and the 
dependence of y on d is associated with resonance 
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radiation. Sparking potential measurements in 
argon made by varying both p and d for values of 
pd corresponding to breakdown for the above range 
of E/p show deviations from Paschen’s law. 
Disregarding the above anomalies, the values 
of a/p are smaller than the earlier measurements 
of Kruithof and Penning by as much as a factor of 
15 at E/p=5 (d=4 cm). At this value of E/p @nd 
of d), the value of exp(ad) at breakdown is only 
1.05, and the value of y is about 20. At larger 
values of E/p, the present values of a/p become 
independent of d and approach theirs. The spark- 
ing potentials obtained are significantly larger 
than those obtained by Kachickas and Fisher. This 
is shown to be due to the condition of the cathode 
surface. 


VELOCITY-DEPENDENT CORRELATIONS IN 
THE STATISTICAL DISTRIBUTION OF THE 
ELECTRIC MICROFIELD IN A PLASMA. 

Amiram Ron, Department of Physics and Faculty 
of Electrical Engineering, Technion-Israel Insti- 
tute of Technology, Haifa, Israel, and G. Kalman, 
Department of Physics, Technion-Israel Institute 
of Technology, Haifa, Israel (Received April 6, 
1961). 


The polarization of a plasma in the neighbor- 
hood of a moving ion depends on the ion velocity. 
This affects the distribution of the stochastic 
field acting upon the ion. Correction to the Holts- 
mark distribution due to the complete test parti- 
cle—field particle correlation including this dy- 
namic effect is calculated up to the order e?. The 
result is: (1) shift toward smaller fields, (2) aniso- 
tropy, and (3) velocity dependence, which is not 
necessarily equal to the zero-velocity effect even 
on the average. 


MAGNETIC ANOMALY IN FeTi0, - aFe,O, 
SYSTEM BY MOSSBAUER EFFECT. S. L. 
Ruby, Radiation and Nucleonics Laboratory, 
Westinghouse Electric Corporation, East Pitts- 
burgh, Pennsylvania, and G. Shirane, Research 
Laboratories, Westinghouse Electric Corpora- 
tion, East Pittsburgh, Pennsylvania (Received 
April 13, 1961). 


The solid solutions (1 -x)FeTiO,-—xFe,O, ex- 
hibit strong ferrimagnetic moments for the com- 
positions x< 0.6, where the Fe and Ti ions are 
ordered in the alternate (111) layers. The anom- 
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aly revealed by the Mossbauer measurements is 
that the ferrimagnetic phase consists of ferri- 
magnetic clusters surrounded by paramagnetic 
media. The size of these clusters decreases 
with increasing temperature or increasing local 
concentration of Ti. The isomer shift of FeTiO, 
is 1.2 mm/sec at room temperature as expected 
for Fe?*, although the quadrupole splitting of 
0.62 mm/sec is smaller than that observed in 
other ferrous environments. 


PRODUCTION OF VACANCIES IN PULSED 
GOLD. James J. Jackson and James S. Koehler, 
University of Illinois, Urbana, Illinois (Received 
April 3, 1961). 


Pulse heating experiments were done on 2-mil 
diameter pure gold wires in an attempt to obtain 
information on the sources for lattice vacancies. 
About a quarter or less of the vacancies are gen- 
erated at the wire surface during the first few 
milliseconds. The fraction is smaller for longer 


times. The majority are produced at dislocations. 


From the high initial rate of production one infers 
that the jog energy in gold is equal to or less than 
1.10 ev. The details of the creation and diffusion 
of vacancies are described. Annealing experi- 
ments on pulsed wires are consistent with the 
idea that dislocations act both as sources and 
sinks for vacancies. 


ELECTRON WAVE FUNCTIONS IN METALLIC 
SODIUM. Joseph Callaway, Department of 
Physics, University of California, Riverside, 
California (Received April 10, 1961). 


Wave functions to order k? are presented for 
electrons in metallic sodium. The calculation 
is an application of the cellular method. The 
empirical potential of Prokofjew was employed. 


ELECTRICAL CONDUCTION IN p-TYPE TITAN- 
IUM SESQUIOXIDE. J. Yahia and H. P. R. Fred- 
erikse, National Bureau of Standards, Washing- 
ton, D. C. (Received March 31, 1961; revised 
manuscript received May 25, 1961). 


The Hall effect, thermoelectric power, and 
electrical conductivity have been measured as 


a function of temperature in crystals of p-type 
titanium sesquioxide. A transition is observed 
at about 450°K. Below this temperature the 
crystals behave like semiconductors, while 
above it the conductivity is apparently metallic. 

The behavior below 450°K is in line with anti- 
ferromagnetic ordering. The effective mass of 
the holes is found to be about 5m, and the aver- 
age Hall mobility 5 cm?/volt-sec at room temper- 
ature. It is claimed that these values, combined 
with the fact that the Hall effect is measurable, 
support the assumption of conduction in a narrow 
3d band. 


HYPERFINE SPECTRUM OF CHROMIUM™*53 IN 
ALO,. R. W. Terhune and J. Lambe, Scientific 
Laboratory, Ford Motor Company, Dearborn, 
Michigan, and C. Kikuchi and J. Baker, Nuclear 
Engineering Department, University of Michigan, 
Ann Arbor, Michigan (Received April 7, 1961). 


Electron nuclear double resonance techniques 
were used to observe the hyperfine spectrum of 
Cr® in Al,O,. Through analysis of the spectrum 
at zero degrees a positive value of 48.5+0.1 
Mc/sec was obtained for the hyperfine coupling 
constant and minus 0.85+ 0.04 Mc/sec for the 
quadrupole coupling constant. From this a value 
of minus 0.03 barn was deduced for the quadrupole 
moment of Cr*S, 


COLOR CENTERS AND RADIATION -INDUCED 
DEFECTS IN Al,O,. Paul W. Levy, Brookhaven 
National Laboratory, Upton, New York (Received 
April 5, 1961). 


The peak energy E, and full width U, in ev, of 
the color centers in Al,O, before irradiation are: 
E,=5.45, U=0.6; E,=4.84, U=0.54; with an in- 
dication of a band at E,>6.2, U>0O.4. Irradiation 
of 3x10* r produces saturation of the gamma-ray- 
induced bands which occur at E,=5.45, U =1.25; 
E,,=3.08, U=1.50, and probably an additional 
band at £,=4.28, U=0.70. Reactor irradiation 
produces bands at E,=6.02, U=0.60; E,=5.35, 
U=0.40; E,=4.85, U=0.54; E,=4.21, U=0.80; 
E,=3.74, U=0.88; E,=2.64, U=0.64; and E, 
=2.00, U=0.44. 

Curves of color center concentration versus 
irradiation time for the reactor -induced bands 
at E,=6.02, E,=5.35, and E,=4.85, can be ac- 
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curately represented by a saturating exponential 
plus a linear increase. This behavior is pre- 
dicted by a simple theory which assumes that 
the color centers are formed by the coloring 

of defects present prior to and increased by 
irradiation and possibly by the coloring of ad- 
ditional centers formed only by radiation damage. 
The measured rate of defect formation could be 
consistent with current radiation damage theo- 
ries; however, several inaccurately known pa- 
rameters preclude a meaningful comparison. 


GYROMAGNETIC RATIO OF PYRRHOTITE. 
G. G. Scott, Research Laboratories, General 
Motors Corporation, Warren, Michigan, and 
André J. P. Meyer, Pierre Weiss Laboratory, 
Institute of Physics, University of Strasbourg, 
Strasbourg, France (Received April 7, 1961). 


The gyromagnetic ratio of pyrrhotite was deter- 
mined by measurements of the Einstein—de Haas 
effect. The value obtained for g’ is 1.9+15%. 


EXTRAORDINARY HALL EFFECT MEASURE - 


MENTS ON Ni, SOME Ni ALLOYS, AND FERRITES. 


Jerome M. Lavine, Research Division, Raytheon 
Company, Waltham, Massachusetts (Received 
April 14, 1961). 


Extraordinary Hall effect measurements over 
a wide range of temperature are reported for 
99.9 Ni, 99.4 Ni, 95 Ni, several 80% Ni alloys, 
Fe,O,, and (NiO), ,,(FeO), .,(Fe,0,). The resis- 
tivity p of all samples over the same tempera- 
ture range is also reported. At high tempera- 
tures, the extraordinary Hall coefficient R, for 
99.9 Ni is proportional to p? as predicted by 
theory. However, the dependence upon p is 
smaller in the 99.9 Ni sample at lower tem- 
peratures and in the samples with higher im- 
purity concentration at all temperatures. The 
80% Ni samples exhibit positive R, while R, in 
the high-concentration Ni samples is negative. 
R, for both ferrites is negative below about 400°C 
and positive above this temperature. 

The theoretical work on metals generally in- 
dicates that R, is proportional to p? and the con- 
stant of proportionality r is weakly temperature - 
dependent and of the order of unity. In the high- 
concentration Ni samples, r is negative, of the 
order of unity, and exhibits a significant tem- 
perature dependence at low temperatures. In 
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the 80% Ni samples, ¢ is positive and of the 
order of 10™*. At room temperature, , is of 

the order of 10~° in Fe,O, and of the order of 
10~° in (NiO), ,,(FeO), ..(Fe,0,). In both ferrites, 
ry is strongly temperature-dependent. The change 
of sign of ry with temperature observed in both 
ferrites has not been previously observed in 
metals or alloys. 


UNPAIRED SPIN DENSITY IN ORDERED Fe, Al. 
S. J. Pickart, U. S. Naval Ordnance Laboratory, 
Silver Spring, Maryland and Brookhaven Nation- 
al Laboratory, Upton, New York, and R. Nathans, 
Brookhaven National Laboratory, Upton, New 
York (Received March 7, 1961). 


A precise determination of the magnetic form 
factor of the ordered alloy Fe,Al has been made 
by diffraction of polarized neutrons from a single 
crystal. Some forty-three reflections in the angu- 
lar range sin@/, <0.9 A“! have been examined 
and show characteristic departures from a smooth 
single-valued function, indicating that the un- 
paired electron density in the unit cell is not 
spherically symmetric. The data are analyzed 
in two ways: by comparison with form factors 
calculated from free-atom Hartree-Fock wave 
functions including crystal-field splitting effects, 
and by two-dimensional Fourier projections of 
the unpaired spin density. The analysis indicates 
that, while the radial spin density of the two types 
of iron atom in the lattice is similar, their orbi- 
tal symmetry is different. The results are dis- 
cussed with reference to various theories of the 
electronic structure in transition metals. 


PRESSURE DEPENDENCE OF THE MICROWAVE 
RESONANCE PROPERTIES OF SOME SPINEL 
AND GARNET FERRITES. I. P. Kaminow* and 
R. V. Jones, Gordon McKay Laboratory, Harvard 
University, Cambridge, Massachusetts (Received 
March 15, 1961). 


The anisotropy field, g value, and linewidth of 
several spinel and garnet ferrites have been meas- 
ured at X band and room temperature as functions 
of hydrostatic pressure to 10* kg/cm*. The crys- 
tals studied include yttrium, ytterbium, and erbi- 
um iron garnet; magnesium ferrite (with differ- 
ent distributions of Mg”* ions on A and B sites); 
and Ni;-,Co, Fe9O4 with z =0, 0.05, and 0.10. 

The pressure dependence of magnetization was 
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measured using magnetostatic mode methods in 
the narrow-linewidth materials, yttrium iron 
garnet, and magnesium ferrite. 

The complexity of the crystal structure and 
magnetic interactions makes any quantitative 
interpretation very difficult. However, the ob- 
servations can be understood qualitatively in 
terms of the volume dependence of the crystal- 
line fields and the exchange interactions. In the 
case of erbium iron garnet, the volume depend- 
ence of the ferric-rare earth exchange constant 
is calculated; and, in the case of nickel cobalt 
ferrite, a simple explanation is offered for the 
observed volume dependence of the Co?* aniso- 
tropy. 

The contribution of thermal lattice vibrations 
to the linewidth in yttrium iron garnet is dis- 
cussed, and the possibility of an anisotropic spin- 
orbit interaction is considered. 


*Now at Bell Telephone Laboratories, Holmdel, New 
Jersey. 


CYCLOTRON RESONANCE AND de HAAS~van 
ALPHEN OSCILLATIONS OF AN INTERACTING 
ELECTRON GAS. Walter Kohn, University of 
California at San Diego, La Jolla, California 
(Received April 5, 1961). 


An electron gas with short-range interactions 
is considered in the presence of a uniform mag- 
netic field. It is shown that (1) the cyclotron 
resonance frequency is independent of the inter- 
action; (2) for a two-dimensional gas, the de 
Haas—van Alphen period is independent of the 
interaction. The low-lying excited states are 
briefly discussed. 


SPECULAR REFLECTION IN THE DIFFRACTION 
OF SLOW ELECTRONS NEAR NORMAL INCI- 
DENCE. J. P. Hobson and I. H. Khan, Radio and 
Electrical Engineering Division, National Re- 
search Council, Ottawa, Canada (Received April 
14, 1961). 


During experiments on the diffraction of slow 
electrons (0-180 ev) at normal incidence, a dif- 
fracted beam was observed corresponding to the 
direction of specular reflection; i.e., straight 
back along the incident direction. This beam had 
properties indicating that it should not be con- 
sidered as a limiting case of beams diffracted 





at other colatitude angles. Maxima in reflected 
intensity of this beam were observed at 4-, 16-, 
and 3l-ev incident energy. The target used was 
a single crystal of tungsten with the (110) face 
exposed. Ultra-high-vacuum techniques were 
employed. 


RESISTIVITY AND THERMOELECTRIC POWER 
OF TRANSITION METALS IN COPPER-GOLD 
ALLOYS. M. D. Blue, Honeywell Research Cen- 
ter, Hopkins, Minnesota (Received April 7, 1961). 


Resistivity and thermoelectric power meas- 
urements on CuAu alloys containing dilute amounts 
of Co, Fe, and Mn are discussed. An analysis of 
these measurements using the thermoelectric pow- 
er formula of Mott leads to values of the scatter- 
ing cross sections of transition elements and 
their derivatives with respect to energy in good 
agreement with the values obtained for these ele- 
ments with Cu or Au as solvent. The results in- 
dicate that contributions to resistance due to dis- 
order scattering and impurity scattering are ad- 
ditive in these alloys. 


SUPERCONDUCTING CRITICAL FIELD OF 
TANTALUM AS A FUNCTION OF TEMPERA- 
TURE AND PRESSURE. C. H. Hinrichs* and 
C. A. Swenson, Institute for Atomic Research 
and Department of Physics, Iowa State Univer- 
sity, Ames, Iowa (Received April 13, 1961). 


The results of precise critical field measure- 
ments on tantalum samples which show “soft” 
superconducting behavior are given along with 
direct measurements of the pressure effect, 

(8H c/®P) 7; as a function of temperature. The 
BCS theory is used as a guide in the extrapola- 
tion of these data to absolute zero from 1.1°K. 
The advantages of using an H? vs T? representa- 
tion for both the critical field and pressure ef- 
fect data are discussed, and it is shown that if 
both sets of data can be represented in terms 

of power series [H? or (8H ,”/aP)7 vs T?] over 

a limited range of temperature, it is then pos- 
sible to write down explicit power-series expres- 
sions for the differences in the thermodynamic 
functions between the normal and superconduct- 
ing states over this same temperature range. 
The electronic contributions to the specific heats 
and the thermal expansions for tantalum are cal- 
culated from the experimental data. The results ~ 
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can be interpreted to show that the thermodynam - 
ic properties of superconducting tantalum cannot 
be expressed in the form of a law of correspond- 
ing states involving the zero-field transition 
temperature. 


*Now at the Linfield Research Institute, McMinnville, 
Oregon. 


SUPERCONDUCTIVITY OF a@ AND 8 MERCURY. 
J. E. Schirber and C. A. Swenson, Institute for 
Atomic Research and Department of Physics, 
Iowa State University, Ames, Iowa (Received 
April 13, 1961). 


Precise critical field measurements and a 
direct measurement of (8H /aP) > as a function 
of temperature were made on physically identical 
samples of a and 8 Hg. The purpose of these 
measurements was to obtain data on the effects 
of crystal structure on the properties of super- 
conductors, and to permit the calculation of 
various thermodynamic quantities difficult to 
obtain in any other way. The critical fields of 
the two phases were found to be identical when 
expressed in terms of the reduced variables 
H/H, and T/T¢. No generalizations of this type 
could be found to explain the pressure effects. 
The advantages of an H? vs T? and 8H?/aP vs 
T? analysis for extrapolation to absolute zero 
are stressed. The critical fields of several 
representative superconductors are compared 
with the critical field predicted by the BCS theory, 
using a plot that emphasizes the detailed shape of 
the curves at low temperatures. This plot also 
can be interpreted in terms of the 6/T. depend- 
ence of the width of the energy gap. The agree- 
ment between calorimetric and critical field 
determinations of the electronic specific heat 
in the normal state is shown to be improved by 
using the H?-T? extrapolation. 


EFFECT OF PERIODIC ADIABATIC TIME 
VARIATIONS ON INTERACTING SYSTEMS. 
H. Suhl, Physics Department, University of 
California, La Jolla, California (Received 
April 3, 1961). 


It is shown that a many-particle system subject 
to periodic adiabatic variation of certain of its 
parameters is to a certain extent equivalent to 
a non-time-varying system with a radically 
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modified interaction between the particles. The 
particular case of an electron gas in a metal is 
discussed in some detail. 


DIFFUSION ALONG SMALL-ANGLE GRAIN 
BOUNDARIES IN SILICON. H. J. Queisser, 

K. Hubner, and W. Shockley, Shockley Transis- 
tor, Unit of Clevite Transistor, Palo Alto, Cal- 
ifornia (Received April 5, 1961). 


Diffusion fronts in samples containing grain 
boundaries are spike-shaped. Velocity of spike 
advance and angle between spike and boundary 
are measured. The “spike velocity method” of 
analysis permits evaluation of two effective 
widths, Wp and Wo» which describe the diffusion 
properties of the boundary. This method has 
been used to analyze data on phosphorus diffusion 
into boron-doped silicon at 1200°C and 1050°C. 

It is concluded that an enhanced diffusion current 
flows along each dislocation of the grain boundary 
over a cross section less than one Burgers-vector 
square. The diffusion current density is about 
300 000 times that of the bulk. This corresponds 
to an energy of 1.5 ev by which grain boundary 
diffusion is favored over the bulk diffusion. This 
enhancement is believed to be caused by enrich- 
ment of phosphorus and also partly by the extra 
concentration of vacancies near the dislocation 
cores. Some possible extensions of the studies 
to include saturation effects at the dislocation 
cores are discussed. 


SOME NEW INTERRELATIONS IN THE PROP- 

ERTIES OF SOLIDS. Johannes N. Plendl, Elec- 
tronics Research Directorate, Air Force Cam- 
bridge Research Laboratories, Office of Aero- 

space Research, Hansom Field, Bedford, Mas- 
sachusetts (Received October 3, 1960; revised 

manuscript received February 13, 1961). 


The lattice vibration spectrum of a solid can 
be characterized by one single frequency, which 
is defined as “the frequency of the center of 
gravity” of this spectrum or simply “center 
frequency.” From its equality with the charac- 
teristic frequency of specific heat we have re- 
cently derived the “center law of the lattice 
vibration spectra.” In assuming now that ina 
lattice, at equilibrium, the hypothetical max- 
imum of vibrational energy (kinetic energy) of 
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an atom (ion) pair equals the total cohesive en- 
ergy (potential energy), and at the same time 
considering the anharmonicity of lattice vibra- 
tions, we derive a basic interrelation between 
center frequency and total cohesive energy. It 
constitutes a substantial extension of the above 
“center law.” Its validity has been illustrated 
for 26 solid compounds of six different lattice 
structures which cover almost the entire range 
of lattice vibration spectra of solids. This inter- 
relation allows a first determination of thus far 
inaccessible data of cohesive energy for solids 
of extremely high sublimation temperatures, 
such as silicon carbide, boron nitride, and the 
two types of diamond. Detailed study of the an- 
harmonicity of lattice vibrations results in addi- 
tional interrelations, such as one between expo- 
nent of repulsion and “related mass” [= reduced 
mass of the vibrating atom (ion) pairs related to 
argon], one between exponent of repulsion and 
“relative compressibility” (change of compres- 
sibility with pressure over compressibility), and 
thus one between “relative compressibility” and 
“related mass.” In combining the two interrela- 
tions of center frequency, namely that with char- 
acteristic temperature and that with cohesive 
energy, we derive an interrelation between char- 
acteristic temperature and total cohesive energy. 
Examination of the relationship between anhar- 
monicity and atomic behavior suggests a classi- 
fication of solids according to their different 
anharmonic force characteristics. The solids 
can have either a soft or linear or hard force 
characteristic, dependent on the configuration 
which they resemble in the periodic chart of the 
atoms. The underlying concept of anharmonicity 
deduced in this paper helps to understand the 
physical properties of solids from an atomistic 
point of view. 


ELASTIC CONSTANTS OF CsBr, CsI, RbBr, 
and RbI. K. Reinitz,* Rensselaer Polytechnic 
Institute, Troy, New York (Received February 
14, 1961; revised manuscript received May 12, 
1961). 


The elastic constants of two body-centered 
halides, CsBr and CsI, were determined as a 
function of temperature from 300°K to 73°K. 
The velocity measurements were obtained 
with an ultrasonic interferometer constructed 
according to the design principles of Williams 
and Lamb. Room-temperature constants of 





RbBr and RbI samples were also measured. 
The values of C,,, C,,, and C,, in units of 10" 
dynes/cm? at 22°C for these salts are given in 
the following table: 





Ci Cie Cu 
CsBr 3.097 0.903 0.7500 
CsI 2.434 0.636 0.6316 
RbBr 3.15 0.493 0.384 
RbI 2.54 0.407 0.276 


The temperature dependence of all the cesium 
salt constants was negative and nearly linear 
over the temperature range investigated. It 
was found that with decreasing temperatures 
C,, increased more rapidly than C,, for the 
two cesium salts examined. 

The elastic constant data of the sodium chlo- 
ride type halides, compiled from the literature, 
are compared with those of the cesium chloride 
type salts. 


*Now at Delco Radio, Kokomo, Indiana. 


INTERPRETATION OF LOW-ENERGY ELECTRON 
DIFFRACTION PATTERNS OF ADSORBED GASES. 
Ernst Bauer, Michelson Laboratory, China Lake, 
California (Received March 27, 1961). 


Recent observations of the adsorption of oxygen 
on nickel are interpreted in terms of double scat- 
tering processes instead of the single scattering 
mechanism which has been used in all previous 
interpretations. This permits a simple explana- 
tion of the observed variation of the intensity of 
the diffraction pattern with electron energy and 
provides a plausible understanding of the early 
stages of adsorption. 


ELECTRONIC STRUCTURE OF THE F CENTER 
IN LiCl. R. F. Wood” and J. Korringa, The Ohio 
State University, Columbus, Ohio (Received April 
14, 1961). 


The electronic structure of the F-center lattice 
defect in LiCl is investigated with calculations 
based on the usual model of the F center proposed 
by de Boer. The ground- and excited-state wave 
functions and energies of the trapped electron 
are determined by two different methods. First, 
the method of linear combination of atomic orbit- 
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als (LCAO) is used. This method is capable of 
yielding good results but the complexity of the 
necessary calculations is great. In an effort to 
avoid this complexity, the method of vacancy- 
centered wave functions is investigated. Very 
simple wave functions are used in this method 
with satisfactory results. The coefficient of the 
hyperfine interaction of the F-center electron 
with the nearest neighbor lithium ion and the 
oscillator strength of the optical transition are 
calculated. The distortions of the lattice in the 
vicinity of the F center are calculated. A very 
small outward movement of the first and second 
nearest neighbors occurs in the ground state. 
The situation in the excited state is complicated 
by the “p-type” symmetry of the F-electron 
wave function. In this case, the two nearest 
neighbor lithium ions located on the symmetry 
axis of the wave function are found to undergo a 
large displacement outward from the vacancy. 
The other four nearest neighbors displace inward 
toward the center of the vacancy. The results 
of the calculations are discussed and detailed 
comparisons with other work of a similar nature 
are given. 


“Present address: Quantum Chemistry Group, 
University of Uppsala, Uppsala, Sweden. 


ANTIFERROMAGNETIC LINEAR CHAIN. L. F. 
Mattheiss,* Solid-State and Molecular Theory 
Group, Massachusetts Institute of Technology, 
Cambridge, Massachusetts (Received February 1, 
1961). 


Many-electron configuration interaction calcu- 
lations have been carried out on a system of six 
hydrogen atoms arranged in a regular hexagonal 
array with a variable lattice spacing. The ap- 
proximate wave functions for this system have 
been expressed as linear combinations of the 
(2x 6)!/(6!)? =924 determinantal functions which 
can be formed from atomic 1s functions. In this 
manner, the effects of ionic configurations con- 
taining as many as three pairs of doubly filled 
orbitals have been introduced into the calcula- 
tions. All three- and four-center integrals have 
been taken into account. The nonorthogonality of 
hydrogenic 1s orbitals localized about different 
atomic sites has been removed by transforming 
to a set of orthonormal Wannier functions. The 
principal result of these calculations is the fact 
that the effects of configuration interaction can 
be represented quite accurately at large inter- 
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nuclear separations in terms of a parameter J’ 
(analogous to a nearest neighbor exchange inte- 
gral) which assumes a negative value in a non- 
ferromagnetic system such as this one. This 
provides a justification for the use of the Heisen- 
berg exchange operator -2J’ (2); 5, *S. 74) to 
describe the magnetic interaction at large sepa- 
rations in this system. In addition, these results 
show that this system of hydrogen atoms is bound 
with respect to six separated atoms, but not with 
respect to three molecules. The ground-state 
wave function is a singlet at ali internuclear 
separations. The general form of the curves 
representing energy as a function of internuclear 
separation shows a striking similarity to those 
obtained for the hydrogen molecule. 


*Now at Cavendish Laboratory, Cambridge Univer- 
sity, Cambridge, England. 


EFFECTIVE EXCHANGE INTEGRAL. L. F. 
Mattheiss,* Solid-State and Molecular Theory 
Group, Massachusetts Institute of Technology, 
Cambridge, Massachusetts (Received February 1, 
1961). 


The magnetic properties of a linear chain of 
monovalent atoms are investigated from the point 
of view of perturbation theory. The many-elec- 
tron wave functions for the system are expanded 
as linear combinations of determinantal functions 
which are eigenfunctions of S* andS,. These 
determinantal functions are constructed from 
orthonormal one-electron orbitals of the Wannier 
type so that the nearest neighbor exchange inte- 
gral is positive definite and approaches zero at 
large lattice spacings. The secular equation is 
set up using the method of the Dirac vector model. 
By means of the Kramers perturbation technique, 
the interaction of ionic states with those arising 
from the ground configuration is represented by 
means of an effective Hamiltonian operator with 
its associated matrix. The results of this treat- 
ment are analogous to those obtained by Paul in 
that an analytic expression is found for an effec- 
tive nearest neighbor exchange integralJ’. This 
quantity is represented as the difference between 
the positive-definite exchange integral and addi- 
tional terms from ionic states. The present treat- 
ment defines in a fairly precise manner the type 
configurations which contribute to this effective 
exchange integral and the limits for which this 
parameterization is valid. The results of this 
analysis are compared with those obtained from 
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recent calculations on a system of six hydrogen 
atoms. 


*Now at Cavendish Laboratory, Cambridge University, 
Cambridge, England. 


STIMULATED OPTICAL EMISSION IN FLUORES- 
CENT SOLIDS. I. THEORETICAL CONSIDERA - 
TIONS. T. H. Maiman,* Hughes Research Labo- 
ratories, a division of Hughes Aircraft Company, 
Malibu, California (Received January 27, 1961; 
revised manuscript received May 17, 1961). 


An analysis of stimulated emission processes 
in fluorescent solids is presented. The kinetic 
equations are discussed and expressions for 
pumping power and effective temperature of the 
exciting source are given in terms of the mate- 
rial parameters. A comparison of excitation 
intensity for three- and four-level systems is 
given. The spectral width of the stimulated 
radiation is discussed with particular attention 
to imperfect crystals. 

*Now at Quantatron, Inc., 2520 Colorado Avenue, 
Santa Monica, California. 


STIMULATED OPTICAL EMISSION IN FLUORES- 
- CENT SOLIDS. II. SPECTROSCOPY AND STIMU- 
LATED EMISSION IN RUBY. T. H. Maiman,* 

R. H. Hoskins,*I. J. D’Haenens, C. K. Asawa, 

and V. Evtuhov, Hughes Research Laboratories, 

a division of Hughes Aircraft Company, Malibu, 
California (Received January 27, 1961). 


Optical absorption cross sections and the fluo- 
rescent quantum efficiency in ruby have been de- 
termined. These data have been used to correlate 
calculations with the analysis of the preceding 
paper. Stimulated emission from ruby under 
pulsed excitation has been studied in some detail; 
the observations are found to depend strongly on 
the perfection of the particular crystal under 
study. A peak power output of approximately 5 
kw, total output energy of near one joule, beam 
collimation of less than 10~? radian, and a spec- 
tral width of individual components in the output 
radiation of about 6x10 A at 6943A have been 
measured, It is suggested that mode instabilities 
due to temperature shifts and a time-varying 
magnetic field are contributing to an oscillatory 
behavior of the output pulse. 


*Now at Quantatron, Inc., 2520 Colorado Avenue, 
Santa Monica, California. 


DIAMAGNETIC SUSCEPTIBILITY OF PYROLYT- 
IC GRAPHITE. D.B. Fischbach, Jet Propulsion 
Laboratory, California Institute of Technology, 
Pasadena, California (Received April 12, 1961). 


The diamagnetic susceptibilities of some pyro- 
lytic graphites deposited at 2100-2300°C have 
been measured at room temperature. As-depos- 
ited sampled had significantly larger suscepti- 
bilities than that of well-graphitized carbons 
or single-crystal graphite. Heat treatment at 
temperatures above 2300°C caused the total 
susceptibility to decrease to a minimum value, 
then rise and level out at a value characteristic 
of graphite, as a function of treatment tempera- 
ture. The relationship of the susceptibility be- 
havior to the structure of the pyrolytic graphite 
is discussed. 


PRESSURE DEPENDENCE OF THE ELASTIC 
SHEAR CONSTANTS OF Li. A. L. Jain, Case 
Institute of Technology, Cleveland, Ohio (Re- 
ceived April 14, 1961). 


The elastic shear constants, C=C,, and C’ 
= $(C,, -C,,), of lithium have been measured asa 
function of pressure at room temperature. The 
measured values of dlnC/d lnr and d1nC’/d Inr 
are -4.1 and -2.8, respectively. This inequality 
of the two pressure variations is contrary to the 
situation found in sodium, where the two coeffi- 
cients were equal. Theoretically, the two coeffi- 
cients are expected to be the same, if only the 
electrostatic interaction between the electrons 
and the ions is responsible for the elastic stiff- 
ness of the metal. The different behavior in 
lithium can be understood in terms of the extra 
contribution to the shear constants arising due 
to the change in the Fermi energy of the elec- 
trons on shearing. 


THERMODYNAMICS OF DIRTY SUPERCONDUC- 
TORS. D. J. Kenworthy* and D. ter Haar, The 
Ciarendon Laboratory, Oxford, Great Britain 
(Received April 7, 1961). 


We use a method due to Tsekhmistrenko to 
eliminate from the Fréhlich Hamiltonian the 
electron-phonon interaction term. We then 
evaluate the thermodynamic properties of a 
superconductor described by this Hamiltonian, 
using a formalism developed by Zubarev and 
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Tserkovnikov which is based on a paper by Bloch 
and De Dominicis. We introduce an extra term 
in the Hamiltonian to take the impurity scattering 
into account and study the effect of this extra 
term on the transition temperature. For the 
product of the mean free path and the relative 
change in the transition temperature we find 
values of 7x10~*, 9x10~, and 8x10~ cm for 
Sn, In, and Al, respectively. 

*Present address: Department of Mathematics, 
McMaster University, Hamilton, Ontario, Canada. 


ENERGY LOSSES IN A MANY-BODY SYSTEM. 
Stanley Engelsberg, Department of Physics, Uni- 
versity of Illinois, Urbana, Illinois (Received 
March 27, 1961). 


The energy loss problem is formulated in such 
a way as to include all losses simultaneously. 
The lifetime and energy losses of a particle in 
a well-defined single-particle state with small 
transition probability are found to be related to 
the self-energy operator. As an illustration of 
the application of the relation obtained, a deri- 
vation of the Bethe sum rule and the Cerenkov 
losses is given for a particle incident on a many- 
body system. 


HARTREE-FOCK ATOMIC WAVE FUNCTIONS 
FROM Cu* TO Kr**. W. W. Piper, General 
Electric Research Laboratory, Schenectady, 
New York (Received April 12, 1961). 


Hartree-Fock atomic wave functions have 
been calculated and tabulated for the ground- 
state configurations of Cu* through Kr**. In- 
terpolation functions for this configuration have 
also been tabulated. 


“REPULSION OF ENERGY LEVELS” IN COM- 
PLEX ATOMIC SPECTRA. R. E. Trees, Nation- 
al Bureau of Standards, Washington, D. C. (Re- 
ceived April 4, 1961). 


Rosenzweig and Porter have shown a “repulsion 
of energy levels” in spacing distributions deter - 
mined from energy levels in complex atomic 
spectra. The present paper extends their work 
by showing that these spacing distributions can 
be determined from calculated positions of the 
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levels in these spectra. Since calculated data 
are available for spectra where the observed 
data are scarce or incomplete, this partially 
overcomes limitations imposed by statistical 
inaccuracy when direct use is made of the ob- 
served data. The equivalence of the two ap- 
proaches is demonstrated by showing that cal- 
culated data for Tall yield the same spacing 
distribution as obtained from observed data for 
Tall and ReI combined. These are complex 
spectra in which a fully developed repulsion 
effect is present. A similar demonstration of 
equivalence is carried out for spectra of RuI 
and Mol, where the repulsion effect is in an 
intermediate state of development. The results 
also indicate that numbers easily evaluated from 
the radial parameters of the theory will indicate 
roughly the degree of repulsion, replacing to 
some extent the need for an explicit calculation 
of the spacing distribution. 


PARAMETERS a@ AND £ IN THE SPECTRA OF 
THE IRON GROUP. Richard E. Trees, National 
Bureau of Standards, Washington, D. C., and Chr. 
Klixbiill Jgrgensen, Cyanamid European Research 
Institute, Cologny, Geneva, Switzerland (Received 
April 4, 1961). 


The parameters a and £8 in the low even config- 
urations of iron-group spectra are qualitatively 
explained as effects of the interaction with con- 
figurations having two 3 p electrons excited from 
the argon-type core to 3d states. Quantitatively, 
the parameters are too large by a factor of two 
when the interaction integrals are taken equal to 
the exchange integrals of Watson’s self-consistent 
field calculations. The parameters are too small 
when the exchange integrals are evaluated from 
observed data in the 3p°3d configuration of Ca Im. 


NEW HAFNIUM ISOTOPE, Hf’. J. Wing, B. A. 
Swartz, and J. R. Huizenga, Argonne National 
Laboratory, Argonne, Illinois (Received April 3, 
1961). 


A new isotope of hafnium, Hf’, has been pro- 
duced by double neutron capture in Hf'® in the 
intense neutron flux of the Materials Testing Re- 
actor (MTR). Mass spectrometric analysis of 
the irradiated hafnium gave a Hf'*/Hf'®° atom 
ratio of 0.00147+ 0.00001. The new isotope de- 
cays with a half-life of (9+ 2)x 10° years by B~ 
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emission predominantly to a 271-kev level in 
Ta!®. The number of 271-kev gamma rays per 
g7 disintegration is 0.84+0.10. The log(/t) for 
the beta transition to the Ta’™ ground state is 
>15, indicating that this transition is at least 
third forbidden. The neutron capture cross sec- 
tion of Hf!** is 40439 barns. 


PARITY CONSERVATION IN NUCLEAR REAC- 
TIONS: SEARCH FOR a DECAY OF THE 8.88- 
Mev STATE IN O**. R. E. Segel, J. W. Olness, 
and E. L. Sprenkel, Aeronautical Research Lab- 
oratory, Dayton, Ohio (Received April 13, 1961). 


A search has been carried out for the parity- 
nonconserving a decay of the 8.88-Mev (2°) state 
in O° by examining the alpha-particle spectrum 
following N*° B decay. An upper limit of [,/T, 
<2x10~° was determined which is shown to lead 
to the estimate that ¢?<2x107!*. The alpha-par- 
ticle group corresponding to disintegration of the 
broad 9.58-Mev (1°) state was observed and the 
log ft for the 8 decay to this state found to be 6.8 
+0.1, the slow transition rate being in accord 
with a shell-model prediction that the 9.58-Mev 
state is due to a three-nucleon excitation. The 
shape of the alpha spectrum was fitted with a 
Breit-Wigner analysis. 


NEW ISOTOPE, Al®°. E. L. Robinson and O. E. 
Johnson, Physics Department, Purdue Universi- 
ty, Lafayette, Indiana (Received April 6, 1961). 


Scintillation measurements were made of the 
beta and gamma radiation from high-purity natu- 
ral silicon targets after bombardment with fast 
neutrons produced by the Li’(d,n)Be® reaction 
(E, < 24 Mev). In addition to well-known radia- 
tions, a beta spectrum with an end point of 
5.05+ 0.25 Mev and two gamma rays with ener- 
gies of 2.26+ 0.03 Mev and 3.52+ 0.03 Mev were 
observed. These gamma rays and the beta group 
decayed, within experimental error, with the 
same half-life, 3.27+ 0.20 sec. The assignment 
of this activity to Al®° and the proposed decay 
scheme are supported by considerations in- 
volving the decay schemes of the well-known 
isotopes produced, half-life studies using por- 
tions of both the beta and gamma spectra, the 
features of experimental beta and gamma spec- 
tra, and nuclear systematics. Strong beta transi- 
tions to the first and second excited states of Si*° 


are inferred by the experimental gamma spec- 
trum and nuclear systematics. A weak beta 
transition (< 2%) to the ground state can not be 
excluded by this investigation. Possible spin 
and parity assignments for the ground state of 
Al®° are 1+, 2+, and 3+. A weak argument is 
made against a spin 1 assignment. The results 
of this investigation can not be used to reduce 
the ambiguity of the spin assignment further. 
The resulting Al°°-Si®° mass difference is 
7.29+0.25 Mev. 


SPECTROSCOPY OF GAMMA RADIATION FROM 
Nd'**, Sr®*, and Pb?”. J. E. Monahan, S. Raboy, 
and C. C. Trail, Argonne National Laboratory, 
Argonne, Illinois (Received March 27, 1961). 


The energies of the cascade gamma rays in 
Nd’ are found to be 1487.04 1.1 kev and 696.7 
+ 0.6 kev and the measured energy of the cross- 
over transition is 2186.0+2.2 kev. The agree- 
ment of these results is used to justify the claim 
of 0.1% accuracy for the scintillation spectrom- 
eter with anticoincidence annulus for the meas- 
urement of gamma-ray energies in the interval 
from 0.5 Mev to roughly 3.0 Mev. An energy of 
570.8+ 0.5 kev is obtained for the low-energy 
radiation from Pb*”, and energies of 1836.2+ 1.7 
kev and 898.74 0.8 kev are reported for two Sr® 
gamma rays. Also measurements are given for 
the relative intensities of the 2.18-Mev, 1.48-Mev, 
and 0.696-Mev gamma rays of Nd’, the relative 
intensities of the 1.8-Mev and 0.898-Mev transi- 
tions in Sr®, and the relative intensities of the 
1.06-Mev and 0.57-Mev transitions in Pb?”. 


MEAN LIFE OF THE 1.61-Mev LEVEL OF Mg”. 
V. K. Rasmussen, F. R. Metzger, and C. P. 
Swann, Bartol Research Foundation of the Frank- 
lin Institute, Swarthmore, Pennsylvania (Re- 
ceived April 10, 1961). 


Nuclear resonance fluorescence techniques 
have been used to measure the mean life of the 
1.61-Mev level of Mg”® and the 1.83-Mev level 
of Mg?®, The exciting y radiation was obtained 
by bombarding metallic Mg?® and Mg”* targets 
with 4,0- and 4.4-Mev protons. For the Mg”® 
level, assumed to be 7/2*, the self-absorption 
of the resonance radiation gives 7 = (2.522:§ 
x107* sec. The angular distribution for the reso- 
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nance scattering was found to be 1+(0.42 + 0.03) 

x P,(cosé@) + (0.03 + 0.003)P,(cosé), where the er- 
rors given are statistical only. For other rea- 
sons it is believed that the correct coefficient of 
the P, term is approximately zero. For the Mg?® 
level, the apparent resonance scattering cross 
section combined with some previous estimates 
of slowing-down times for the excited nuclei 
gives r=(7+3)x107" sec. Further evidence as 
to the collective nature of these nuclei and of 
Al*”’ is discussed. Support is given to the sug- 
gestion of the Chalk River group that the 1.61- 
Mev Mg” and the 2.21-Mev Al”’ levels are the 
7/2* second members of K = 5/2* rotational bands 
based on the ground states. For the Mg” level, 
a spin and parity assignment of 7/2* is required 
to obtain agreement between the quadrupole tran- 
sition probability from these measurements and 
that found by Coulomb excitation. 


TEST OF THE APPROXIMATION IN SOME 
FIRST-FORBIDDEN 2~— 2+ 8 TRANSITIONS. 
R. M. Steffen, Department of Physics, Purdue 
University, Lafayette, Indiana (Received April 
11, 1961). 


The £-y directional correlations of the first- 
forbidden nonunique 2~ ~ 2* g transitions of K*®, 
Sb’”*, and Au’® were investigated and compared 
with the predictions of the approximation, whose 
range of applicability is discussed. Upper limits 
for the contribution of the tensor-type matrix ele- 
ment {B;; to the 8 transitions were estimated on 
the basis of the modified B;; approximation. 

The anisotropy coefficient A,(W) in the B-y 
directional correlation involving the 1.98-Mev 
8 transition of K* varies from A,(1.66) = -0.009 
+ 0.002 to A, (4.60) = -0.049+ 0.002. The energy 
dependence of A,(W) deviates from the predic- 
tions of the — approximation by about 40% over 
the measured energy range. A rough estimate of 
the upper limit for the contribution from the 
{Bjj component is IC4 IB | <0.3(IVo! +1¥,1) 

(in the notation of Kotani). 

The anisotropy factor A,(W) of the B-y direc- 
tional correlation involving the 1.40-Mev 8£ tran- 
sition of Sb’ varies from A,(1.96) = +0.035+ 0.003 
to A,(3.5) = -0.081+ 0.004. The energy dependence 
of A,(W) is well represented by the factor ,(Z, W) 
x (W? -1)/W as predicted by the ~ approximation. 
The upper limit of the / B;; contribution to this 
6 transition is estimated as ICasB, | <0.15/Y,] 
or IC4JBy|<0.21Vol. 
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The anisotropy factor A,(W) of the Au’ g-y 
directional correlation involving the 0.96-Mev 
8 transition varies between A,(1.39) =+ 0.0076 
+ 0.0010 and A,(2.78) =+ 0.02864 0.0010, and its 
energy dependence agrees very well with the 
predictions of the approximation. The upper 
limit for the [Bj; matrix element is estimated 
as IC 4 [By | <0.1/¥,1- 


HIGH-RESOLUTION MEASUREMENTS OF THE 
O**(p, a)N** EXCITATION FUNCTION. H. A. Hill, 
E. L. Haase, and D. B. Knudsen, Palmer Physi- 
cal Laboratory, Princeton University, Princeton, 
New Jersey (Received April 10, 1961). 


The O'*(p, a)N’® activation cross section was 
measured from 12 to 18 Mev with an energy reso- 
lution of 30 kev. The results are essentially the 
same as obtained by Whitehead and Foster and 
Rouse using poorer resolution except for one 
very narrow resonance. As the proton energy 
increases through this resonance, the cross sec- 
tion first rises 7% to a maximum, then drops 
33% to a minimum, and finally rises 6%. The 
peak and valley have a width at half maximum of 
30 kev and the peak and valley are separated by 
60 kev. The proton energy at-halfway down the 
drop was determined using the limp-wire tech- 
nique to be 14.600+ 0.020 Mev. This resonance is 
well suited for calibration purposes. A qualitative 
interpretation of the results is made using the 
“cluster model.” 


O**( p, a)N?* ANGULAR DISTRIBUTIONS AT 13.5- 
18.1 Mev. Donald R. Maxson,* Palmer Physical 
Laboratory, Princeton University, Princeton, 
New Jersey (Received April 10, 1961). 


Angular distributions of alpha particles from 
the O'*(p, a)N'* ground-state reaction were meas- 
ured with an ionization chamber at 10 bombarding 
energies from 13.5 to 18.1 Mev. The angular dis- 
tributions are oscillatory but not of the form pre- 
dicted by the plane wave pickup or knock-on theo- 
ries, and the variation with energy is more pro- 
nounced than would be expected for a simple di- 
rect reaction. The excitation curve has a mini- 
mum at Ep~ 16.5 Mev, and the angular distri- 
butions are markedly different above and below 
that energy. The O'*(p, a)N'**(2.4-Mev) reaction 
is also strongly energy dependent, and the 
O'*(p, p)O** elastic scattering cross section is 
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quite energy sensitive at large angles. The 
energy dependence of the scattering cross 
section at 125 degrees appears to be correlated 
with the O'°(p, a)N* excitation function. 


*Present address: Department of Physics, Brown 
University, Providence 12, Rhode Island. 


EXCITATION FUNCTIONS FOR LITHIUM-6 IN- 
DUCED REACTIONS ON ALUMINUM-27. Inge- 
Maria Ladenbauer,* Ivor L. Preiss, and Carl E. 
Anderson,f Yale University, New Haven, Con- 
necticut (Received April 5, 1961). 


Excitation functions for a number of Li®-in- 
duced reactions on Al®’ have been studied using 
stacked-foil techniques and the Li® ion beam 
from the Yale heavy-ion accelerator. Excitation 
functions corresponding to radioactive residual 
nuclei P*?, P%°, si®*, Al?®, Al®, Mg’’, Na”, and 
Na”? have been measured in the Li® range from 
1 to 63.3 Mev. The data strongly suggest that 
the P**, P®*°, and Si* result from compound- 
system processes and the Na” and Na™ from a 
predominant direct-knockout process. In the 
cases of Al” and Al”® both compound-system and 
direct-pickup reaction amplitudes contribute to 
the reaction yield. 


*On leave from the University of Vienna, Analytical 
Chemical Institute, Vienna, Austria. 

t Present address: Applied Physics Laboratory, The 
Johns Hopkins University, Baltimore, Maryland. 


NEUTRONS AND GAMMA RAYS FROM THE 
BOMBARDMENT OF OO" by He’. K. L. Dunning 
and J. W. Butler, Nucleonics Division, U. S. 
Naval Research Laboratory, Washington, D. C. 
(Received April 11, 1961). 


The threshold energy for the O'*(He*, n)Ne’® 
reaction has been measured. The value obtained, 
3.811+ 0.015 Mev, determines the mass of Ne™® 
to be 18.011446+ 0.000014 amu (O”* standard, 
1960 mass tables). The slow-fast ratio method 
for the observation of neutron thresholds was 
employed at bombarding energies from the 
ground-state threshold to 5.6 Mev, corresponding 
to a region of excitation in the residual nucleus 
from zero to 1.5 Mev. No excited states in Ne’® 
were identified. The bombardment of O'* by He® 
also produced the reactions O'*(He°, p)F"® and 
O'*(He®, @)O®. Energy spectra were obtained by 


means of a scintillation spectrometer for gamma 
rays resulting from certain transitions in F** and 
o'’. For 4.5-Mev He® particles impinging on a 
1-Mev thick target of TiO,, gamma rays of the 
following energies were observed and attributed 
to F**; 0.652+0.007, 0.939+ 0.005, 1.0414 0.005, 
1.17+ 0.01, 1.61+0.02, 1.68+0.02, 2.09+0.01, 
2.51+0.01, 2.65+0.05, 3.06+0.05, 3.35+0.10?, 
and 3.84+0.10 Mev. The following gamma rays 
were also observed and attributed to O°: 5.25 
+0.05, 6.22+0.10, and 6.87+0.10 Mev. 


C4(d,n)N*® REACTION. Ren Chiba,* University 
of Wisconsin, Madison, Wisconsin (Received 
January 25, 1961). 


Differential cross sections for ground-state 
neutrons from the C'*(d,n)N*® reaction have been 
measured at Eg=3.53 Mev and 2.786 Mev by a 
neutron spectrometer. A stripping peak implying 
ly =1 is observed. The excitation function of the 
ground-state neutrons has been measured from 
Eq=1.2 Mev to Eg=3.53 Mev. A number of reso- 
nances were found corresponding to virtual states 
of N**, The angular distributions of neutrons as- 
sociated with the 5.28-, 5.31-Mev doublet (unre- 
solved) and 6.33-Mev levels of N'® were meas- 
ured by photographic emulsion technique. A 
stripping peak characteristic of 1, =0 corre- 
sponds to the unresolved doublet, and lp =1 to the 
6.33-Mev state. The excitation function for all 
neutrons from C**(d,n)N*® was measured both by 
“slow” and by “fast” neutron counters. Several 
possible slow-neutron thresholds were found cor- 
responding to excited states of N’®. The sensi- 
tivity of the slow-neutron threshold technique was 
checked by the O**(d,n)F"”. Accurate values of 
these thresholds are reported. 


*Present address: Nippon Atomic Industry Group 
Company, Ltd., Tokyo, Japan. 


EVIDENCE FOR AN ISOMERIC STATE OF Y®. 
W. L. Alford, D. R. Koehler, and C. E. Mande- 
ville,* Army Rocket and Guided Missile Agency, 
Redstone Arsenal, Alabama (Received April 6, 
1961). 


The recently reported activity induced by neu- 
tron bombardment of niobium has been produced 
by 14-Mev neutrons and by neutrons of energy 
less than 6 Mev on both niobium and zirconium. 
In each case, chemical separation showed the 
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activity to be due to an isotope of yttrium. Two 
coincident gamma rays having energies of 0.200 
and 0.485 Mev and a half-life of 3.14 0.1 hours 
were observed; these observations were in agree- 
ment with earlier results. The activity appears 
very similar to that which has been previously 
attributed to the decay of Y*. However, thresh- 
old considerations and the failure to observe by 
means of a thin-window Geiger counter any beta 
emission associated with this gamma activity 
point to an isomeric state of Y®. Experiments 
with separated isotopes of Zr® and Zr™ support 
this assignment. 


* Present address: Physics Department, University 
of Alabama, University, Alabama. 


FISSION OF GOLD WITH 112-Mev C” IONS: 

A YIELD-MASS AND CHARGE -DISTRIBUTION 
STUDY. H. Marshall Blann,* Lawrence Radia- 
tion Laboratory, University of California, 
Berkeley, California (Received October 4, 1960; 
revised manuscript received May 12, 1961). 


Fission-product cross sections have been 
measured radiochemically and mass-spectro- 
metrically for gold bombarded with 112-Mev C™” 
ions. Cross sections for 43 nuclides have been 
measured for elements from nickel to barium. 
Thirty-six yields are either primary fission- 
product yields (independent yields) or have been 
corrected (with less than 25% correction) so as 
to represent independent yields. The independent 
yields have been empirically systematized, and 
a yield- mass curve has been constructed. The 
yield-mass curve is compared with the yield- 
mass curves obtained from the fission of bismuth 
with 22-Mev and 190-Mev deuterons. The yield 
systematics indicate that the sum of the mass 
numbers of complementary fission products is 
13+1 amu less than that of the compound nucleus, 
and the sum of the charges of complementary fis- 
sion products is two units less than that of the 
compound nucleus. By thermodynamic arguments 
it is shown that the loss of charge was carried by 
an alpha particle, not by protons. 

The most probable charge of the fission prod- 
ucts as a function of mass number has been de- 
termined empirically and compared with theo- 
retical prediction. The charge-dispersion curve 
(fraction chain yield vs Z-Z,) may be fitted well 
by the Gaussian y =exp[-(Z 2 ,)*/0.9)/(0.9 n)¥?, 
Experimental yields on both sides of Z-Z, =0 
support the symmetry of the charge-dispersion 
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curve that many workers have assumed. 


*Present address: University of Rochester, Roch- 
ester, New York. 


CONVERSION OF MUONIUM INTO ANTIMUO- 
NIUM. G. Feinberg, Columbia University, New 
York, New York, and S. Weinberg, University 
of California, Berkeley, California (Received 
April 3, 1961). 


A detailed analysis is made of the possible 
conversion of muonium into antimuonium in 
various environments. An assumed flefie weak 
interaction of the usual form and strength gives 
a probability of 2.5x107° in vacuum, even in 
the presence of reasonable external electric 
fields. In a solid the probability is less by at 
least 10, and probably 20, orders of magnitude. 
In an inert gas the probability is roughly to be 
divided by the numbers of collisions during a 
muon lifetime, and hence is quite small unless 
the pressure at room temperature is less than 
about 10~* atmosphere. Lowering the tempera- 
ture does not help. A possible experiment is 
suggested. 


INELASTIC ELECTRON-DEUTERON SCATTER- 
ING CROSS SECTIONS AT HIGH ENERGIES. 

Il. FINAL-STATE INTERACTIONS AND RELA- 
TIVISTIC CORRECTIONS. Loyal Durand, -III, 
Brookhaven National Laboratory, Upton, New 
York (Received March 30, 1961; revised manu- 
script received May 11, 1961). 


Measurements of the cross section d*o/(dQ,dE ,') 
for the inelastic electron-deuteron scattering proc- 
ess e+d-—e+n+phave been used to determine the 
electromagnetic structure of the neutron. The ef- 
fects on the theoretical cross section of interac- 
tions between the outgoing nucleons are examined 
in detail using the methods of a previous paper. 
The transition matrix elements connecting the 
initial state of the two-nucleon system (the deu- 
teron) to a final state with specified total, orbital, 
and spin angular momenta are calculated using 
approximate wave functions which are matched 
to the experimentally determined neutron-proton 
scattering phase shifts. While individual matrix 
elements may be drastically changed by the dis- 
tortion of the final-state wave functions by the 
neutron-proton interaction, the over-all correc- 
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tions to the peak value of the cross section are 
found to be small (-1 to -2%) for electron mo- 
mentum transfers in the range q=3.4 -2.6 f. 
The precise magnitude of the corrections is 
somewhat uncertain because of the approximate 
nature of the wave functions, but it is unlikely 
either that they are large, or that the correc- 
tions could become positive. The effects of 
final-state interactions on the cross section 
d’q/(d2,dE,') are also examined for final elec- 
tron energies near the upper limit of the inelastic 
continuum. In this region, the nucleons emerge 
with low relative momenta, and, in agreement 
with the predictions of Jankus, the cross section 
is found to be drastically changed by the strong 
interactions in the final S states. However, it 
is shown that the presence in the neutron-proton 
interaction of a strongly repulsive core results 
in a considerable diminution of the cross section 
relative to the predictions of Jankus for large 
values of g. This lowering of the cross section 
has been observed by Kendall et al. Results 
obtained with approximate repulsive core wave 
functions provide a reasonable fit both to the 
inelastic cross section near the end point, and 
to the deuteron electromagnetic form factor 
obtained from elastic electron-deuteron scat- 
tering. Finally, the relativistic theory of ine- 
lastic electron-deuteron scattering is examined 
using the methods of dispersion relations. It is 
found that in the region of the large peak, the 
cross section d*0/(dQ,dE,') is given essentially 
correctly by a nonrelativistic calculation using 
a modified Hamiltonian, provided the results 
are interpreted correctly with respect to the 
kinematics. The approximations inherent in 
the calculation are examined in detail. The 
resulting cross section differs significantly 
from the modified Jankus cross section which 
has been used in the analysis of the high-energy 
electron-deuteron scattering data obtained by 
the Stanford group. It is found that the apparent 
values of the neutron charge form factor Fj, are 
reduced essentially to zero for q¢’ in the range 

5 f~? <q? < 20 f~* when relativistic correction, 
the effects of the deuteron D-state scattering, 
and the effects of final-state interactions are 
taken into account. Corresponding reductions 
in the value of the neutron anomalous magnetic 
moment form factor F9, range up to about 30%, 
and bring F5,, into closer agreement with Fo,. 
A complete reanalysis of the experimental data 
will be necessary. 


S-WAVE PION-NUCLEON SCATTERING. Jack 

L. Uretsky, Physics Department, Purdue Univer- 
sity, Lafayette, Indiana (Received November 17, 
1960; revised manuscript received May 18, 1961). 


The Mandelstam relations for pion-nucleon 
scattering are used to obtain equations for the 
S-wave partial wave amplitudes in the two isotop- 
ic spin states. The solutions of these equations 
are investigated in the approximation where only 
the one-nucleon contributions and the unitarity 
integral are kept. It is found that there are no 
solutions of the form N/D without complex zeros, 
and that this is a consequence of the large size 
of the one-nucleon terms. A comparison with 
experiment is made which suggests that the 
dominant contribution to the T=3/2 s-wave am- 
plitude (other than the one-nucleon contribution) 
comes from a region of the complex energy plane 
that is outside the physical region for the related 
processes (7-7 into NN and “crossed” m-N scat- 
tering). An Appendix is devoted to discussing 
the available experimental data and they are 
found to be consistent with a scattering length 
(6/k at threshold) of 0.098+ 0.004 in the T=3/2 
state. 


X-RAY DEFICIENCY IN MESONIC ATOMS. 

Y. Eisenberg, The Weizmann Institute of Sci- 
ence, Rehovoth, Israel, and D. Kessler, Israel 
Atomic Energy Establishment, Rehovoth, Israel 
(Received February 23, 1961; revised manuscript 
received May 9, 1961). 


We present an analysis of 7-mesonic atoms, 
based upon cascade calculations taking into ac- 
count the known processes of radiation, Auger 
transitions, and nuclear absorption. This anal- 
ysis, together with the previous one on p -mes- 
onic atoms, is intended to provide a deeper in- 
sight into the unsolved problem of the deficiency 
of x rays in mesonic atoms. It is shown that the 
m-mesonic L x-ray yields (for Z < 20) are quite 
insensitive to the strength of nuclear absorption 
and depend only upon the chosen initial meson 
population of the higher levels. Similarly, the 
ratios of basic (Kg, Lg, etc.) to higher x-ray 
yields, both for » and a mesons, depend strongly 
on the initial distribution. The best agreement 
between the calculations and experiment was ob- 
tained for a “modified statistical” initial popula - 
tion of the form (2/+1)e?!, with a=0.2, in the 
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n=14 level. From the existing experimental data 
on m-mesonic K x rays, the mean life of the 7 
meson in nuclear matter was deduced: To =2.75 
x10~*% sec. Within the framework of the present 
theory we are still unable to account for the x-ray 
deficiency in the light atoms. However, it is 
shown that the quantum loss as a function of en- 
ergy is different for 7- and .-mesonic.atoms, 
and therefore it is very probably due to a real 
physical effect. Furthermore, by comparing our 
predicted Auger electron yields with the experi- 
mental data, we can rule out any hypothetical 
simple Auger process in which the full energy 

of the “missing” quantum is given to a single 
electron. 


SECONDARY COSMIC-RAY PHOTONS BELOW 
CASCADE ENERGY. Kinsey A. Anderson, De- 
partment of Physics, University of California, 
Berkeley, California (Received April 7, 1961). 


Investigations with small unshielded scintilla- 
tion crystals carried through the atmosphere by 
balloons show large fluxes of photons in the en- 
ergy region 30 to 300 kev in equilibrium with the 
primary cosmic-ray beam. At 90 g cm~? depth 
the flux is about 22 photons cm~’ sec™ compared 
with a charged particle flux determined from a 
Geiger tube of 1.9 cm™sec™ at this same depth. 
The photon flux at zero depth, taken to be the 
albedo of this secondary cosmic-ray component, 
has been estimated by extrapolation to be 8 pho- 
tons cm’ sec™ greater than 30 kev. 


PHOTOPRODUCTION OF CHARGED 7 MESONS 
FROM NUCLEI. W. Melville McClelland,* Lab- 
oratory of Nuclear Studies, Cornell University, 
Ithaca, New York (Received March 6, 1961). 


The photoproduction of charged 7 mesons by a 
1000-Mev bremsstrahlung beam has been studied 
for the elements Be, C, Al, Cu, and Pb. Mesons 
with energies in the range 100 to 400 Mev emerg- 
ing from the targets at angles of 58° and 115° 
were detected, and absolute measurements for 
the cross section are given. An optical model for 
the nucleus was employed to predict absolute up- 
per and lower limits for the nuclear cross sec- 
tion, and reasonable agreement with the data was 
obtained. The measured cross sections had a de- 
pendence on the target atomic weight of A™ and 
this result lay between the limits predicted by 
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the model. The experimental nuclear 7~/7* ratio 
exhibited the general behavior of this quantity for 
deuterium, but the model could make no predic- 
tion here. The results seem to be consistent with 
an optical-model treatment of an assumed initial 
production of mesons throughout the nuclear vol- 
ume, and no recourse to a surface production 
mechanism was found to be necessary. 


*Present address: Lawrence Radiation Laboratory, 
University of California, Livermore, California. 


ANGULAR DISTRIBUTION OF PROTONS FROM 
1~-p SCATTERING AT 900 Mev. Bogdan C. 
Magli¢, Bernard T. Feld, and Carol A. Diffey, 
Physics Department and Laboratory for Nuclear 
Science, Massachusetts Institute of Technology, 
Cambridge, Massachusetts (Received April 10, 
1961). 


The shape of the 1~-p differential scattering 
cross section in the backward hemisphere should 
be sensitive to the nature of the “resonances” as- 
sumed to be responsible for the peaks in the total 
cross section at 600 and 900 Mev. The angular 
distribution of protons scattered in the forward 
hemisphere by pions of kinetic energy around 
925 Mev, corresponding to pion c.m. angles 
from 65 deg to 150 deg, was obtained by placing 
nuclear emulsions close to liquid hydrogen and 
by measuring the direction angle and the grain 
count of every proton track. It is shown that the 
sensitivity of emulsions in the temperature re- 
gion 22°K< 7 <90°K does not drop below 85% of 
the sensitivity at 300°K. The resulting distribu- 
tion is consistent with the assignment of D,,.“” 
and F,,.“*, respectively, for the 600- and 900- 
Mev levels. 


FORMATION OF N* IN HIGH-ENERGY NUCLEAR 
REACTIONS. I. Dostrovsky, The Weizmann In- 
stitute of Science, Rehovoth, Israel and Chemistry 
Department, Brookhaven National Laboratory, 
Upton, New York, Z. Fraenkel, The Weizmann 
Institute of Science, Rehovoth, Israel, and J. 
Hudis, Chemistry Department, Brookhaven Na- 
tional Laboratory, Upton, New York (Received 
April 10, 1961). 


Experimental cross sections are reported for 
the formation of N** in the bombardment of Zn, 
In, Pb, and U with protons of 0.9-, 1.9-, and 2.9- 
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Bev energy. These values are compared with 
theoretical N** emission cross sections for pro- 
ton energies of 0.84 and 1.84 Bev. The calcula- 
tions are based on the evaporation model. The 
previously described Monte Carlo procedure was 
modified in order to obtain better statistical ac- 
curacy for the calculated N** cross sections. 
Previously computed emission cross sections 

for He®, Li®, and Be’ were also recomputed using 
the modified Monte Carlo procedure. The cross 
sections were computed for three different formu- 
lations of the interaction radius. Good fit with 
the experimental He®, Li®, and Be’ cross sec- 
tions is obtained when the smaller values for the 
interaction radius are used. However, the fit 
with the experimental N** values is not good 
enough to exclude processes other than evapora- 
tion as contributing to the experimentally ob- 
served cross sections. 


PRODUCTION OF STRANGE PARTICLES IN p-p 
COLLISIONS AT 2.85 Bev. R. I. Louttit, T. W. 
Morris, D. C. Rahm,* R. R. Rau, A. M. Thorn- 
dike, and W. J. Willis, Brookhaven National Lab- 
oratory, Upton, New York, and R. M. Lea, The 
City College of New York, New York, New York 
(Received April 12, 1961). 


From a sample of 98 hyperon production events 
observed in a liquid hydrogen bubble chamber the 
partial cross sections for various final states are 
found to be: £*K*n—0.047, =*+K°p—0.030, D°K*p— 
0.013, A°K*tp—0.051, =“K*tpat—0.003, D*KNn—- 
0.004, (A°D°)Ktpn°—0.011, (A°D°)K°pn* —0.014, 
(A°D°)K*nn*+—0.002, all in millibarns. For the 
first four processes the values are in general 
agreement with those calculated by Ferrari using 
a one-pion-exchange model. Only one exampie of 
K-pair production was observed, indicating a 
cross section less than 0.01 millibarn. 


*Now on leave at Centre d’ Etudes Nucléaires de 
Saclay, Saclay, France. 


EVOLUTION OF A QUASI-STATIONARY STATE. 
Rolf G. Winter, Department of Physics, The 
Pennsylvania State University, University Park, 
Pennsylvania (Received April 10, 1961). 


To elucidate the time development of quasi- 
stationary states, a simple barrier penetration 
problem has been studied. Both approximate 
expressions and numerical results for some 


parameters were obtained for the decay rate. 
First, irregular oscillations occur for a short 
time. Second, the exponential region follows. 
Third, further oscillations occur during which 
the decay rate dips to negative values, so that 
the probability of finding the undecayed system 
increases briefly at several times. Fourth and 
finally, the decay rate decreases like an inverse 
power of the time. 


DYNAMICAL THEORY FOR STRONG INTERAC- 
TIONS AT LOW MOMENTUM TRANSFERS BUT 
ARBITRARY ENERGIES. Geoffrey F. Chew and 
Steven C. Frautschi, Lawrence Radiation Labora- 
tory and Department of Physics, University of 
California, Berkeley, California (Received March 
30, 1961). 


Starting from the Mandelstam representation, 
it is argued on physical grounds that “strips” 
along the boundaries of the double spectral re- 
gions are likely to control the physical elastic 
scattering amplitude for arbitrarily high energies 
at small momentum transfers. Pion-pion scat- 
tering is used as an illustration to show how the 
double spectral functions in the nearest strip re- 
gions may be calculated, and an attempt is made 
to formulate an approximate but “complete” set 
of dynamical equations. The asymptotic behavior 
of the solutions of these equations is discussed, 
and it is shown that if the total cross section is 
to approach a constant at large energies then at 
low energy the S-dominant mz solution is inad- 
missible. A principle of “maximum strength” 
for strong interactions is proposed, and it is 
argued that such a principle will allow large 
low-energy phase shifts only for / </,,,,, where 


es 


THEORY OF 2-N SCATTERING IN THE STRIP 
APPROXIMATION TO THE MANDELSTAM REP- 
RESENTATION. Virendra Singh* and B. M. Ud- 
gaonkar,! Lawrence Radiation Laboratory, Uni- 
versity of California, Berkeley, California (Re- 
ceived April 7, 1961). 


The strip approximation to the Mandelstam 
representation is applied to the 7-N problem, 
and the basic equations given. The asymptotic 
behavior of the invariant amplitudes in the physi- 
cal regions is discussed in terms of the unitarity 
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condition on partial-wave amplitudes, the con- 
stancy of high-energy scattering cross sections, 
and the Pomeranchuk theorem, and it is shown 
to imply that no subtractions should be necessary 
except in the J = 1/2 wave of the 7-N channel and 
the J =0 wave of the 7+1~N+N channel. This 
obviates the difficulties encountered by earlier 
workers when they subtracted higher waves. 


*On deputation from Tata Institute of Fundamental 
Research, Bombay, India. 


tOn deputation from Atomic Energy Establishment, 
Trombay, India. 


THEORY OF STRONG-INTERACTION SYMME- 
TRY. Murray Gell-Mann, California Institute 
of Technology, Pasadena, California (Received 
March 27, 1961). 


We attempt once more, as in the global sym- 
metry scheme, to treat the eight known baryons 
as a supermultiplet, degenerate in the limit of 
a certain symmetry but split into isotopic spin 
multiplets by a symmetry-breaking term. Here 
we do not try to describe the symmetry violation 
in detail, but we ascribe it phenomenologically to 
the mass differences themselves, supposing that 


there is some analogy to the p-e mass difference. 


The symmetry is called unitary symmetry and 
corresponds to the “unitary group” in three di- 
mensions in the same way that charge independ - 
ence corresponds to the “unitary group” in two 
dimensions. The eight infinitesimal generators 
of the group form a simple Lie algebra, just like 
the three components of isotopic spin. In this 
important sense, unitary symmetry is the sim- 
plest generalization of charge independence. The 
baryons then correspond naturally to an eight- 
dimensional irreducible representation of the 
group; when the mass differences are turned on, 
the familiar multiplets appear. The pion and K 
meson fit into a similar set of eight particles, 
along with a predicted pseudoscalar meson ° 
having J=0. The pattern of Yukawa couplings of 
na, K, and yx is then nearly determined, in the 
limit of unitary symmetry. 

The most attractive feature of the scheme is 
that it permits the description of eight vector 
mesons by a unified theory of the Yang- Mills 
type (with a mass term). Like Sakurai, we have 
a triplet p of vector mesons coupled to the isotop- 
ic spin current and a singlet vector meson w® 
coupled to the hypercharge current. We also 
have a pair of doublets M and M, strange vector 
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mesons coupled to strangeness-changing cur- 
rents that are conserved when the mass dif- 
ferences are turned off. There is only one coup- 
ling constant, in the symmetric limit, for the 
system of eight vector mesons. There is some 
experimental evidence for the existence of w°® 
and M, while p is presumably the famous /=1, 
J=1, m-7 resonance. 

A ninth vector meson coupled to the baryon 
current can be accommodated naturally in the 
scheme. 

The most important prediction is the qualitative 
one that the eight baryons should all have the 
same spin and parity and that the pseudoscalar 
and vector mesons should form “octets,” with 
possible additional “singlets.” 

If the symmetry is not too badly broken in the 
case of the renormalized coupling constants of 
the eight vector mesons, then numerous detailed 
predictions can be made of experimental results. 

The mathematics of the unitary group is de- 
scribed by considering three fictitious “leptons,” 
v, e~, and y~, which may or may not have some- 
thing to do with real leptons. If there is a con- 
nection, then it may throw light on the structure 
of the weak interactions. 


HYPERVIRIAL THEOREMS FOR VARIATIONAL 
WAVE FUNCTIONS. Saul T. Epstein, Department 
of Physics, University of Nebraska, Lincoln, 
Nebraska, and Joseph O. Hirschfelder, Theoret- 
ical Chemistry, University of Wisconsin, Madi- 
son, Wisconsin (Received March 24, 1961). 


It is shown that a sufficient condition for an op- 
timal energy variational wave function y, to sat- 
isfy the hypervirial relation (¥,,[H, W]y,) =0 is 
for the trial function y to admit variations of the 
form 8y/8a=(i/f)Wy. Here H is the Hamiltonian, 
W is a Hermitian operator, and a is a variational 
parameter. Explicit forms of such trial functions 
are exhibited for several W’s. The case in which 
W generates a point transformation of the coor- 
dinates is discussed in detail. Conditions are 
given for the existence of simultaneous hyper- 
virial theorems. 


FURTHER CONSIDERATIONS ON ELECTROMAG- 


NETIC POTENTIALS IN THE QUANTUM THEORY. 


Y. Aharonov, Department of Physics, Brandeis 
University, Waltham, Massachusetts, and D. 
Bohm, H. H. Wills Physics Laboratory, Univer- 
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sity of Bristol, Bristol, England (Received April 
6, 1961). 


In this article we discuss in further detail the 
significance of potentials in the quantum theory, 
and in so doing, we answer a number of argu- 
ments that have been raised against the conclu- 
sions of our first paper on the same subject. We 
then proceed to extend our treatment to include 
the sources of potentials quantum mechanically, 
and we show that when this is done, the same 
results are obtained as those of our first paper, 
in which the potential was taken to be a specified 
function of space and time. In this way, we not 
only answer certain additional criticisms that 
have been made of the original treatment, but 
we also bring out more clearly the importance 
of the potential in the expression of the local 
character of the interaction of charged particles 
and the electromagnetic field. 


TRANSFER OF HELICITY IN RADIATION AND 
ABSORPTION OF HIGH-ENERGY PHOTONS. 
R. H. Pratt, Institute of Theoretical Physics, 
Department of Physics, Stanford University, 
Stanford, California (Received April 7, 1961). 


Nearly complete transfer of momentum be- 
tween a high-energy electron (or positron) and 
a photon in a Coulomb field implies that helicity 
is also transferred. This is not a consequence 
of conservation of total angular momentum but, 
rather, of spin angular momentum, and follows 
from a demonstration that it is possible to use 
free-particle spinors (though not free wave 
functions) for the high-energy particles. Polar- 
ization correlations of the lower energy particle 
in such a process are discussed. Applications 
are made to bremsstrahlung, pair production, 
photoeffect, and one-photon pair annihilation. 





ANNOUNCEMENTS 


Supplement 
Beginning with the issue of July 1, the “Abstracts of Articles 
to be published in The Physical Review” will be printed as a loose 
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